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Executive  Summary 


This  report  documents  the  latest  version  of  the  Stochastic  Lanchester  Air-to-Air 
Campaign  Model  (SLAACM),  developed  by  LMI  for  the  Tactical  Air  Forces  Di¬ 
vision  of  the  Program  Analysis  and  Evaluation  Directorate  of  the  Office  of  the 
Secretary  of  Defense  (OSD  PA&E  TACAIR). 

SLAACM  is  a  fast  and  flexible  model  designed  for  analysts  who  need  to  consider 
many  combat  scenario  options  quickly  and  who  wish  to  have  indications  of  the 
uncertainties  in  military  outcomes.  SLAACM  models  the  defensive  counter-air 
battle,  including  order  of  battle  optimization  by  both  Red  attackers  and  Blue  de¬ 
fenders.  The  current  version  can  address  issues  of  fighter  combat  effectiveness, 
Blue-side  battle  management  capabilities,  time -phasing  of  offensive  and  defen¬ 
sive  forces,  use  of  indigenous  defensive  forces,  and  basing  options  for  deployed 
defensive  forces.  Unique,  efficient  mathematical  algorithms  and  fast  integer  pro¬ 
gramming  tools  allow  SLAACM  to  generate  optimized  results  including  statisti¬ 
cal  uncertainties  for  numerous  types  and  large  numbers  of  combatants  for 
campaigns  of  many  days,  with  sufficient  speed  to  allow  statistical  variation  of  in¬ 
put  parameters.  Included  tables  and  charts  display  the  day-by-day  development  of 
the  campaign.  In  addition,  the  Microsoft  Excel  implementation  of  the  model  al¬ 
lows  wide  flexibility  for  storing  and  displaying  results  in  user-preferred  formats. 

This  report  includes  descriptions  of  the  fundamental  mathematics  and  analytical 
structure  of  SLAACM  and  a  user’s  guide.  The  report  also  includes  discussions  of 
SLAACM  extensions  to  address  issues  such  as  weapons  expenditure  (missile  count¬ 
ing)  and  engagements  involving  radar  tracking  and  electronic  countermeasures. 

Planned  future  work  includes  application  of  SLAACM  methodology  to  optimized 
offensive  air  campaigns,  including  suppression  of  enemy  air  defense  and  the  im¬ 
pact  of  electronic  warfare  and  low  observable  technology. 

The  authors  wish  to  thank  the  Director  of  OSD  PA&E,  Mr.  Bradley  Berkson  for 
his  support  of  the  development  of  SLAACM. 
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Chapter  1 

Introduction 


The  Stochastic  Lanchester  Air- to- Air  Campaign  Model  (SLAACM),  developed 
by  LMI,  is  a  parsimonious  probabilistic  defensive  counter-air  model  intended  to 
be  a  flexible  tool  for  analysts  who  need  to  consider  many  cases  quickly  in  the  PC 
environment  and  who  wish  to  have  indications  of  the  uncertainties  in  military  out¬ 
comes.  SLAACM  outputs  in  benchmarking  tests  compare  closely  with  those  of 
some  large-scale  simulation  models,  but  SLAACM  is  intended  to  complement, 
rather  than  replace,  such  models.  SLAACM  typically  runs  a  fully  optimized  10- 
day  campaign  on  a  standard  PC  in  less  than  2  minutes.  This  report  describes  the 
analytical  and  statistical  methods  used  to  generate  SLAACM  outputs.  It  also  gives 
detailed  instructions  for  operating  SLAACM,  with  illustrative  examples. 

SLAACM  Overview 


SLAACM  treats  campaigns  between  two  opponents,  called  Red  and  Blue.  The  Red 
side  uses  fighter-escorted  bombers  to  attack  assets  defended  by  the  Blue  side’s 
lighters.  In  each  day’s  operations,  the  Red  side  detennines  optimal  allocations  of 
its  fighters  and  bombers  to  attack  packages,  maximizing  a  payoff  function  that  con¬ 
siders  the  bombers’  payload  and  weapon-delivery  accuracy,  the  value  to  Red  of 
destroying  Blue  fighters,  and  the  penalty  to  Red  of  losing  fighters  and  bombers. 

Blue’s  fighters  respond.  Each  Blue  fighter  can  be  “smart”  or  not.  Smart  fighters 
can  determine  the  makeup  of  oncoming  Red  attack  packages  before  engaging  Red 
and  can  coordinate  their  operations  to  make  optimal  defenses,  maximizing  a  pay¬ 
off  function  that  considers  the  value  to  Blue  of  destroying  Red  fighters  and  bomb¬ 
ers,  and  the  penalty  to  Blue  of  losing  a  fighter. 

Blue  fighters  that  are  not  smart  encounter  Red’s  packages  randomly.  An  adjust¬ 
able  parameter  varies  the  effectiveness  of  Blue’s  battle  management  for  not-smart 
fighters,  reducing  the  probability  of  engagements  below  that  implied  simply  by 
the  numbers  of  attacking  packages  and  defending  flights. 

SLAACM  does  not  use  iterative  simulation.  Rather,  it  calculates  the  outcome 
probabilities  of  Blue  fighters  versus  Red  fighter-plus-bomber  packages  using  ana¬ 
lytic  probability  algorithms.  With  those  results,  SLAACM  determines  statistics  of 
the  results  of  a  day’s  combat.  These  are  used  to  determine  Red  and  Blue  orders  of 
battle  for  the  next  day’s  operations. 1 


1  The  term  “day”  here  refers  to  a  single  set  of  campaign  sorties.  There  could  be  more  than  one 
such  set  on  a  single  calendar  day  or  there  could  be  calendar  day  gaps  between  sets. 
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SLAACM  tracks  the  total  tons  of  bombs  delivered,  as  well  as  the  tons  of  smart 
bombs  (guided  munitions  with  higher  destructive  capability  per  ton).  SLAACM 
also  tracks  losses  to  both  sides,  which  Blue  aircraft  killed  which  Red  aircraft  and 
vice  versa,  and  standard  deviations  of  Blue  and  Red  losses. 

SLAACM ’s  worksheets  and  analytic  structure  provide  complete  audit  trails  for 
the  input  parameters  used  in  the  calculations.  The  kill  rate  ratios  used  in  the  calcu¬ 
lations  are  derived  directly  from  input  loss  ratios  that  are  displayed  in  the  model. 
In  classified  versions  of  SLAACM,  the  loss  ratios  are  based  on  documented  his¬ 
torical  data,  simulations,  and  military  expertise.  Further,  the  analytic  nature  of  the 
model  ensures  that  the  effects  of  parameter  changes  are  consistently  and  individu¬ 
ally  reflected  in  the  output  and  are  not  convolved  with  other  parameter  effects  as 
happens  in  iterative  simulations.  These  features  can  be  quite  helpful  to  clarify  the 
implications  of  the  model’s  results. 

Report  Organization 


This  report  is  organized  as  follows: 

♦  Chapters  2  through  5  describe  how  SLAACM  works.  Specifically,  those 
chapters,  respectively,  describe  how  the  model  is  used  to 

>  calculate  air-to-air  engagement  results, 

>  calculate  bomber  effectiveness, 

>  specify  the  Red  and  Blue  payoff  logic  and  calculate  optimal  attacks 
and  defenses,  and 

>  calculate  campaign  results. 

♦  Chapter  6  is  a  users  manual  for  operating  SLAACM. 

The  appendixes  contain  supplemental  information  of  interest  for  combat  modeling: 

♦  Appendix  A  discusses  the  mathematics  of  the  fast  engagement  analysis  al¬ 
gorithm. 

♦  Appendix  B  addresses  the  use  of  NASA-developed  Markov  tools  for  aux¬ 
iliary  engagement  analysis  and  scenario  development. 

♦  Appendix  C  discusses  SLAACM  alternatives. 

♦  Appendix  D  discusses  advanced  probabilistic  engagement  models. 

♦  Appendix  E  discusses  modeling  of  electronic  attack 
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♦  Appendix  F  discusses  a  mathematical  analysis  of  phased  attack  and  de¬ 
fense 
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Chapter  2 

How  SLAACM  Works: 

Engagement  Models  and  Calculations 


SLAACM  does  not  simulate;  rather,  it  makes  complete  calculations  of  outcome 
probability  distributions  for  specific  probabilistic  engagement  models.  This  is  a  sig¬ 
nificant  difference  between  SLAACM  and  other  campaign  models,  so  it  is  appropri¬ 
ate  to  begin  our  discussion  of  how  SLAACM  works  by  explaining  that  approach. 

The  Basic  M  vs.  N  Probabilistic 
Engagement  Model 


The  heart  of  SLAACM’s  calculations  is  a  set  of  probabilistic  models  of  air-to-air 
engagements.  The  simplest  engagement  model  used  in  SLAACM  is  the  M  vs.  N 
probabilistic  engagement  model.  We  describe  that  model  in  considerable  detail,  to 
clarify  some  fundamental  aspects  of  SLAACM’s  operation. 

In  the  basic  M  vs.  N  probabilistic  model,  M  Blue  aircraft  engage  N  Red  aircraft. 
Throughout  the  engagement,  the  time  between  kills  by  each  Blue  aircraft  is  taken 
to  be  identically  and  independently  distributed  exponentially  with  parameter  kb, 
and  the  time  between  kills  by  each  Red  aircraft  is  assumed  to  be  identically  and 
independently  distributed  exponentially  with  parameter  kr.  Specifically, 

xr  ~  kre-krTr  [Eq.  2-1] 

and 

ib~kbe-kbI\  [Eq.  2-2] 

where  the  means  “distributed  as”  and  xr  and  ib  are,  respectively,  the  time  be¬ 
tween  kills  made  by  each  Red  aircraft  and  the  time  between  kills  made  by  each 
Blue  aircraft. 

If  the  time  between  kills  by  an  engaged  aircraft  is  exponentially  distributed  with 
parameter  k,  then  the  probability  that  the  aircraft  makes  a  kill  in  a  time  interval 
(t,  t  +  5t)  is  k5t  +  0(8t  ),  independent  of  t.  (This  property  is  sometimes  called  the 
“memoryless”  property  of  exponential  distributions.) 


1  0(5t2)  signifies  “higher  order  terms”  of  order  5t2  and  higher,  which  are  very  small  and  can  be 
ignored. 
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If  two  such  aircraft  are  present,  then  by  virtue  of  the  assumption  that  their  times 
between  kills  are  statistically  independent,  the  probability  that  the  pair  makes  a 
kill  in  the  interval  (t,  t  +  5t)  is  k5t  +  k5t  +  0(8t2)  =  2k8t  +  0(8t2).  If  m  such  aircraft 
are  present,  the  probability  that  the  set  of  them  makes  a  kill  in  the  interval 
(t,  t  +  8t)  is  mk8t  +  0(8t2). 

If,  at  a  given  time  during  the  engagement,  m  Blue  aircraft  and  n  Red  aircraft  are 
present,  we  will  say  that  the  engagement  is  in  state  (m,  n).  Figure  2-1  diagrams 
the  possible  transitions  that  the  engagement  can  make  into  and  out  of  state  (m,  n). 


Figure  2-1.  Engagement  State  Transition  Diagram 


Now  let  P m,n(t)  be  the  probability  that  the  engagement  is  in  state  (m,  n)  at  time  t. 
Let  us  consider  the  probability  Pm,n(t  +  8t).  The  engagement  can  reach  state  (m,  n) 
at  time  (t  +  8t)  in  just  three  ways: 

♦  The  engagement  was  in  state  (m,  n)  at  time  t,  and  no  kills  happened  dur¬ 
ing  the  interval  (t,  t  +  St). 

♦  The  engagement  was  in  state  (m,  n  +  1)  at  time  t,  and  the  set  of  Blue  air¬ 
craft  made  a  kill  during  the  interval  (t,  t  +  St). 

♦  The  engagement  was  in  state  (m  +  1,  n)  at  time  t,  and  the  set  of  Red  air¬ 
craft  made  a  kill  during  the  interval  (t,  t  +  St). 

By  the  assumptions  that  all  the  kill  events  are  independent,  these  three  events  are 
statistically  independent.  Accordingly,  the  probability  that  one  of  them  occurs  is 
the  sum  of  their  individual  probabilities,  so  that,  neglecting  terms  of  0(5t2), 

pm,n  (t  +  St)  =  Pm  n  (t)[l  -  mk  b8t  -  nk  r8t]+  Pm  n+1  (t)mk  b8t  +  Pm+I  n  (t)nk  r8t  .[Eq.  2-3] 
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How  SLAACM  Works:  Engagement  Models  and  Calculations 


Subtracting  Pm,n(t)  from  each  side  of  Equation  2-3,  dividing  by  8t,  and  taking  the 
limit  as  St  — »  0  gives 

Pm,n (t)  =  - (mkb  +  nkr  )Pm  n  +  mkbPm  n+1  +  nkrPm+u .  [Eq.  2-4] 


With  the  result  (Equation  2-4),  we  can  write  down  a  general  initial  value  problem 
describing  the  evolution  of  the  statistics  of  an  M  vs.  N  engagement,  in  which  the 
Blue  forces  break  away  when  they  have  fewer  than  bmin  aircraft,  and  the  Red 
forces  break  away  when  they  have  fewer  than  rmjn  aircraft: 


^m,n  ^^Pr)Pm,n  ~*~^kbPm  n+j  +TlkrPm+j  n/ 

vbmin  <  m  <  M  and  rmin  <  n  <  N; 
po,n  —  r|krPbmjnn  ,  rmin  <  n  <  N 
Pm,0  —  ^kbPm,rmin'  bmin  —  ni  <  M 

pm,n  (0)  =  0  except  PM/N  (0)  =  1 


[Eq.  2-5] 


(The  probabilities  that  the  system  has  more  than  M  Blues,  or  more  than  N  Reds, 
are  of  course  always  zero.) 

The  initial  value  problem  (Equation  2-5)  presents  a  system  of  linear  ordinary  dif¬ 
ferential  equations.  The  number  of  equations  is  equal  to  (M  -  bmin  +  1)(N  -  rmjn  + 
1)  +  (M  -  bmin  +1)  +  (N  -  rmin  +1).  The  equations  for  the  probabilities  of  the  ab¬ 
sorbing  boundary  states,  Po,n  and  Pm,o,  decouple  from  those  for  the  transient  states 
(those  whose  probabilities  tend  to  zero  for  long  times),  however.  Therefore,  Equa¬ 
tion  2-5  may  be  treated  by  solving  only  the  (M  -  bmin  +  1)(N  -  rmjn  +  1)  equations 
for  the  transient  states  and  computing  the  probabilities  of  the  absorbing  boundary 
states  by  integration. 

In  any  case,  Equation  2-5  offers  no  difficulties  to  numerical  solution.  Moreover, 
introducing  the  non-dimensional  time  parameter  x  =  krt  allows  us  to  write  Equa¬ 
tion  2-5  in  a  form  that  involves  only  the  non-dimensional  parameter  Kill  Rate  Ra¬ 
tio  =  KRR  =  k  =  kb/kr: 


Pm,n  =  -(mK  +  n)pm.n  +  mKPm,„+l  +  nPm+l,„> 

vb,„m  <  rn  <  M  and  rmin  <  n  <  N; 

P0,n  =  nPb„„n,n  »  fn.n  <  n  <  N 
Pm,0  =  mKPm,rmi„>  b,n,n  <™<M 
Pm,n(°)  =  0  eXCePt  PM,N(°)  =  1 


[Eq.  2-6] 


The  state  probabilities  in  Equation  2-6  are  functions  of  x,  not  t,  and  so  are  differ¬ 
ent  functions  from  those  giving  the  state  probabilities  as  functions  of  t.  Also  the 
derivatives  in  Equation  2-6  are  with  respect  to  x,  not  t.  In  the  interest  of  unclut- 
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tered  notation,  we  will  not  introduce  new  symbols  for  these  new  functions;  rather, 
in  Equation  2-6  and  in  all  following  parts  of  this  chapter,  the  probabilities  are 
functions  of  x,  and  dots  denote  differentiation  with  respect  to  x. 

Writing  the  initial  value  problem  for  the  M  vs.  N  probabilistic  engagement  in  non- 
dimensional  form  gives  useful  insights  into  the  nature  of  the  problem’s  solutions. 
First,  because  the  only  parameter  appearing  in  Equation  2-6  is  k,  we  see  that  the 
probabilities  Pm  n.  which  in  general  depend  on  time  t  and  the  two  parameters  kr 
and  kb,  can  always  be  written  as  functions  of  the  non-dimensional  time  x  and  just 
one  parameter,  the  kill-rate  ratio  k. 

Second,  the  fact  that  the  solutions  of  the  initial  value  problem  depend  only  on  k 
and  x  leads  to  the  conclusion  that  the  limiting  values  of  the  probabilities  as  time 
tends  to  infinity  (which  of  course  is  also  the  limit  as  x  tends  to  infinity)  are  func¬ 
tions  only  of  the  single  parameter  k. 

Figure  2-2  illustrates  these  ideas.  It  shows  the  evolution  of  selected  probabilities, 
obtained  by  numerical  integration  for  a  4  vs.  4  engagement  with  kb/kr  =  k  =  1.46, 
and  in  which  the  combatants  fight  to  annihilation.2 


0  0.5  1  1.5  2 

kr  t 


— • —  P4,4  — P0,3  -A— P0,4  - P3,0  — P4,0 


To  find  the  value  of  one  of  the  probabilities  shown  in  Figure  2-2  at  a  given  time  t 
for  a  given  engagement,  one  would  find  the  value  of  x  corresponding  to  that  t  by 
computing  x  =  krt  and  reading  the  desired  value  from  the  figure.  For  example,  at 
t  =  20  time  units  in  an  engagement  in  which  the  mean  time  for  the  Red  aircraft  to 
make  a  kill  is  180  time  units,  the  value  of  x  is  20/180  ~  0.1 1;  thus  the  desired 
value  ofP4j4  is  0.4. 

2  This  engagement  also  includes  absorbing  states  (2,  0),  (1,  0),  (0,  2),  and  (0,  1)  that  are  not 
shown  in  the  figure. 
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How  SLAACM  Works:  Engagement  Models  and  Calculations 


We  will  be  more  interested  in  the  long-time  limiting  values  of  the  probabilities. 
Figure  2-2  indicates  that  for  x  >  2,  the  absorbing  boundary  states  have  essentially 
reached  their  long-time  limiting  values.  The  numerical  solution  confirms  that  for 
x  >  2,  the  total  probability  of  the  absorbing  boundary  states  is  larger  than  95  per¬ 
cent.  The  physical  time  required  to  reach  the  long-time  limit  will  of  course  vary 
with  the  value  of  the  dimensional  parameter  kr.  But  the  values  of  the  probabilities 
depend  only  on  the  single,  non-dimensional  parameter  k. 

Relating  Loss  Data  to  Kill  Rate  Ratios 

The  engagement  models  underlying  SLAACM,  like  all  Lanchester-type  models, 
employ  kill  rates  for  the  opposing  forces.  In  detenninistic  Lanchester  models,  the 
parameters  are  kills  per  (i rer  per  unit  time.  In  stochastic  Lanchester  models  like 
SLAACM,  the  parameters  are  the  reciprocals  of  mean  times  between  kills  by  a 
single  H  rer.  For  simplicity,  we  refer  to  parameters  of  stochastic  Lanchester  mod¬ 
els  as  “kill  rates.” 

Observed  engagement  data,  from  simulations  or  from  observations  of  actual  com¬ 
bat,  are  rarely,  if  ever,  rate  data.  Rather,  they  are  typically  loss  ratios  for  specific 
M  vs.  N  engagements.  To  apply  empirical  loss  ratio  data  to  rate-based  models,  we 
need  a  method  to  convert  from  loss  ratio  data  to  kill  rate  data.  To  do  this,  we 
make  the  significant  assumption  that  the  long-time  limits  of  SLAACM’s  engage¬ 
ment  models  are  appropriate  models  for  the  observations. 

Thus,  with  the  basic  M  vs.  N  engagement  model,  we  assume  that  the  observed 
engagements  end  only  with  one  side’s  reaching  its  specific  breakaway  condi¬ 
tion — a  certain  number  of  losses,  or  annihilation — rather  than  with  a  clock-time- 
dependent  condition  such  as  fuel  exhaustion,  or  an  event-dependent  condition 
such  as  exhaustion  of  munitions.  (More  advanced  engagement  models,  described 
in  Appendixes  B  and  D,  do  account  for  missiles,  and  can  be  used  in  SLAACM.) 

As  noted  above,  the  non-dimensional  kill  rate  ratio  parameter,  k,  is  the  single  pa¬ 
rameter  of  interest  in  the  basic  engagement  model.  To  infer  a  value  of  k  from  loss 
ratio  data  for  an  M  vs.  N  engagement,  we  find  the  value  of  k  for  which  the  long¬ 
time  limiting  value  of  the  ratio  of  expected  Blue  loss  to  expected  Red  loss  in  the 
basic  M  vs.  N  probabilistic  model  is  equal  to  the  observed  loss  ratio.  In  this  way, 
we  infer  a  non-dimensional  quantity,  k,  from  a  non-dimensional  quantity,  the  ob¬ 
served  loss  ratio. 

Suppose,  for  example,  that  the  available  data  are  for  a  2  vs.  2  engagement  to  anni¬ 
hilation  and  that  they  show  a  loss  ratio  of  1  Blue  to  10  Reds.  Calculations  show 
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that,  for  the  basic  probabilistic  engagement  model,  a  kill  rate  ratio  k  of  8.28 

-5 

makes  the  ratio  of  expected  Blue  loss  to  expected  Red  loss  equal  to  1:10. 

Although  we  usually  use  a  finite-step  iterative  method  to  compute  long-time  lim¬ 
its  directly,  it  is  of  course  possible  to  find  long-time  limits  with  time -based  meth¬ 
ods,  such  as  numerical  integration.  In  time-based  Markov  calculations,  the  time 
required  to  reach  the  long-time  limit  will  depend  on  the  value  chosen  for  kr,  where 
kb  is  equal  to  k  times  kr.  For  example,  analysis  shows  that  the  kill  rate  ratio  of 
8.28  results  in  a  1:10  loss  ratio  in  approximately  0.22  time  units  for  kr=  1,  in 
roughly  2.2  time  units  iterations  for  kr=  0. 1,  and  in  about  22  time  units  for 
kr  =  0.01.  The  0.2,  2,  and  22  time  units  represent  adequate  model  time  to  bring  the 
solution  to  the  long-time  limit,  for  their  corresponding  kill  rates.  Thus,  in  a  time- 
unit-based  model,  selection  of  the  individual  rates  is  arbitrary  as  long  as  the  ratio 
k  is  preserved  and  the  calculations  continue  to  a  time  on  the  order  of  twice  the  in¬ 
verse  of  the  sum  of  the  parameters  kb  and  kr.  (Appendix  A  describes  the  general 
analytic  method  that  produces  these  numerical  results.) 

Another  major  assumption  in  our  engagement  modeling  is  that  the  kill  rate  ratio  is 
independent  of  the  Blue:Red  aircraft  ratio.  We  assume  that  kill  rate  ratios  derived 
from  4  vs.  8  data  can  be  applied  directly  to  4  vs.  4  and  other  engagement  combi¬ 
nations.  A  corollary  assumption  is  that  kill  rate  ratios  are  independent  of  the  state 
of  an  engagement.  These  are  significant  and  powerful  simplifying  assumptions, 
which  we  are  always  eager  to  test  with  new  combat  data. 

Available  data  from  TAC  BRAWLER  engagement  simulations  give  expected  val¬ 
ues  of  Red  losses  and  Blue  losses.  These  data,  some  of  them  classified,  are  avail¬ 
able  for  several  Red  and  Blue  aircraft  types  of  interest.  Government  subject 
matter  experts  have  scaled  these  results  to  model  additional  aircraft  combinations. 
A  SLAACM  utility  uses  bisection  iteration  and  the  long-time  limit  algorithms 
(discussed  below  and  in  Appendix  A)  to  calculate  the  values  of  k  from  the  loss 
ratios  for  each  of  the  Red  vs.  Blue  combinations.  The  calculations  are  repeated 
until  successive  iterations  agree  to  four  significant  figures. 

Calculations  for  the  Long-Time  Limit 

If  the  engagement  continues  to  completion — that  is,  until  the  probabilities  that  the 
system  is  not  in  one  of  the  absorbing  boundary  states  (0,  n)  or  (m,  0)  are  negligi¬ 
bly  small — then  the  probabilities  of  these  states  are  functions  only  of  k  and  the 
initial  values  M  and  N.  These  conditions  allow  significant  simplification  in  the 
calculation  of  state  probabilities.  One  can  obtain  values  of  the  long-time  limits  of 
the  probabilities  of  absorbing  boundary  states  by  a  simple,  finite-step  iteration. 


3  The  kill  rate  ratio  reflects  the  performance  of  individual  Blue  vs.  Red  pairs  and  is  numeri¬ 
cally  equal  to  the  loss  ratio  for  a  1  vs.  1  engagement.  A  10:1  kill  rate  ratio  would  be  required  to 
produce  a  1:10  loss  ratio  for  a  1  vs.  1  engagement,  but  only  an  8.28:1  kill  rate  ratio  is  necessary  for 
the  same  loss  ratio  in  a  2  vs.  2  engagement,  because  all  available  aircraft  are  contributing  to  each 
kill. 
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This  speeds  up  the  calculations  tremendously,  compared  to  time-based  Markov 
methods,  and  is  a  key  to  SLAACM  utility. 

An  Example 

As  an  example  of  the  use  of  the  initial  value  problem  (Equation  2-5),  we  consider 
four  Blues  engaging  four  Reds.  Both  sides  break  away  on  sustaining  2  losses,  that 
is?  bmin  rmm  3. 

Figure  2-3  illustrates  the  states  in  the  engagement  and  the  transitions  between 
them.  In  that  figure,  the  absorbing  boundary  states  are  labeled  as  having  0  aircraft 
for  the  side  that  has  broken  away. 

Figure  2-3.  Engagement  Diagram 


Red  Win  States 


3kr 


Blue  Win  States 


3,0  <*- 


3kb 


3kb 

3,3  * - ( 3,4 


3kr 
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4kb 


4kr 


4,3  4»4,4 


For  this  case,  the  general  initial  boundary  value  problem  (Equation  2-6)  take  this 
specific  form: 


P4/4  =  -(4k  +  4)P4  4 
P4  3  =-(4k  +  3)P4  3  +  4kP4  4 
1*3,4  =  -(3k  +  4)P4  3  +  4P4  4 
P},3  =  _(3k  +  3)P3  3  +  3P4  3  +  3kP3  4 

P4,0  =  4kP4,3  >  [Eq.  2-7] 

P3,0  =  3^3,3 
P<),3  =  3P33 
P0,4  =  4P3,4 

pm,n  (0)  =  0,  except  P4  4  (0)  =  1 
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where  k  =  kb/kr  and  the  (m,  0)  and  (0,  n)  states  correspond  to  the  break  condition 
of  2  aircraft. 

The  Red  side  wins  the  engagement  if  all  the  Blues  are  gone,  that  is,  if  the  en¬ 
gagement  is  in  a  state  (0,  n),  rrain  <  n  <  N.  The  Blue  side  wins  if  the  engagement  is 
in  a  state  (m,  0),  bmin  <  m  <  M. 

Note  that  the  first  four  equations  in  Equation  2-7  describe  the  evolution  of  the 
probabilities  of  transient  states,  whose  probabilities  will  tend  to  zero  at  long  times. 
The  remaining  four  equations  describe  the  variation  of  the  absorbing  boundary 
states,  or  outcome  states.  At  long  times,  the  system  will  be  found  in  one  of  these 
states. 

The  equations  for  the  transient  states  decouple  from  those  of  the  outcome  states. 
That  is,  the  equations  involve  only  the  four  transient-state  probabilities.  The  evo¬ 
lution  of  the  transient  states  can  thus  be  detennined  independently  of  the  outcome 
states. 

The  outcome  states’  evolution  is,  by  contrast,  completely  determined  by  the  tran¬ 
sient  states  and  the  initial  conditions.  For  example,  given  P4j3(t),  the  entire  history 
of  P4jo(t)  follows  from 

JP4  0(x)dx  =  P4/o(t)-P4,o(0)  =  P4,o(t)  =  4KjP4  3(x)dx  [Eq.  2-8] 

0  0 

where  we  have  used  the  initial  condition  P4,o(0)  =  0. 

Letting  t  — >  go  in  Equation  2-8,  we  see  that  the  long-time  limiting  value  of  P4,o(t) 
is  given  by 


OO 

lim  P4 ,0(t)  =  4k|P4 /3(x)dx  .  [Eq.  2-9] 


We  now  show  how  to  exploit  characteristics  of  the  engagement  we  are  consider¬ 
ing,  shown  by  Figure  2-3  and  reflected  in  properties  of  Equation  2-7,  to  allow  cal¬ 
culation  of  long-time  limiting  probabilities  of  outcome  states  with  a  simple,  finite 
iterative  scheme.  This  engagement  is  an  example  of  what  have  been  called  “pure 
death”  processes:  the  states  always  transition  to  states  with  fewer  participants.  Re¬ 
flecting  this,  Figure  2-3  is  an  “acyclic”  graph.  That  is,  the  graph  gives  no  path  by 
which  to  return  to  a  previously  occupied  state. 

A  consequence  of  the  acyclic  character  of  Figure  2-3  is  found  in  the  structure  of 
Equation  2-7.  Defining  the  vector  x  as  (P44,  P43,  P34,  P33),  we  see  that  the  four 
equations  for  the  transient  states  can  be  written  as 
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x  =  Ax 


[Eq.  2-10] 


where  the  matrix  A  is  given  by 


r  —  (4k +  4)  0  0  0  N 

4k  -(4k +  3)  0  0 

4  0  -(3k +  4)  0 

0  3  3k  -(3k  +  3)v 


[Eq.  2-11] 


The  acyclic  character  of  Figure  2-3  is  reflected  in  the  lower-triangular  character  of 
the  matrix  A.  This  lower-triangular  feature  allows  calculation  of  the  integrals 
from  0  to  go  of  each  of  the  transient-state  probabilities  with  a  simple  iterative 
scheme.  Defining  the  vector  x  by 


x  =  Jx(t)dt,  [Eq.  2-12] 

o 


we  see,  on  integrating  Equation  2-10  from  t  =  0  to  t  =  qo  and  remembering  that 
P4,4(0)  =  1  and  all  other  transient-state  probabilities  are  zero  at  t  =  0,  that 

-ej=Ax,  [Eq.  2-13] 


where  ei  is  the  vector  (1,  0,  0,  0).  Thus  x  is  detennined  by  the  solution  of  a  sys¬ 
tem  of  linear  algebraic  equations.  Moreover,  by  virtue  of  the  lower- triangular 
character  of  A,  that  system  is  easily  solved  with  an  obvious  iterative  scheme.  On 
writing  out  Equation  2-13  in  detail,  we  find 


-1  =  -(4k  +  4)X[ 

0  =  4kx1  -(4k  +  3)x2 
0  =  4xj  -  (3k  +  4)x3 
0  =  3x2  +  3x3  -  (3k  +  3)x4 


[Eq.  2-14] 


To  solve  those  equations,  one  may  obtain  xx  from  the  first  equations,  then  x2  and 
x3  follow  from  the  second  and  third  equations,  and  finally,  x4  follows  from  the 
last  equation.  The  lower-triangular  character  of  the  matrix  A  guarantees  that  the 
nth  equation  in  the  system  of  Equation  2-14  involves  xn  and  does  not  involve  Xj 

for  any  j  >  n.  Another  way  to  view  the  benefits  of  the  lower-triangular  character 
of  the  matrix  A  is  to  observe  that  this  makes  Equation  2-13  have  the  fonn  of  a 
system  of  equations  after  the  forward  course  of  Gaussian  elimination  has  been 
performed,  so  that  only  the  rapid  back  course  remains  to  be  done. 
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Why  were  we  justified  in  integrating  the  x  from  0  to  oo?  Texts  on  differential 
equations4  show  that  solutions  of  equations  such  as  Equation  2-10  are  always  lin¬ 
ear  combinations  of  functions  exp(Ajt),  where  a:!  is  an  eigenvalue  of  the  matrix  A, 
when  these  eigenvalues  are  distinct.  The  eigenvalues  of  a  lower-triangular  matrix 
are  the  elements  on  the  main  diagonal.  Inspection  of  Equation  2-1 1  shows  that 
these  elements  are  distinct,  and  negative.  Thus  solutions  of  Equation  2-7  decrease 
exponentially  with  time,  and  the  integrals  we  use  exist. 

In  Appendix  A,  we  show  that  the  approach  to  calculating  long-time  limiting  prob¬ 
abilities  illustrated  here  can  be  applied  to  many  engagement  models.  We  use  this 
method,  which  from  now  on  we  will  call  “the  method  of  Appendix  A,”  frequently 
in  the  present  SLAACM  and  in  planned  extensions  to  it. 

For  SLAACM,  the  key  result  from  engagement  calculations  is  the  discrete  bivari¬ 
ate  loss  distribution.  This  gives  the  loss  statistics  that  are  propagated  day-by-day 
in  a  campaign.  Figure  2-4  shows  an  example  of  the  discrete  bivariate  loss  distribu¬ 
tion,  for  a  4  vs.  4  engagement  in  which  both  sides  fight  to  annihilation  (breakpoint 
=  0),  and  k  =  1.46.  In  the  figure,  blue  bars  indicate  Blue  win  states  (states  with 
four  Red  losses),  and  red  bars  indicate  Red  win  states  (states  with  four  Blue 
losses).  The  figure  shows  that,  as  should  be  the  case  when  k  >  1,  the  Blues  are 
likely  to  win  the  engagement  and  the  total  probability  of  the  blue  bars  is  distinctly 
larger  than  the  total  probability  of  the  red  bars.  There  is  significant  dispersion  in 
both  Red  and  Blue  losses. 

Figure  2-4.  Bivariate  Loss  Distribution 


The  marginal  distributions  of  Blue  and  of  Red  losses  and  the  total  win  probabili¬ 
ties  are  shown  in  Figure  2-5.  The  marginal  distributions  show  that  the  marginal 


4  W.E.  Boyce  and  R.C.  DiPrima,  Elementary  Differential  Equations ,  Third  Edition  (New 
York:  Wiley,  1976),  Sections  7-6  through  7-9. 
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probabilities  are  not  necessarily  monotone.  They  also  show  that,  in  this  case,  the 
probability  of  4  losses  is  the  largest  Blue  probability. 


Figure  2-5.  Marginal  Loss  Distributions  with  KRR=1.46 
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Such  a  “toe-up”  at  the  greatest  possible  loss  may  seem  counter-intuitive.  In  fact, 
that  feature  is  certain  to  happen  for  kill-rate  ratios  k  sufficiently  close  to  1.  To  see 
this,  consider  a  4  vs.  4  engagement  to  annihilation,  with  k  =  1 .  The  probability  of 
a  Blue  win  must  equal  the  probability  of  a  Red  win,  and  these  probabilities  must 
both  equal  0.5.  The  probability  of  a  Red  win  is  the  probability  of  4  Blue  losses,  so 
the  probability  of  4  Blue  losses  must  be  0.5. 

The  sum  of  the  marginal  probabilities  of  Blue  losses  must  be  1.  Therefore,  none 
of  the  marginal  probabilities  of  0,  1,  2,  or  3  Blue  losses  can  exceed  0.5.  If,  as  in 
fact  is  the  case,  all  the  marginal  probabilities  of  0,  1,  2,  and  3  Blue  losses  are  posi¬ 
tive,  each  of  them  must  be  less  than  0.5,  so  that,  in  this  case  of  equal  strength 
forces  (k  =  1),  the  marginal  probability  of  4  Blue  losses  must  be  the  largest  mar¬ 
ginal  loss  probability. 

The  marginal  loss  probabilities  are  continuous  functions  of  k  (they  are  in  fact  ra¬ 
tional  functions  of  k,  as  may  be  seen  by  applying  the  method  of  Appendix  A),  so 
that,  as  k  increases  from  1,  at  least  for  a  while,  the  most  probable  Blue  loss  must 
be  4.  Figure  2-6  shows  that  the  Blue  toe-up  is  negligible  for  k  =  5. 
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Figure  2-6.  Marginal  Loss  Distribution 


Blue  and  Red  Loss  Probabilities 
4v4  Blue  :  Red,  kb/kr  =  5,  Fight  to  Death 


Scalability  and  Coping  with  Large- 
Dimensioned  Problems 


The  examples  illustrate  that  initial  value  problems  (Equations  2-5  and  2-6)  are 
tractable,  and  that  their  solutions  give  results  that  military  planners  may  find  use¬ 
ful.  These  examples  also  show  why  computational  probabilistic  models  have  not 
been  widely  used:  Even  moderate  values  of  M  and  N  generate  a  lot  of  equations. 
For  example,  an  engagement  of  4  Blues  with  8  Reds,  not  at  all  uncommon,  leads 
to  a  system  of  44  differential  equations.  The  more  sophisticated  models  that  we 
consider  in  Appendixes  B  and  D  can  lead  to  systems  of  thousands  of  equations. 
Computational  probabilistic  models  have  significant  scalability  concerns. 

Until  fairly  recently,  it  was  not  practical  to  deal  with  such  large  state  spaces.  Ad¬ 
vances  on  two  fronts  have  changed  that,  however.  First,  memory  and  speed  of 
widely  available  PCs  now  match  or  exceed  those  of  the  supercomputers  of  only  a 
decade  or  so  ago.  Second,  recently  available  applications  packages  make  it  possi¬ 
ble  to  treat  the  large-dimensioned  probabilistic  models  and  the  many-variable  in¬ 
teger  programming  problems,  to  which  computational  probabilistic  models 
sometimes  lead. 

In  the  development  of  SLAACM,  and  for  offline  analyses,  we  have  used  Markov 
analysis  tools  developed  by  NASA  for  analysis  of  high-reliability  electronics  that 
greatly  facilitate  treating  the  large-dimensioned  initial  value  problems  to  which 
computational  probabilistic  models  often  lead.  Appendix  B  discusses  how  we 
have  applied  these  tools  to  air  combat  analysis. 
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Composite  Engagements 


In  the  current  version  of  SLAACM,  engagements  are  always  composite  engage¬ 
ments  of  a  particular  kind,  in  which  a  flight  of  4  defending  fighters  engages  an 
attack  package  made  up  of  4  advanced  escorts,  4  close  escorts,  and  4  bombers. 

The  defenders  first  engage  the  advanced  escorts.  If  the  defenders  win  that  en¬ 
gagement,  the  surviving  defenders  engage  the  close  escorts.  If  the  defenders  again 
prevail,  defenders  who  survive  the  engagement  with  close  escorts  engage  the 
bombers.  Results  of  this  final  engagement  detennine  the  distribution  of  the  num¬ 
ber  of  bombers  “leaking”  through  the  defense.  Figure  2-7  illustrates  the  opposing 
forces  for  this  composite  engagement. 

Figure  2-7.  Opposing  Forces  in  SLAACM  Composite  Engagement 


) 

) 

Sixteen  final  outcome  or  absorbing  states  are  possible  when  both  sides  fight  to 
annihilation.  The  outcome  states  may  be  described  by  the  ordered  quadruple  (A, 

C,  B,  D),  where  A  is  the  number  of  surviving  Red  advanced  escorts,  C  the  number 
of  Red  close  escorts  that  survive,  B  the  number  of  surviving  Red  bombers,  and  D 
the  number  of  surviving  Blue  defenders.  Table  2-1  lists  the  16  possible  absorbing 
quadruples  for  the  SLAACM  composite  engagement.  Table  2-2  shows  the  64  total 
transition  and  absorbing  states  possible  for  the  scenario. 
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0 

0 

1 

4 

0 

0 

0 

4 

0 

0 

0 

3 

0 

0 

0 

2 

0 
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Table  2-1.  Absorbing  States, 
SLAACM  Composite  Engagement 


A 

Red  LE 

C 

Red  CE 

B 

Red  B 

D 

Blue 

0 

0 

1 

0 

0 

0 

0 

4 

0 

0 

0 

3 

0 

0 

0 

2 

0 

0 

0 

1 

Table  2-2.  Transition  and  Absorbing  States, 
SLAACM  Composite  Engagement 


A 

Red  LE 

B 

Red  CE 

C 

Red  B 

D 

Blue 

4 

4 

4 

4 

4 

4 

4 

3 

4 

4 

4 

2 

4 

4 

4 

1 

4 

4 

4 

0 

3 

4 

4 

4 

3 

4 

4 

3 

3 

4 

4 

2 

3 

4 

4 

1 

3 

4 

4 

0 

2 

4 

4 

4 

2 

4 

4 

3 

2 

4 

4 

2 

2 

4 

4 

1 

2 

4 

4 

0 

1 

4 

4 

4 

1 

4 

4 

3 

1 

4 

4 

2 

1 

4 

4 

1 

1 

4 

4 

0 

0 

4 

4 

4 

0 

4 

4 

3 

0 

4 

4 

2 

0 

4 

4 

1 

0 

4 

4 

0 

0 

3 

4 

4 
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Table  2-2.  Transition  and  Absorbing  States, 
SLAACM  Composite  Engagement 


A 

Red  LE 

B 

Red  CE 

c 

Red  B 

D 

Blue 

0 

3 

4 

3 

0 

3 

4 

2 

0 

3 

4 

1 

0 

3 

4 

0 

0 

2 

4 

4 

0 

2 

4 

3 

0 

2 

4 

2 

0 

2 

4 

1 

0 

2 

4 

0 

0 

1 

4 

4 

0 

1 

4 

3 

0 

1 

4 

2 

0 

1 

4 

1 

0 

1 

4 

0 

0 

0 

4 

4 

0 

0 

4 

3 

0 

0 

4 

2 

0 

0 

4 

1 

0 

0 

4 

0 

0 

0 

3 

4 

0 

0 

3 

3 

0 

0 

3 

2 

0 

0 

3 

1 

0 

0 

3 

0 

0 

0 

2 

4 

0 

0 

2 

3 

0 

0 

2 

2 

0 

0 

2 

1 

0 

0 

2 

0 

0 

0 

1 

4 

0 

0 

1 

3 

0 

0 

1 

2 

0 

0 

1 

1 

0 

0 

1 

0 

0 

0 

0 

4 
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Table  2-2.  Transition  and  Absorbing  States, 
SLAACM  Composite  Engagement 


A 

Red  LE 

B 

Red  CE 

C 

Red  B 

D 

Blue 

0 

0 

0 

3 

0 

0 

0 

2 

0 

0 

0 

1 

The  first  engagement  is  4  defenders  vs.  4  advanced  escorts.  Eight  outcome  states 
are  possible  from  this  engagement:  four  Red  win  states — (4,  4,  4,  0),  (3,  4,  4,  0), 
(2,  4,  4,  0),  and  (1,  4,  4,  0) — in  which  the  advanced  escorts  defeat  the  four  de¬ 
fenders  with  0,  1,  2,  or  3  losses.  Four  Blue  win  outcomes  are  possible:  (0,  4,  4,  4), 
(0,  4,  4,  3),  (0,  4,  4,  2),  and  (0,  4,  4,  1),  when  the  Blues  defeat  the  advanced  es¬ 
corts  with  0,  1,  2,  or  3  losses.  We  use  the  method  of  Appendix  A  to  determine  the 
probabilities  of  these  states.  If  the  advanced  escorts  win  this  engagement,  the 
composite  engagement  ends. 

The  Blue  win  states  cause  the  engagement  to  continue  into  engagements  between 
1,  2,  3,  or  4  defenders  and  the  4  close  escorts.  Four  possible  outcome  states  are 
Red  win  states:  (0,  4,  4,  0),  (0,  3,  4,  0),  (0,  2,  4,  0),  and  (0,  1,  4,  0),  corresponding 
to  close  escort  losses  of  0,  1 ,  2,  or  3  aircraft. 

The  engagement  between  1  defender  and  the  4  close  escorts  has  one  Blue  win 
state:  (0,  0,  4,  1).  The  engagement  between  2  defenders  and  4  close  escorts  has 
two  Blue  win  states:  (0,  0,  4,  1)  and  (0,  0,  4,  2).  The  engagement  between  3  de¬ 
fenders  and  4  close  escorts  has  three  Blue  win  states — (0,  0,  4,  1),  (0,  0,  4,  2),  and 
(0,  0,  4,  3) — and  the  engagement  between  4  defenders  and  4  close  escorts  has  four 
Blue  win  states:  (0,  0,  4,  1),  (0,  0,  4,  2),  (0,  0,  4,  3),  and  (0,  0,  4,  4).  Again,  we  use 
the  method  of  Appendix  A  to  evaluate  the  probabilities  of  all  these  outcomes. 

Thus  the  engagements  with  close  escorts  lead  to  four  Red  win  states  and  to  four 
cases  of  engagements  between  defenders  and  bombers:  1  vs.  4,  2  vs.  4,  3  vs.  4, 
and  4  vs.  4.  These  engagements  lead  to  the  Red  win  states  (0,  0,  1,  0),  (0,  0,  2,  0), 
(0,  0,  3,  0),  and  (0,  0,  4,  0)  and  to  the  Blue  win  states  (0,  0,  0,  1),  (0,  0,  0,  2), 

(0,  0,  0,  3),  and  (0,  0,  0,  4).  Using  the  method  of  Appendix  A  to  evaluate  the  prob¬ 
abilities  of  the  outcome  states  in  the  defenders  versus  bomber  engagements  com¬ 
pletes  determination  of  the  probabilities  of  all  16  outcome  states  in  Table  2-1. 

SFAACM  contains  individual  kill-rate  ratios  for  each  Blue -Red  aircraft  pair  in  the 
Blue  and  Red  inventories.  SFAACM  uses  these  to  calculate  probabilities  for  each 
of  the  16  possible  outcomes  for  all  the  feasible  4  vs.  4+4+4  combinations.  In  a 
campaign,  Red  and  Blue  use  the  resulting  probabilities  to  construct  optimal  attack 
packages  and  optimal  defenses. 
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Extensions  to  Engagement  Models 


We  have  now  reviewed  the  basic  engagement  model,  and  the  composite  engage¬ 
ment  model,  presently  used  in  SLAACM.  The  same  methods  used  here  can  also 
be  used  to  extend  SLAACM  by  using  more  sophisticated  engagement  models, 
such  as  the  ones  described  in  Appendix  D.  Given  more  data  than  the  presently 
available  loss  ratios,  for  example,  we  could  introduce  the  more  realistic  two-phase 
kill  engagement  model  of  Appendix  D.  When  breaking  radar  lock  is  not  consid¬ 
ered,  the  long-time  limiting  values  of  outcome  probabilities  for  this  model  are  de¬ 
termined  by  three  non-dimensional  ratios: 

♦  Ratio  of  the  mean  time  for  lock-on  by  the  Blue  aircraft  to  the  mean  time 
for  lock-on  by  the  Red  aircraft 

♦  Ratio  of  the  mean  time  for  the  Blue  aircraft  to  make  a  kill  after  lock-on  to 
the  mean  time  for  lock-on  by  the  Red  aircraft 

♦  Ratio  of  the  mean  time  for  a  Red  aircraft  to  make  a  kill  after  lock-on  to 
the  mean  time  for  lock-on  by  the  Reds. 

When  breaking  radar  lock  is  considered,  two  more  ratios  are  added. 

Given  at  least  three  outcome  data,  we  could  infer  the  non-dimensional  ratios 
needed  for  the  two-phase  kill  model  without  lock-breaking.  With  at  least  five  out¬ 
come  data,  we  could  infer  the  ratios  needed  for  the  two-phase  kill  model  with 
lock-breaking. 

SLAACM ’s  composite  engagements  also  can  be  extended,  for  example,  by  in¬ 
cluding  the  effects  of  SAM  and  AAA  defenses  (defensive  counter  air  effects  are 
already  included).  Engagement  models  of  the  kind  described  above  for  air-to-air 
engagements  can  also  be  made  for  suppression  of  enemy  air  defense  (SEAD)  ac¬ 
tions;  those  engagements  could  also  be  included  in  SLAACM’s  composite  en¬ 
gagements. 
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Chapter  3 

How  SLAACM  Works: 

Bomber  Effectiveness  Parameter 


The  basic  campaign  modeled  in  SLAACM  is  a  Red  bomber  attack  against  Blue 
defenders.  Red’s  bombers  have  different  levels  of  effectiveness,  based  on  the  at¬ 
tack  package.  The  offensive  effectiveness  of  a  bomber  depends  on  two  factors: 
the  payload  of  bombs  carried,  and  the  accuracy  with  which  those  weapons  can  be 
delivered  to  targets.  We  developed  a  measure  of  bomber  effectiveness  based  on 
these  factors. 

We  considered  first  the  relation  between  blast  overpressure  p  and  the  properties  of 
a  structure  designed  to  resist  that  overpressure.  We  used  the  simple  model  of  rein¬ 
forcing  structure  shown  in  Figure  3-1. 

Figure  3-1.  Reinforcing  Structure  Model 


— H  a  b — 


In  this  simple  model,  two  columns  reinforce  the  structure.  Their  thicknesses  are  A, 
their  heights  are  h,  and  the  spacing  between  their  centers  is  d.  Their  dimension 
into  the  page  is  w. 

The  reinforcing  columns  act  as  cantilevers  against  lateral  loads.  A  pressure  p  act- 

pwh 

mu  on  the  left  face  of  this  simple  structure  will  induce  axial  loads -  in  each 

2d 

column,  tension  in  the  left  and  compression  in  the  right.  The  pressure  will  also 
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pwh 

induce  shear  loads - in  each  column.  Since  both  loads  are  proportional  to  p, 

2 

the  maximum  stress  induced  in  the  columns  will  also  be  proportional  to  p.  If  the 
columns’  dimensions  are  adjusted  so  that  a  material  of  given  maximum  allowable 
stress  can  support  the  loads  induced  by  the  overpressure,  the  area  of  the  columns’ 
bases  wA  will  be  proportional  to  the  maximum  stress  and,  therefore,  also  to  p. 

Thus,  the  volume  of  the  reinforcing  structure,  2wAh,  will  be  proportional  to  p. 
Assuming  that  the  cost  of  the  reinforcing  structure  is  proportional  to  its  volume 
leads  to  the  conclusion  that  the  cost  of  reinforcement  will  be  proportional  to  over¬ 
pressure  p. 

Now  let  us  consider  how  blast  overpressure  varies  with  explosive  yield  E.  Dy¬ 
namic  overpressure  p  is  proportional  to  the  square  of  the  air  velocity  v  immedi¬ 
ately  behind  the  blast  wave  produced  by  the  explosion.  Dimensional  analysis 
shows  how  v  varies  with  E. 1  In  addition  to  E,  the  parameters  affecting  v  are  the 
undisturbed  air  density  p  and  the  time  t  required  for  the  blast  wave  to  reach  the 
locations  of  interest. 

According  to  the  principles  of  dimensional  analysis,  v  can  be  expressed  as  the 
product  of  a  term  E“p^t7  and  a  dimensionless  function  of  all  possible  dimensionless 
combinations  of  E,  p,  and  t.  If  Eup*V  is  equal  to  a  velocity,  that  term  must  have  the 
dimension  L/T,  where  L  represents  dimension  length  and  T  dimension  time.  Now, 
energy  E  has  dimension  ML  / T  ,  where  M  denotes  dimension  mass;  p  has  dimen¬ 
sion  M/L3,  and,  of  course,  t  has  dimension  T.  Thus 


V 


ML 


2  A 


fMY 

T  3 

vL  y 


[Eq.  3-1] 


Equation  3-1  implies  that 

a  +  P  =  0 

2a  -  3P  =  1  .  [Eq.  3-2] 

-2a  +  y  =  -1 


The  unique  solution  to  Equation  3-2  is  a  =  1/5,  P  =  -1/5,  y  =  -3/5.  Similar 
arguments  to  those  leading  to  Equations  3-1  and  3-2  show  that  there  are  no  non- 
dimensional  combinations  of  E,  p,  and  t,  other  than  a  constant.  Therefore, 


V  oc 


( -c  A 


P) 


[Eq.  3-3] 


1  L.  Sedov,  Similarity’  and  Dimensional  Methods  in  Mechanics  (New  York:  Academic  Press, 
1959). 
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and  the  dynamic  overpressure  p  will  satisfy 

2  3  _6 

p  oc  pv2  oc  E5p5t  5.  [Eq.  3-4] 


Similar  dimensional  considerations  lead  to  the  conclusion  that  the  radius  r  of  the 
blast  wave  satisfies  the  proportionality 


r  oc 


E 

\PJ 


i 

^  2 

t5. 


[Eq.  3-5] 


From  proportionalities  (Equations  3-4  and  3-5),  it  follows  that  the  radius  tl  at 
which  the  overpressure  falls  to  Pl  satisfies 


rL 


"en 

v  P  y 


1 

3 

Pl  3. 


[Eq.  3-6] 


Thus,  the  lethal  radius  of  a  blast  wave  increases  as  the  cube  root  of  the  blast’s  en¬ 
ergy  and  decreases  as  the  cube  root  of  the  lethal  overpressure. 


Since  in  our  model  the  value  of  a  structure  is  proportional  to  the  overpressure  re¬ 
quired  to  destroy  it,  proportionality  (Equation  3-6)  implies  that  the  lethal  radius  of 
a  blast  with  energy  E  varies  directly  with  the  cube  root  of  E,  and  inversely  with 
the  cube  root  of  the  value  of  the  structures  it  affects.2 


We  now  use  the  information  about  the  variation  of  lethal  radius  with  explosion 
energy  and  structure  value  to  see  how  bombs’  single  shot  kill  probabilities 
(sspk’s)  vary  with  explosion  energy  and  value  of  structures  affected. 

The  sspk  of  a  bomb  with  lethal  radius  rL,  delivered  with  circular  error  probable 
(CEP),  is 


sspk  =  1  - 


rn 


rL 

CEP2 


[Eq.  3-7] 


Introducing  the  variation  of  tl  with  explosive  energy  E  and  value  V  obtained 
above,  we  have 


sspk  =  1 


kE2/3y-2/3 


CEP2 


[Eq.  3-8] 


2  These  conclusions,  obtained  with  simple  dimensional  analyses,  agree  with  those  from  the 
exact  Taylor-Sedov  solution  [Sedov,  1.  c.  ante]  for  values  immediately  behind  the  shock. 
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where  k  is  a  constant. 


Now, 


kE2/3v-2/3  =  kE2/3y-e2/3 


VEref  2 


V, 


ref 


V 


CEP^ 

CEP 


[Eq.  3-9] 


Defining  the  parameter  ki  by 


u  =1 <F2/3V-2/3 

K1  -KEref  Vref  ’ 


[Eq.  3-10] 


we  have 


sspk  =  1 


f E 

'VrefV 

V  Eref  J  l 

CEPrefC 
CEP  J 


[Eq.  3-11] 


Let  Eref  be  the  energy  of  a  500  lb  bomb.  Let  Vref  be  the  value  of  a  target  hardened 
to  withstand  a  500  lb  bomb.  We  will  consider  500  lb,  1,000  lb,  and  2,000  lb 
bombs  and  assume  that  bombs  are  always  matched  to  targets,  so  that  a  bomber  on 
a  mission  against  500  lb  targets  carries  500  lb  bombs,  and  one  on  a  mission 
against  1,000  lb  targets  carries  1,000  lb  bombs,  and  that  one  going  against  2,000 
lb  targets  carries  2,000  lb  bombs. 

With  those  assumptions,  the  product  of  the  term  involving  E  and  the  tenn  involv¬ 
ing  V  in  Equation  3-11  is  always  1,  and 


sspk  =  1 


CEPrefC 
CEP  J 


12  J 


[Eq.  3-12] 


We  arbitrarily  set  the  constant  ki  to  1 ;  we  only  want  a  systematic  way  to  compute 
the  value  to  Red  of  a  bomber  with  given  payload  P  that  is  capable  of  delivering 
bombs  with  relative  accuracy  CEP/CEPref. 

With  this  assumption,  the  value  of  the  expected  number  of  targets  killed  by  such  a 
bomber  carrying  500  lb  bombs  is 


Vi 


P 

500 


sspkgVjoo  i 


[Eq.  3-13] 
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where 


sspk0  =  1 


rn 


CEPreff 
CEP  J 


U j 


[Eq.  3-14] 


and  where  V500  is  the  value  associated  with  a  target  hardened  to  withstand  500  lb 
bombs. 

When  the  bomber  is  loaded  with  1 ,000  lb  bombs,  the  value  of  the  expected  num¬ 
ber  of  targets  killed  is 


sspk0(2)^5V5M,  [Eq.  3-15] 

since  the  value  of  a  target  goes  as  the  2/5  power  of  the  explosive  yield  required  to 
destroy  it. 

Finally,  if  the  bomber  is  loaded  with  2,000  lb  bombs,  the  value  of  the  expected 
number  of  targets  killed  is 

sspk0(4)2/5V500.  [Eq.  3-16] 

Using  this  information,  SLAACM  assigns  the  average  of  Vi,  V2,  and  V3,  divided 
by  V500,  as  the  value  of  the  bomber. 

SLAACM  currently  assigns  CEP/CEPref  ratio  to  1  if  the  bomber  cannot  deliver 
precision-guided  (smart)  bombs.  With  payloads  of  smart  bombs  the  CEP/CEPref 
will  be  less  than  1  and  the  (CEPref/CEP)  tenn  in  the  sspko  equation  will  be  greater 
than  1 .  If  specific  values  of  CEP  are  available,  they  can  be  used  in  the  calculation 
of  bomber  value. 

Although  bombers  are  considered  to  have  poor  defensive  capability,  they  do  have 
kill-rate  ratios  assigned  and  their  defensive  capability  is  considered  along  with 
their  offensive  capability  in  the  payoff  calculations  (discussed  in  the  next  chap¬ 
ter). 


V2  = 


2000 
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Chapter  4 

How  SLAACM  Works: 

Campaign  Logic  and  Optimization 


This  chapter  discusses  the  Red  and  Blue  payoff  logic  and  optimization  algorithms 
used  in  SLAACM. 

Specifying  Red  and  Blue  Payoff  Functions 

SLAACM  determines  optimal  Red  attacks  and  optimal  Blue  responses  using  pay¬ 
off  functions  for  each  side.  This  section  describes  those  calculations. 

Red’s  Payoff  Function 

Red  values  an  attack  package  by  considering  the  benefit  of  getting  bombers 
through  and  of  downing  Blue  aircraft.  Red  may  also  consider  penalties  for  losing 
its  own  aircraft,  assigning  different  values  to  losing  a  fighter  and  to  losing  a 
bomber.  Specifically,  a  package’s  value  to  Red  is  detennined  in  this  way:  A  Red 
package  is  defined  by  the  triple  (a,  c,  b),  where  a  is  the  advanced  escort  type,  c  the 
close  escort  type,  and  b  the  bomber  type  (the  bombers  may  be  fighter/bomber  air¬ 
craft).  The  outcome  statistics  for  an  engagement  between  the  package  and  the 
“planning  Blue”  aircraft  (the  user  identifies  a  Blue  type  for  Red’s  planning)  are 
determined,  and,  from  these,  the  expected  number  Bs  of  surviving  bombers,  the 
expected  number  RFk  of  lost  fighters,  and  the  expected  number  Dk  of  downed 
Blues  are  computed.  Then  VR,  Red’s  value  for  the  package,  is  given  by 

vr  =  ciBs  +  c2Dk  -c3RFk  -c4(4-Bs)  .  [Eq.  4-1] 

The  positive  parameter  ci  is  currently  set  at  the  value  of  a  Red  bomber,  computed 
by  the  method  described  in  Chapter  3.  The  value  of  C2  is  set  at  20,  reflecting  the 
fact  that  Red  believes  Blue  to  be  loss  averse  and  likely  to  abandon  the  campaign  if 
their  side  incurs  many  losses.  Currently,  C3  and  C4  are  both  set  at  zero.  When  the 
Blue  aircraft  are  significantly  stronger  than  the  Red  aircraft  (values  of  k  signifi¬ 
cantly  greater  than  1),  even  values  of  C3  and  C4  that  are  small  compared  to  1  may 
cause  the  Reds  not  to  dispatch  any  packages  at  all. 

Blue’s  Payoff  Function 

Blue  values  engaging  a  given  Red  package  by  considering  the  benefit  of  downing 
Red  fighters  and  bombers,  and  the  penalty  of  losing  its  own  aircraft.  Specifically, 
the  value  to  Blue  of  engaging  a  Red  package  with  fighters  of  a  given  type  is  cal- 
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culated  by  computing  the  outcome  statistics  for  that  composite  engagement  (de¬ 
scribed  in  Chapter  3)  and,  from  them,  the  expected  number  Bk  of  bombers  killed, 
the  expected  number  RFk  of  Red  fighters  killed,  and  the  expected  number  Dk  of 
defenders  killed.  Then  the  value  Vb  of  the  engagement  to  Blue  is  evaluated  as 

VB=diBk-d2Dk+d3RFk.  [Eq.  4-2] 

Currently,  di  is  set  to  Red’s  value  for  the  bomber  type  in  the  package,  d2  is  set  to 
20,  and  d3  is  set  to  1. 

Values  of  ci,  c2,  C3,  C4,  di,  d2,  and  d3  are  not  currently  user  inputs  to  SLAACM; 
they  may,  however,  readily  be  changed  by  simple  modifications  to  the  code. 

Determining  Attacks  and  Defenses 


A  key  part  of  SLAACM ’s  operation  is  determining  the  engagements  that  occur 
each  day.  This  is  done  by  finding  optimal  sets  of  attack  packages  dispatched  by 
Red  and  by  determining  Blue’s  response  in  two  ways:  the  response  of  smart  Blues 
who  know  the  makeup  of  Red’s  packages  before  intercepting  them  and  make  op¬ 
timal  defenses,  and  the  response  of  other  not-smart  Blues  who  encounter  Red 
packages  randomly.  This  section  explains  these  calculations. 

Determining  Red’s  Attack  Packages 

As  the  aggressor,  Red  always  plans  its  attacks  in  advance.  This  leads  to  an  integer 
programming  (IP)  problem,  which  SLAACM  treats  in  one  of  two  ways:  a  classi¬ 
cal  “greedy”  heuristic  or  an  exact  solution  using  the  commercial  general  optimiza¬ 
tion  package  LINGO™.  This  section  describes  those  calculations. 

Red’s  Integer  Programming  Problem 

Each  day,  Red  has  a  set  of  fighter,  bomber,  and  fighter/bomber  aircraft  to  deploy. 

In  SLAACM,  Red  chooses  the  set  of  attack  packages  to  maximize  the  total  of  the 
expected  Red  payoff  of  the  packages  dispatched  (evaluating  payoff  as  described 
above  in  “Red’s  Payoff  Function”),  constrained  by  the  available  forces. 

Here  we  fonnalize  the  IP  problems  perfonned  in  the  LINGO™  model  within 
SLAACM.  The  main  idea  behind  using  an  IP  tool  is  to  solve  the  following  optimiza¬ 
tion  problem:  Red  wishes  to  come  with  optimal  attack  packages  based  on  a  determined 
payoff  function  (this  is  the  first  optimization).  Then,  the  Blue  aircraft  that  are  sophisti¬ 
cated  enough  to  have  a  priori  knowledge  of  what  Red  attack  packages  are  coming  may 
optimize  their  defending  packages  accordingly  (the  second  optimization). 

A  Red  attack  package  consists  of  4  advanced  escorts,  4  close-in  escorts,  and  4 
bombers.  Since  some  fighters  may  be  bombers,  we  could  have  a  set  of  possible 
attack  packages  T  ,  where,  given  R  types  of  Red  fighters  and  bombers,  we  could 
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have  |T|  =  R3 .  In  reality,  however,  most  mathematically  feasible  combinations  of 

packages  are  unrealistic  for  warfare  (one  would  not,  for  example,  dispatch  a  pack¬ 
age  of  heavy  bombers  escorting  fighters).  Since  the  number  of  reasonable  attack 
packages  is  much  smaller  than  the  complete  enumeration,  we  assume  those  unde¬ 
sirable  packages  are  removed  in  advance  of  the  optimization;  then,  we  limit  our¬ 
selves  to  those  remaining,  reasonable  possibilities  in  the  integer  program.  The 
degree  to  which  one  wishes  to  limit  the  number  of  potential  packages  depends 
primarily  on  the  size  of  the  integer  program,  which  is  a  function  of  the  size  of  R  . 

We  define  a  set  of  attack  packages  for  Red  as  IcT.  We  define  variable  ry  as  the 
number  of  Red  aircraft  of  type  j  used  in  attack  package  i .  To  calculate  its  payoff,  Red 

plans  its  attack  assuming  each  package  will  be  confronted  by  one  specific  Blue  de¬ 
fender  type,  usually  the  most  numerous  of  the  Blues.  Based  on  the  payoff  functions 
described  previously,  each  attack  package  for  Red  has  a  certain  payoff,  denoted  as  p, . 

We  denote  the  number  of  Red  aircraft  of  type  j  as  n  j .  At  this  point,  Red  can  solve  for 

its  optimal  attack  strategy,  denoted  by  the  following  integer  optimization  problem: 

i 

max^p.x, 

i=l 

I 

st*  Sruxi  - nj  vj’eJ- 

i=l 

Xj  e  {0,1,2,...} 

We  solve  for  the  values  of  x(  ’s,  which  represent  the  optimal  number  of  packages 

of  type  i  sent.  This  problem  is  a  variant  of  the  classic  knapsack  problem.  The  so¬ 
lution  to  this  integer  programming  problem  provides  the  optimal  set  of  attack 
packages  dispatched  by  Red. 

Heuristic  Solution 

The  exact  integer  programming  problem  can  be  solved  using  an  IP  software  tool 
such  as  LINGO.  In  SLAACM,  the  user  can  specify  that  the  exact  solution  tech¬ 
nique  to  be  used,  provided  the  user’s  PC  has  a  LINGO  license.  Although  utilizing 
an  optimization  package  is  the  ideal  circumstance — it  can  guarantee  optimal  solu¬ 
tions  if  run  to  completion — many  analysts  will  not  have  LINGO  available  on  their 
PCs.  For  this  reason,  SLAACM  includes  a  greedy  heuristic  using  Red  and  Blue 
pay-off  functions  that  attempts  to  approximate  the  optimum  solution.  The  heuris¬ 
tic  should  yield  a  near-optimum  solution  when  the  Red  and  Blue  aircraft  types 
and  numbers  are  limited.  However,  as  the  choices  grow,  and  particularly  when 
Red  has  options  of  using  aircraft  as  fighters  or  fighter/bombers,  the  packages  se¬ 
lected  by  the  heuristic  will  deviate  significantly  from  an  IP  solution.  This  differ¬ 
ence  will  be  compounded  when  Blues  are  smart  and  also  optimize.  (When  Blue 
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aircraft  are  dominant,  the  different  package  choices  may  produce  only  small  dif¬ 
ferences  in  the  outcome  statistics.)1 

The  greedy  heuristic  does  what  might  be  considered  the  “obvious”  answer.  It  is 
considered  greedy  because  it  takes,  without  more  holistic  considerations,  pack¬ 
ages  with  the  largest  payoffs  first.  This  approach  can  be  suboptimal  because  it 
does  not  consider  (or  mathematically  eliminate  from  consideration)  all  combina¬ 
tions  of  package  types. 

As  we  noted  earlier,  the  payoff  to  Red  for  attack  package  i  is  denoted  as  p;.  The 
heuristic  first  orders  the  packages  from  greatest  to  least.  Without  loss  of  general¬ 
ity,  we  denote  the  largest  payoff  as  package  pi,  the  second  largest  as,  P2,  etc. 

Thus,  we  have  Pj  <  p2  <  ...  <  p, .  Again,  we  define  rij  as  the  number  of  Red  air¬ 
craft  of  type  j .  The  heuristic  can  then  be  described  by  the  following  procedure. 

Step  0:  Set  nf  =  n^  for  all  j . 

Initialize  x;  =  0  for  alii. 

Set  i  =  1 . 

Step  1 :  If  n  j  - 1  >  0  for  all  j  such  that  r,j  >  0,  then  { 

n  j  =  n  •  - 1  for  all  j  such  that  ry  >  0 
x,  =Xi+l 

Go  to  Step  1 
}  Else  { 

Go  to  Step  2  } 

Step  2:  i  =  i  +  l. 

If  Pi  <  0 ,  quit. 

If  i  =  I  + 1 ,  quit. 

Else,  Go  to  Step  1 . 

In  essence,  this  heuristic  keeps  selecting  the  package  with  the  greatest  payoff  until  it 
exhausts  one  of  the  package’s  aircraft  types  (i.e.,  “hits  a  constraint”).  It  then  looks  for 
the  package  with  the  next  greatest  payoff  (assuming  all  its  aircraft  types  are  still  avail¬ 
able).  It  keeps  selecting  that  package  until  one  of  its  aircraft  types  is  exhausted,  and 
so  forth.  It  continues  this  procedure  until  there  are  no  more  packages  with  positive 


1  For  further  information  on  integer  programming  and  greedy  heuristics,  the  interested  reader 
should  consider  one  of  the  many  textbooks  on  integer  programming  and  optimization,  such  as 
L.  Wolsey  and  G.  Nemhauser,  Integer  and  Combinatorial  Optimization  (New  York:  Wiley,  1988). 
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payoffs  for  which  the  advanced  escort,  close-in  escort,  or  bomber  aircraft  types  have 
not  been  exhausted. 

It  is  important  to  reemphasize  that  this  approach  will  not  necessarily  yield  an  op¬ 
timal  solution.  Indeed,  for  sample  campaigns  modeled  in  SLAACM,  the  above 
heuristic  typically  does  not  yield  the  same  solution  for  Red  as  the  IP  solution 
solved  to  optimality  in  LINGO.  For  the  purposes  of  SLAACM,  it  appears  that  it  is 
rare  for  the  heuristic  solutions  to  differ  by  more  than  a  few  percent  from  optimal¬ 
ity.  Nevertheless,  theoretical  bounds  on  such  a  heuristic  for  a  general  IP  problem 
can  be  rather  far  from  optimality. 

Exact  Solution 

Red’s  IP  problem  may  be  solved  exactly  using  an  integer  programming  solver. 
LINGO  is  a  commercially  available  general  linear,  nonlinear,  and  integer  pro¬ 
gramming  solver  tool.2  The  version  embedded  in  SLAACM  is  Extended 
LINGO™  9.0,  which  does  not  limit  the  number  of  variables  (integer  or  otherwise) 
or  constraints.  Of  course,  for  a  large  enough  problem,  even  sophisticated  IP 
solvers  will  have  trouble  producing  the  optimal  solution  in  reasonable  amounts  of 
time.  The  example  problems  run  in  SLAACM,  but  rarely  take  more  than  a  couple 
of  seconds  on  a  modem  desktop  PC. 

The  IP  model  below  is  written  in  LINGO’s  modeling  language  and  implemented 
in  SLAACM  for  Red’s  optimization.  The  data  and  variables  are  defined  elsewhere 
in  SLAACM.  For  the  sake  of  clarity,  we  omit  portions  of  the  data,  variable  defini¬ 
tions,  and  input/output  declarations. 

MODEL: 

SETS: 

RPACKAGE/1  .rposize/:  rpo,  X; 

RED  AC/1  .rac/:  rob; 

RMATRIX(RPACKAGE,REDAC):  rcnst; 

ENDSETS 

!  INTEGER  PROGRAMMING  PROBLEM; 

MAX  =  VALUE; 

VALUE  =  @SUM(RPACKAGE:  rpo*X); 

@FOR(  REDAC(J): 

@SUM(RPACKAGE(I):  rcnst(l,J)*X(l))  <=  rob(J) ); 

@SUM(RPACKAGE(I):  X(l))  <=  (redLim/12); 

@FOR(  R PACKAGE:  @GIN(X)); 

END 


2  Lindo  Systems,  Inc.,  http://www.lindo.com/. 
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Determining  Blue’s  Response 

Certain  Blue  aircraft  are  identified  by  the  user  as  smart.  They  are  assumed  to 
know  the  composition  of  Red’s  attack  packages  before  intercepting  them  and  to 
coordinate  their  defense  to  maximize  its  total  value  (value  computations  are  de¬ 
scribed  above  in  the  section  “Blue’s  Payoff  Function”).  The  remaining  Blue  air¬ 
craft  encounter  Red  packages  randomly.  The  following  subsections  describe  how 
SLAACM  determines  the  set  of  engagements  made  by  smart  Blues  and  by  the 
remaining  Blue  aircraft. 

Blue’s  Integer  Programming  Problem 

Given  the  set  of  Red  attack  packages,  the  smart  Blues  detennine  which  types  to 
engage  by  maximizing  the  total  Blue  payoff,  subject  to  constraints  imposed  by 
their  order  of  battle.  As  stated  previously,  Blue  responds  with  optimal  defenses  of 
Red’s  attack  if  a  given  aircraft  is  smart  enough  to  have  a  priori  knowledge  of 
Red’s  attack  package  composition.  We  assume  Red  has  obtained  its  solution  be¬ 
fore  Blue  calculates  its  optimal  response.  We  denote  the  optimal  solution  for 
Red — that  is,  (al,a2,...,aI)  =  {x[  ,x2,...,x])  — as  constraints  for  Blue’s  optimiza¬ 
tion.  In  other  words,  a;  is  the  number  of  Red  attack  packages  of  type  i  that  Red 
elected  to  send. 

Blue  has  its  own  payoff  function  for  each  type  of  potential  Red  attack  package. 
We  denote  pik  as  the  payoff  to  Blue  for  intercepting  a  Red  package  i  with  the 

Blue  4-ship  package  of  aircraft  type  k  .  Again,  the  components  of  the  payoff 
function  have  been  described  previously;  they  include  a  loss-aversion  factor  for 
Blue,  along  with  positive  payout  for  expected  Red  kills  and  bomb  damage  pre¬ 
vented.  We  denote  the  set  of  Blue  aircraft  types  with  a  priori  knowledge  of  Red’s 
attack  packages  as  K'cK.  We  denote  the  number  of  Blue  aircraft  of  type 
keK'  as  mk  .  We  then  formulate  and  solve  the  following  IP  problem  for  Blue. 

In  short,  we  want  to  solve  for  the  optimal  number  of  Blue  4-ship  package  of  type 
k  that  intercepts  Red  attack  package  i ,  denoted  by  yik  . 

I  K1 

max  ElXy* 

i=i  j=i 

s.t.  X4ylk<mk  VkeK’ 


K’ 

IX 


k=l 


<  a; 


Vi  e  I 
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yik  e  {0,1,2,...}. 


[Eq.  4-3] 


The  IP  model  written  in  LINGO’s  modeling  language  is  embedded  into 
SLAACM.  LINGO  solves  the  IP  problem  and  returns  the  solution  to  variables 
defined  in  SLAACM.  At  this  point,  the  remaining  not-smart  Blue  aircraft,  denoted 
by  K'cK,  where  K' u K"  =  K  and  K'  n  K"  =  0 ,  engage  the  Red  packages  that 
were  not  selected  by  the  smart  Blues.  That  is,  the  Blue  aircraft  without  a  priori 
knowledge  engage  randomly  the  Red  packages  that  remain  after  the  above  IP 
problem  is  solved,  i.e.,  those  Blue  aircraft  engage  a  set  of  attack  packages 

f  K  V 


,  where  a'  =  max  a;  ~^yik,0 


.  Those  encounters  are  described  in 


V  k=l  J 


a  later  section. 


Heuristic  Solution 


As  stated  previously,  SLAACM  allows  the  user  to  specify  whether  to  use  LINGO 
to  solve  the  IP  formulation  or  to  use  a  greedy  heuristic.  If  SLAACM  uses  the  heu¬ 
ristic  (or  solves  the  IP  problem  in  LINGO),  it  must  do  so  for  both  Red’s  and 
Blue’s  optimization. 

The  greedy  heuristic  for  Blue  works  similarly  to  the  heuristic  for  Red.  Following 
the  notation  from  the  IP  formulation,  we  denote  pik  as  the  payoff  to  Blue  for  in¬ 
tercepting  a  Red  package  i  with  the  Blue  4-ship  package  of  aircraft  type  k  .  We 
then  order  the  pik  from  greatest  to  least  over  all  i  and  k  .  We  denote  that  new  set 

of  payoffs  as  Pj  >  p2  >...  >  pB  ,  and  note  their  associated  i  and  k  . 

It  is  worth  noting  that  Blues  that  are  not  members  of  the  local  air  defense  (LAD) 
are  presumed  to  be  loss  averse,  and  it  is  not  uncommon  for  some  pik  <  0 ,  i.e.,  in¬ 
tercepting  certain  Red’s  attack  packages  may  have  negative  payoffs  for  Blue.  This 
result  occurs  because  loss-averse  Blues  suffer  a  large  penalty  for  their  expected 
losses.  Neither  the  heuristic  nor  the  LINGO  model  uses  any  aircraft  whose  payoff 
values  result  in  pik  <0,  since  selecting  those  combinations  would  of  course  re¬ 
duce  Blue’s  objective  function. 

Again,  we  denote  the  number  of  Red  packages  of  type  i  sent  as  .  As  before,  we 
denote  the  number  of  Blue  aircraft  of  type  keK'  as  mk  .  The  set  of  decision 

variables  (number  of  Blue  4-ship  packages  of  type  k  that  intercept  Red  attack 
package  i )  is  denoted  by  yik  .  As  with  Red,  we  wish  to  solve  a  multidimensional 
knapsack  problem. 
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We  can  then  describe  the  optimization  heuristic  we  employ  by  the  following  algo¬ 
rithm: 


Step  0:  Set  mk  =  mk  for  all  k  . 

Seta)  =  a;  for  alii. 

Initialize  yik  =  0  for  all  i,k. 

Begin  with  b  =  1 . 

Step  1:  Let  i,k  be  the  values  for  p.-  such  that  p.k  =  pb 
If  fhk  -  4  >  0  and  a)  >  0 ,  Then  { 
fhk  =ihk  -4 
a,  =  as  -1 
yik  =  ylk+1 

Go  to  Step  1 
}  Else  { 

Go  to  Step  2} 

Step  2:  b  =  b  +  l. 

If  pb  <  0,  quit. 

If  b  =  B  + 1 ,  quit. 

Else,  Go  to  Step  1 . 

In  essence,  this  heuristic  finds  the  Blue  k  vs.  Red  i  solution  with  the  highest  pay¬ 
off  to  Blue;  notes  the  associated  Red  and  Blue  package  types  (i.e.,  i  and  k  ,  re¬ 
spectively);  checks  to  see  if  there  are  4  Blue  aircraft  of  type  k  available;  checks  to 
see  if  Red  sent  any  packages  of  type  i;  and  chooses  those  Blue  aircraft  to  send  af¬ 
ter  that  Red  package.  It  sends  the  lesser  of  the  number  of  Blue  aircraft  of  that  type 
still  available  or  the  number  of  Red  packages  of  that  type  that  were  sent.  Then,  the 
heuristic  looks  for  the  next  highest  payoff  and  repeats  the  above  process  until  all 
smart  Blue  aircraft  with  positive  payoffs  are  assigned  or  until  there  are  no  more 
Red  packages  to  assign  a  smart  Blue  4-ship  package. 

As  this  process  is  similar  to  the  one  for  the  Red  heuristic,  it  needs  to  be  stated 
again  that  this  method  does  not  guarantee  optimality  for  Blue,  because  it  does  not 
consider  the  optimization  of  the  entire  Blue  fleet  over  all  possibilities,  but  rather 
looks  for  a  greedy  local  solution.  That  said,  both  Blue’s  and  Red’s  heuristics  are 
computationally  easy  to  implement.  Indeed,  without  the  ability  to  solve  moder¬ 
ately  complex  integer  programs,  these  techniques  can  fonn  good  approximations. 
However,  neither  heuristic  is  able  to  provide  tight  bounds  on  its  performance,  at 
least  for  a  general  knapsack  problem.  Therefore,  the  preferred  method,  given  the 
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choice  of  the  two  methods,  is  to  solve  both  optimizations  exactly  using  an  IP 
solver, 

Exact  Solution 

Blue’s  IP  problem  may  be  solved  exactly  using  LINGO.  In  general,  Blue’s  opti¬ 
mization  is  a  little  more  complex  than  Red’s,  because  it  requires  assigning  Blue’s 
defending  packages  to  Red’s  attack  packages.  Red  simply  creates  optimized  at¬ 
tack  packages,  so  Blue  is  effectively  optimizing  over  two  indices  (which  packages 
should  be  sent,  and  whom  they  should  intercept).  In  most  cases,  however,  the 
number  of  variables  in  Blue’s  optimization  is  not  larger  than  the  number  of  deci¬ 
sion  variables  in  Red’s  optimization,  because  Red  sends  only  a  small  fraction  of 
its  potential  package  types.  In  both  cases,  LINGO  obtains  the  optimal  solution  to 
the  IP  problem  in  a  few  seconds  on  a  standard,  modern  desktop  PC  for  most  sam¬ 
ple  problems  in  SLAACM. 

The  IP  model  written  below  in  LINGO’s  modeling  language  solves  Blue’s  opti¬ 
mization  problem.  It  is  similar  in  syntax  to  Red’s  LINGO  model.  Again,  the  data 
and  the  variable  names  are  defined  elsewhere  in  SLAACM.  For  the  sake  of  clar¬ 
ity,  we  omit  most  data  and  variable  definitions,  as  well  as  input/output  calls. 

MODEL: 

SETS: 

PPACKAGE/1  .bsize/:  sboob; 

REDAC/1  .rsize/:  rsent; 

PMATRIX(REDAC,  PPACKAGE):  bpo,  X; 

ENDSETS 

!  INTEGER  PROGRAMMING  PROBLEM; 

!  Objective  Function; 

MAX  =  VALUE; 

VALUE  =  @SUM(PMATRIX:  bpo*X); 

!  subject  to; 

@FOR(  PPACKAGE(J):  @SUM(REDAC(I):  X(I,J)*4)  <=  sboob(J) ); 

@FOR(  REDAC(I):  @SUM(PPACKAGE(J):  X(I,J))  <=  rsent(l) ); 

@FOR(PMATRIX(l,J):  @GIN(X(I,J))); 

END 

Random  Engagements  by  “Not-Smart”  Blue  Aircraft 

The  optimal  engagements  made  by  smart  Blue  aircraft  may  leave  some  Red  pack¬ 
ages  unengaged.  For  determining  the  packages  that  need  to  be  intercepted  by  not 
smart  Blues,  we  assume,  optimistically  for  Blue,  that  a  Red  package  intercepted 
by  a  smart  Blue  flight  does  not  continue  the  engagement.  This  assumption  is  rea- 
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sonable  when  the  smart  Blues  are  much  stronger  than  any  of  Red’s  fighters  and 
bombers. 

Other  SLAACM  versions  have  considered  the  opposite  polar  case,  which  may  be 
implemented  in  SLAACM  with  simple  reprogramming,  where  the  Red  forces  that 
survive  the  smart  Blue  defense  reorganize  into  optimal  attack  packages  and  con¬ 
tinue  the  assault. 

Blue  aircraft  not  identified  as  smart  encounter  Red  packages  randomly.  They  do 
not  conduct  optimized  defensive  encounters.  There  is  still  an  element  of  battle 
management  for  these  aircraft,  however,  and  that  is  the  efficiency  with  which  they 
can  locate  and  engage  the  Red  attack  packages. 

With  perfect  battle  management,  if  the  number  of  defending  flights  is  at  least  as 
large  as  the  number  of  attack  packages,  every  attack  package  would  be  inter¬ 
cepted.  The  intercepting  aircraft  types  would  be  random,  however,  rather  than  op¬ 
timal.  If  the  number  of  attack  packages  is  larger  than  the  number  of  defending 
flights,  every  defending  flight  would  engage  an  attack  package. 

With  less-than-perfect  battle  management,  not  every  attack  package  would  be  in¬ 
tercepted  in  the  former  case,  and  not  every  defending  flight  would  be  engaged  in 
the  latter  case.  SLAACM ’s  battle  management  feature  allows  the  user  to  enter  a 
“goodness”  parameter  that  characterizes  the  effectiveness  of  battle  management 
for  Blue  aircraft  that  are  not  smart.  The  following  paragraphs  explain  how  that 
parameter  affects  SLAACM’ s  calculations. 

Let  K  distinct  types  of  defender  flights  deal  with  an  attack  by  J  distinct  types  of 
attack  packages.  The  defenders  are  assumed  not  to  know  the  makeup  of  individual 
attack  packages  before  interception,  so  that  the  type  of  attack  package  engaged  by 
a  given  defending  flight  is  the  result  of  random  selection  from  the  set  of  attack 
packages. 

Let  m;  be  the  number  of  defender  flights  of  type  i,  and  let  nj  be  the  number  of  at¬ 
tack  packages  of  type  j .  Then  the  total  number  of  attack  packages,  N,  and  the  total 
number  of  defending  flights,  M,  are  given  respectively  by 

N  =  2>j;M  =  Xmi.  [Eq.  4-4] 

1  1 

First  we  consider  the  case  of  perfect  battle  management.  If  M  >  N,  then  every  at¬ 
tack  package  will  be  intercepted.  Not  all  defending  flights  engage;  in  this  simple 
analysis,  we  assume  that  the  fraction  N/M  of  each  defending  flight  type  engages.3 


3  We  recognize  that  this  simple  analysis  treats  a  complex  combinatoric  problem  crudely. 
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Then,  a  measure  of  E- ,  the  central  tendency  of  the  number  E,j  of  mi  vs.  nj  en¬ 
gagements,  is 


N  nj  mjiij 

Ej;  =  - ITl,  -  =  - 

1]  M  1  N  M 


[Eq.  4-5] 


For  this  simple  analysis,  we  take  E-  to  be  the  number  of  m,  vs.  nj  engagements. 

Summing  E^  over  j  shows  that  the  fraction  of  each  defending  flight  type  engaged 

is  N/M,  as  it  should  be;  summing  over  i  shows  that  all  attacking  flights  of  each 
type  are  engaged. 

When  M  <  N,  every  defending  flight  engages,  but  not  all  attack  packages  can  be 
engaged.  The  estimate  for  E-  analogous  to  the  one  given  in  Equation  4-5  is 


Ejj  =  nij  — 

11  N 


n  j  irqrij 


N 


[Eq.  4-6] 


Summing  this  estimate  for  E^  over  i  shows  that  the  fraction  M/N  of  each  attack 

package  type  is  engaged;  summing  over  j  shows  that  every  defending  flight  of 
each  type  is  engaged. 

The  Battle  Management  Efficiency  Factor 

In  some  cases,  M  and  N  are  both  0(10“).  For  such  large  M  and  N,  providing  this 
“perfect”  defender’s  battle  management  could  exceed  the  capabilities  of  available 
systems.  We  account  for  the  limitations  of  the  defender’s  battle  management  with 
a  simple  adaptation  of  the  “perfect”  case.  We  make  the  adaptation  by  introducing 
“ghost”  attack  packages  or  defender  flights.  A  defender  flight  that  engages  a  ghost 
attack  package  does  not  actually  engage;  an  attack  package  engaged  by  a  ghost 
defender  flight  is  not  actually  intercepted. 

When  M  >  N,  we  introduce  mg  defending  flights  of  type  “ghost,”  and  proceed  as 
in  the  above  analysis  (where  the  defender  has  perfect  battle  management).  This 
gives  us  a  new  estimate  for  Ey : 


_  rruri:  nuri: 

E  = _  1  = _ - _ Lb 

1]  M  +  nw  M  1 

o 


where 


bi 


M 

M  +  mg  ’ 


[Eq.  4-7] 


[Eq.  4-8] 
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When  M  <  N,  we  introduce  ng  attack  packages  of  type  “ghost,”  and  find 


E  -^L 
1]  N  +  ns 


minj 

M 


'2  5 


[Eq.  4-9] 


where 


b7  =  — — —  .  [Eq.  4-10] 

N  +  ng 

Noting  the  similarity  of  Equations  4-8  and  4-10,  we  introduced  into  SLAACM  the 
simple  model 


-  mini 

E::  =  — — J-b  [Eq.  4-11] 

]  M 

and  allow  the  user  to  choose  the  “battle  management  efficiency  factor”  b,  where  0 
<  b  <  1. 

Once  seen,  Equation  4-11  is  so  simple,  and  seemingly  obvious,  that  it  is  fair  to  ask 
why  we  did  not  simply  introduce  it  into  SLAACM  without  analysis.  The  reason  is 
that  we  wanted  to  understand  what  that  simple  choice  implied  about  how  the  bat¬ 
tle  proceeded.  The  analysis  given  here  supplies  that  understanding. 
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Chapter  5 

How  SLAACM  Works: 

SLAACM  Campaign  Calculations 


This  chapter  describes  the  calculation  protocols  used  by  SLAACM  for  campaign 
analysis. 

Determining  and  Propagating 
Loss  and  Destruction  Statistics 


The  calculations  described  in  Chapter  4,  under  “Determining  Attacks  and  De¬ 
fenses,”  generate  the  sets  of  engagements  that  take  place  in  each  day’s  fighting. 
There  are  then  specified  numbers  of  engagements  between  a  specific  Red  attack 
package  (with  its  specified  aircraft  types  for  the  advanced  escorts,  close  escorts, 
and  bombers)  and  a  flight  of  4  defenders  of  a  specific  type.  SLAACM  uses  the 
method  described  in  Chapter  2,  under  “Composite  Engagements,”  to  generate 
probabilities  for  all  16  outcome  states  of  each  kind  of  engagement  that  occurs. 

Then,  to  evaluate  Red  losses,  SLAACM  calculates  expected  values  of  the  losses 
for  advanced  escorts,  close  escorts,  and  bombers.  These  expected  values  are  ac¬ 
cumulated  for  each  Red  aircraft  type,  for  engagements  between  smart  Blues  and 
the  attack  packages,  and  for  engagements  between  non-smart  Blues  and  the  attack 
packages  not  engaged  by  smart  defender  aircraft. 

The  total  expected  losses  for  each  Red  type  are  rounded  to  integers  and  subtracted 
from  the  Red  order  of  battle.  SLAACM  then  generates  the  Red  order  of  battle  for 
the  following  day  using  the  surviving  Reds  plus  any  reinforcements  scheduled  to 
arrive  for  that  day.  Reinforcements  are  a  user  input  to  SLAACM. 

Expected  values  of  surviving  bombers  are  multiplied  by  bomber  payloads  and 
rounded  to  the  nearest  integer,  to  give  tons-of-bombs  results.  Tons  delivered  by 
bombers  with  relative  CEPs  less  than  a  user-specified  value  are  identified  as 
“smart”  and  accounted  for  separately  from  tons  delivered  by  other  bombers. 

In  treating  Blue  losses,  SLAACM  considers  measures  of  dispersion  as  well  as 
measures  of  central  tendency.  The  probability  distributions  for  outcome  states  of 
the  engagements  occurring  each  day  are  used  to  calculate  both  means  and  vari¬ 
ances  of  the  losses  of  each  type  of  Blue  aircraft.  Mean  losses,  rounded  to  the 
nearest  integer,  are  subtracted  from  Blue’s  order  of  battle  and — after  increments 
from  any  scheduled  reinforcements — are  posted  as  Blue’s  order  of  battle  for  the 
next  day. 
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Variances  in  the  losses  of  each  type  of  Blue  aircraft  are  accumulated  day  by  day, 
and  the  total  variances  are  reported  to  give  a  measure  of  dispersion  in  Blue  losses. 

This  procedure  neglects  a  potentially  important  source  of  dispersion:  the  effect  of 
propagating  losses  other  than  the  means  of  both  Red  and  Blue  losses.  In  the  cases 
that  we  have  considered,  the  Blue  aircraft  are  always  much  stronger  than  the  Red 
aircraft.  Consequently,  there  is  little  dispersion  in  Red  losses — nearly  all  Red  air¬ 
craft  that  are  engaged  are  lost,  nearly  all  the  time — so  the  significant  dispersion  in 
each  day  is  the  dispersion  in  Blue  losses.  But,  again  by  virtue  of  Blue’s  much 
greater  strength,  the  dispersions  never  lead  to  significant  probabilities  of  Blue’s 
losing  substantial  fractions  of  their  forces,  except  near  the  ends  of  campaigns  in 
which  Blue  is  losing.  The  present  dispersion  calculations  are  useful  for  checking 
the  validity  of  using  only  mean  values.  We  intend  to  treat  the  impact  of  dispersion 
of  day-to-day  losses  in  future  work. 

SLAACM’s  Sequence  of  Calculations 


This  section  describes  SLAACM’s  step-by-step  calculations.  It  explains  how 
SLAACM  uses  the  methods  described  above  to  model  air-to-air  campaigns. 

User  Inputs 

The  SLAACM  user  inputs  Red  and  Blue  orders  of  battle  for  the  first  day  of  fight¬ 
ing  and  lists  any  scheduled  reinforcements. 1  The  following  subsections  summa¬ 
rize  the  inputs,  focusing  on  their  meanings  and  use.  (Chapter  6  details  the 
methods  for  entering  data  into  SLAACM.) 

Inputs  for  the  Blue  Force 

The  user  may  identify  certain  Blue  aircraft  as  smart.  Smart  Blue  aircraft  can  de¬ 
termine  the  composition  of  Red  attack  packages  before  intercepting  them,  and  can 
coordinate  their  defensive  response  among  themselves  to  provide  optimal  de¬ 
fenses. 

Blue  aircraft  based  long  distances  from  the  battle  space  will  have  to  be  rationed 
and  rotated  to  maintain  continuous  defensive  coverage.  The  user  may  select  a 
combat  air  patrol  (CAP)  factor  for  any  Blue  aircraft  type,  so  that  only  a  specified 
fraction  of  the  inventory  of  that  Blue  type  is  available  to  meet  a  Red  attack. 

The  ability  of  Blue  aircraft  that  are  not  smart  to  find  and  reach  the  Red  packages 
is  detennined  by  the  battle  management  efficiency  factor,  b,  discussed  in  Chapter 
4.  This  factor  has  a  user-defined  value  of  0  <  b  <  1,  where  b  corresponds  to  the 
fraction  of  the  assigned  Blue  flights  able  to  intercept  Red  packages. 


1  “Day”  is  used  to  signify  a  set  of  attack  and  defense  sorties.  There  could  be  more  than  one  set 
of  sorties  on  a  single  calendar  day,  or  there  could  be  gaps  of  calendar  days  between  attacks. 


5-2 


How  SLAACM  Works:  SLAACM  Campaign  Calculations 


The  user  may  identify  certain  Blue  aircraft  as  local  air  defenders  (LADs). 
SLAACM  treats  LADs  differently  from  other  Blue  types  when  considering  Blue 
dispatch  options.  Normally  Blue  values  the  life  of  it  forces,  and  will  not  dispatch  a 
Blue  type  if  engaging  a  Red  package  gives  a  negative  value  of  Blue’s  payoff  func¬ 
tion.  Blue  types  designated  as  LADs,  however,  will  be  dispatched  in  all  cases. 

Inputs  for  the  Red  Force 

The  user  inputs  the  bomb  payload  of  each  Red  aircraft  and  specifies  whether  the 
aircraft  is  to  be  used  as  a  fighter  only,  or  as  a  bomber  only,  or  as  either  a  fighter  or 
a  bomber.  The  user  also  inputs  the  ratio  of  the  circular  error  probable  (CEP)  with 
which  bomber  aircraft  can  deliver  bombs,  to  the  reference  CEP.  As  discussed  in 
Chapter  3,  the  default  ratio  is  1.0  for  dumb  bombs  and  less  than  1.0  for  smart 
bombs.  If  delivery  of  cruise  missiles  is  of  interest,  the  user  inputs  the  number  of 
cruise  missiles  that  can  potentially  be  carried  by  each  aircraft. 

Red’s  order  of  battle  may  be  so  large  that  airfield  capacity  and  battle  management 
capability  may  preclude  using  all  forces.  The  user  may  specify  a  maximum  num¬ 
ber  of  Red  aircraft  that  can  be  dispatched  in  a  given  attack. 

Red  optimizes  its  dispatch  choices  with  the  assumption  that  each  package  encoun¬ 
ters  a  specified  Blue  fighter  type.  The  user  identifies  this  type. 

Input  for  Relative  Capabilities  of  Red  and  Blue  Aircraft 

Currently,  SLAACM  models  engagements  with  the  M  vs.  N  probabilistic  en¬ 
gagement  model,  and  it  assumes  that  outcome  probabilities  have  their  long-time 
limiting  values  as  described  in  Chapter  2  under  “Calculations  for  the  Long-Time 
Limit.”  Thus,  only  one  aircraft  parameter — the  ratio  k  of  the  parameter  of  the 
Blue  time-between-kill  distribution  to  that  parameter  for  the  Red  time-between- 
kill  distribution — affects  outcomes. 

Values  of  k  parameters  may  not  be  readily  available,  and  SLAACM  includes  a 
utility  to  calculate  k  from  the  Blue  to  Red  “loss  ratio,”  that  is  the  generally  avail¬ 
able  output  from  simulations.  (Some  simulations  provide  exchange  ratios,  which 
are  ratios  of  Blue  relative  losses  to  Red  relative  losses.  Multiplying  an  exchange 
ratio  by  the  ratio  of  the  initial  number  of  Blues  to  the  initial  number  of  Reds  gives 
the  loss  ratio.)  The  user  inputs  the  loss  ratios  and  the  source  scenario,  such  as 
4  Blues  vs.  8  Reds  with  Blue  breaking  away  after  2  losses.  SLAACM  then  deter¬ 
mines  the  value  of  k  that  gives  the  observed  loss  ratio  and  uses  it  in  subsequent 
calculations.  SLAACM  uses  bisection  to  find  k  values  correct  to  four  significant 
figures. 

With  all  these  inputs  completed,  the  user  enters  the  number  of  days  to  run  the 
campaign,  specifies  the  battle  management  efficiency  parameter  as  described  in 
Chapter  4  under  “The  Battle  Management  Efficiency  Parameter”  and  instructs 
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SLAACM  to  use  either  heuristic  optimization  or  exact  optimization  with  LINGO. 
SLAACM  then  performs  the  sequence  of  calculations  described  below. 

Steps  in  SLAACM’s  Calculations 

SLAACM  first  calculates  values  of  the  Red  payoff  function,  as  described  in 
Chapter  4  under  “Red’s  Payoff  Function,”  for  all  possible  choices  of  attack  escort, 
close  escort,  and  bomber  that  can  make  up  a  Red  assault  package.  There  will  often 
be  hundreds  of  such  choices.  SLAACM  stores  the  payoff  values  for  use  in  later 
optimizations. 

SLAACM  then  begins  the  day-by-day  calculations  of  outcomes  in  the  campaign. 
Using  either  the  heuristic  or  the  exact  solution,  as  specified  by  the  user,  SLAACM 
determines  the  set  of  attack  packages  that  maximizes  the  total  Red  payoff  func¬ 
tion,  subject  to  constraints  imposed  by  Red’s  available  forces  and  to  the  maxi¬ 
mum  number  of  Red  aircraft  that  can  be  dispatched  in  one  day  as  specified  by  the 
user.  These  calculations  are  described  in  Chapter  4  under  “Determining  Attacks 
and  Defenses.” 

SLAACM  considers  this  set  of  attack  packages  to  determine  the  optimal  defense  by 
the  smart  Blue  aircraft.  After  calculating  the  Blue  payoff,  as  described  in  Chapter  4 
under  “Blue  Payoff  Function,”  for  each  possible  combination  of  smart  Blue  defense 
with  a  Red  package  type  actually  dispatched,  SLAACM  determines  the  Blue  dis¬ 
patch  option  that  maximizes  the  total  Blue  payoff,  constrained  by  the  available  Blue 
forces  and  the  CAP  values  input  by  the  user.  SLAACM  uses  either  the  heuristic  or 
the  exact  solution  of  the  resulting  IP  problem,  as  directed  by  the  user. 

As  described  above  under  “Determining  and  Propagating  Loss  Statistics,” 
SLAACM  then  evaluates  the  losses  to  all  forces  from  engagements  between  the 
Red  attack  packages  and  the  smart  Blue  defenders.  Then  SLAACM  calculates  the 
results  of  engagements  among  not-smart  Blue  defenders  and  the  Red  packages 
that  were  not  engaged  by  smart  Blue  aircraft.  Only  those  Red  packages  not  en¬ 
gaged  by  smart  defenders  are  available  for  engagements  with  the  not-smart  Blue 
aircraft. 

SLAACM  calculates  the  Blue  payoff  for  each  combination  of  not-smart  Blue  air¬ 
craft  and  a  Red  attack  package  type  that  was  not  intercepted  by  a  flight  of  smart 
Blue  defenders.  Then,  following  the  procedure  described  in  Chapter  4  under 
“Random  Engagements  by  Not-Smart  Blue  Aircraft,”  and  considering  the  user’s 
input  of  battle  management  efficiency  factor,  SLAACM  determines  the  numbers 
of  engagements  between  the  several  types  of  not-smart  Blues  and  the  Red  pack¬ 
age  types  available  to  them.  As  described  above  under  “Determining  and  Propa¬ 
gating  Loss  Statistics,”  SLAACM  updates  the  losses  of  all  aircraft  types, 
calculates  the  tons  of  “smart”  and  “not-smart”  bombs  delivered,  and,  in  view  of 
any  scheduled  reinforcements,  updates  Red  and  Blue  orders  of  battle  for  the  next 
day  of  the  campaign. 
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After  calculating  outcomes  for  the  number  of  days  specified  by  the  user, 
SLAACM  produces  daily  records  of  Red  and  of  Blue  orders  of  battle,  of  Red 
losses,  of  Blue  losses  and  their  variances  due  to  engagements,  of  tons  of  smart  and 
of  not-smart  bombs  delivered,  and  of  cruise  missiles  delivered.  As  described  in 
Chapter  6,  SLAACM  also  produces  graphs  of  Red  and  of  Blue  drawdowns,  the 
quantities  and  types  of  bombs  delivered,  and  infonnation  on  “who  shot  whom.” 
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Chapter  6 

Users  Guide  to  SLAACM 


SLAACM  is  housed  in  a  Microsoft  Excel  workbook  with  extensive  Visual  Basic 
code.  In  addition,  it  has  an  option  to  conduct  integer  programming  by  linking  to 
an  external  optimization  program  called  LINGO. 1  Using  Excel  to  house 
SLAACM  has  several  advantages.  Most  analysts  are  familiar  with  the  Excel  envi¬ 
ronment  and  will  be  able  to  navigate  it  easily.  Other  than  acquiring  the  SLAACM 
workbook,  users  do  not  need  to  install  additional  software  unless  they  want  to  run 
integer  programming  optimizations  in  LINGO.  Finally,  users  can  easily  copy  and 
paste  output  values  and  charts  into  their  preferred  presentation  formats  for  dis¬ 
playing  results  to  others. 

This  chapter  describes  how  to  operate  SLAACM.  It  begins  with  a  description  of 
the  model  inputs.  It  then  describes  how  to  run  the  model  once  all  of  the  input  pa¬ 
rameters  are  set  up.  Next,  the  chapter  describes  the  outputs  reporting  SLAACM’s 
results.  The  last  section  addresses  how  the  outputs  are  analyzed  to  determine  the 
effects  of  the  various  engagement  scenarios. 


Inputs 


SLAACM’s  input  parameters  are  defined  in  five  worksheets: 

♦  BlueSupply,  which  is  used  to  specify  the  quantities  of  Blue  aircraft  avail¬ 
able  on  each  day  of  the  campaign 

♦  RedSupply,  which  is  used  to  specify  the  quantities  of  Red  aircraft  avail¬ 
able  on  each  day  of  the  campaign 

♦  BlueOOB,  which  is  used  to  define  Blue’s  order  of  battle  (SLAACM  also 
has  a  RedOOB  worksheet,  but  it  is  purely  an  output  worksheet  and  does 
not  receive  any  user  inputs) 

♦  ExRatios,  which  is  used  to  specify  engagement  parameters  such  as  loss  ra¬ 
tios  and  breakpoints 

♦  CM,  which  is  used  to  specify  the  number  of  cruise  missiles  that  a  Red  air¬ 
craft  could  carry. 

We  discuss  each  sheet  below  and  describe  the  proper  way  to  complete  each  one. 


1  For  more  information  on  LINGO,  see  http://www.lindo.com/products/lingo/. 
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BlueSupply  and  RedSupply  Worksheets 

For  each  force,  Blue  and  Red,  the  user  will  need  to  set  up  the  initial  quantities  of 
aircraft  available  and  the  reinforcements  available  on  each  day  of  the  campaign. 
The  user  will  need  to  type  those  quantities  into  the  appropriate  worksheets,  called 
BlueSupply  and  RedSupply. 

Figure  6-1  is  an  example  of  the  BlueSupply  worksheet.  As  shown  in  the  figure, 
the  worksheet  has  a  row  for  each  Blue  aircraft  type  and  a  column  for  each  day  of 
the  campaign.  Initial  aircraft  quantities  are  entered  into  the  column  labeled  0.  Re¬ 
inforcements  may  be  provided  for  subsequent  days  of  the  campaign.  In  the  figure, 
we  see  that  reinforcements  are  entered  for  Day  2  and  Day  3. 2 

Figure  6-1.  BlueSupply  Worksheet 


The  RedSupply  worksheet  is  designed  in  the  same  way  as  the  BlueSupply  sheet. 
An  example  is  shown  by  Figure  6-2.  In  the  figure,  we  see  that  only  initial  aircraft 
quantities  are  provided  and  no  reinforcements  are  specified. 


2  Note  that  a  day  is  used  in  the  model  as  a  convenient  designator  to  represent  an  individual  at¬ 
tack/defense  sortie.  In  an  actual  campaign  there  may  be  more  than  one  sortie  on  a  particular  cal¬ 
endar  day,  or,  conversely,  sorties  may  be  distributed  among  several  calendar  days. 
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Figure  6-2.  RedSupply  Worksheet 


The  BlueSupply  and  RedSupply  worksheets  can  accommodate  up  to  20  aircraft 
types  each.  If  new  rows  are  needed  to  accommodate  additional  aircraft  types, 
they  can  be  added  using  the  “Add  a  New  Row”  button. 


Caution! 

Any  time  a  row  is  added  or  deleted  it  is  imperative  that  the  buttons  on  the 
BlueSupply  and  RedSupply  worksheets  are  used  to  add  or  delete  the  rows. 

Do  not  simply  insert  or  delete  rows  and  columns. 


Several  of  SLAACM ’s  worksheets  reference  the  aircraft  types  and  the  “Add  a 
New  Row”  and  “Delete  a  Row”  buttons  must  be  used  to  keep  the  worksheets  syn¬ 
chronized,  Using  typical  worksheet  approaches  to  inserting  and  deleting  rows  will 
likely  result  in  model  errors  or  misinterpretation  of  results. 

BlueOOB  Worksheet 

The  worksheet  named  BlueOOB  (which  stands  for  Blue  order  of  battle)  is  the 
main  worksheet  for  input  parameters  for  the  Blue  aircraft.  It  is  also  the  worksheet 
used  to  launch  campaign  analyses.  BlueOOB  also  serves  as  an  output  worksheet, 
as  does  the  RedOOB  worksheet;  the  output  portion  of  BlueOOB  will  be  discussed 
in  the  output  section  of  this  chapter.  Figure  6-3  shows  the  BlueOOB  worksheet. 
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Figure  6-3.  BlueOOB  Worksheet 


Caution! 


Rows  or  columns  should  never  be  deleted  from  this  worksheet! 


The  aircraft  listed  in  the  BlueOOB  worksheet  are  linked  to  those  that  were  shown 
in  the  BlueSupply  worksheet  of  Figure  6-1.  The  Day  0  quantities  are  also  linked 
to  those  in  the  BlueSupply  worksheet.  The  user  who  attempts  to  click  the  aircraft 
names  or  Day  0  quantities  on  the  BlueOOB  sheet  will  receive  the  warning  shown 
in  Figure  6-4. 


Figure  6-4.  Warning:  Do  Not  Alter  Aircraft  of  BlueOOB  Worksheet 


I  lot  iditilblr' 


V  Change  Blue  aircraft  quantities  on  BlueSupply  worksheet. 
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The  first  four  columns  on  the  BlueOOB  worksheet  define  aircraft-specific  input 
parameters  for  Blue: 

♦  “Smart?  ”  column.  The  user  may  use  this  column  to  identify  specific  Blue 
aircraft  as  smart.  Smart  Blue  aircraft  can  determine  the  composition  of 
Red  attack  packages  before  intercepting  them  and  can  coordinate  their  de¬ 
fensive  response  among  themselves  to  provide  optimal  defenses.  Coordi¬ 
nation  and  optimization  are  discussed  in  Chapter  5.  Multiple  Blue  aircraft 
types  may  be  designated  as  smart.  In  Figure  6-3,  the  first  two  aircraft  are 
selected  to  be  smart. 

♦  “LAD?  ”  column.  The  user  may  use  this  column  to  identify  certain  Blue 
aircraft  as  local  air  defenders.  SLAACM  treats  LADs  differently  from 
other  Blue  types  when  considering  Blue  dispatch  options.  In  general,  a 
non-LAD  Blue  has  a  penalty  for  its  own  losses  that  is  significantly  greater 
than  a  LAD  Blue,  which  presumably  is  less  loss-averse  due  to  defending 
its  own  territory.  Some  non-LAD  Blues  may  have  a  negative  payoff  func¬ 
tion  for  engaging  certain  Red  packages  dispatched  and  Blue  will  not  dis¬ 
patch  that  type.  LADs,  however,  are  dispatched  in  all  cases.  Multiple  Blue 
aircraft  types  may  be  designated  as  LADs.  In  Figure  6-3,  the  eighth  air¬ 
craft  is  designated  as  such.  Note  that  SLAACM  also  has  a  separate  aircraft 
type  labeled  BlueLAD,  which  is  intended  to  represent  an  indigenous 
“Green”  defender.  In  earlier  versions  of  SLAACM,  this  Blue  LAD  was 
the  only  aircraft  in  the  model  that  placed  no  value  on  its  own  survival.  The 
current  version  allows  all  aircraft  to  be  optionally  “fearless.”  This  facili¬ 
tates  the  analysis  of  supplying  different  types  of  aircraft  to  the  indigenous 
forces.  (We  plan  to  address  any  confusion  in  the  LAD  and  Blue  LAD 
tenninology  in  future  versions  of  SLAACM.) 

♦  “CAP  Factor”  column.  This  factor  addresses  situations  in  which  Blue 
fighters  are  stationed  at  remote  airfields  and  must  divide  their  forces  to 
maintain  continuous  combat  air  patrol  coverage  over  the  battle  space.  The 
user  may  select  a  CAP  factor  for  any  Blue  aircraft  type,  so  that  only  a 
specified  fraction  of  the  inventory  of  that  Blue  type  is  available.  CAP  fac¬ 
tor  values  are  chosen  from  a  drop-down  box  for  each  aircraft  type.  Values 
range  from  1  to  10.  A  CAP  factor  of  1  indicates  that  no  CAP  limitation  is 
placed  on  that  aircraft  type.  In  Figure  6-3,  several  aircraft  have  CAP  fac¬ 
tors  greater  than  1.  For  example,  the  first  aircraft  has  been  assigned  a  CAP 
factor  of  4.  That  aircraft,  BlueFl,  also  begins  with  an  initial  quantity  of 
24  aircraft.  Therefore,  Blue_Fl  has  only  24/4  =  6  available  aircraft  on  the 
first  day,  which  means  that  only  one  4-ship  package  of  BlueFls  can  be 
deployed.  In  addition,  if  BlueFl  losses  on  subsequent  campaign  days 
were  to  exceed  8  aircraft,  Blue  could  no  longer  dispatch  Blue  Fls,  since 
15/4  <  4.  Thus,  Blue  Fl  would  be  unable  to  maintain  a  4-ship  package  at 
all  times.  This  assumption  is  fairly  conservative:  if  there  are  not  enough 
aircraft  to  fulfill  the  combat  air  patrol  at  all  times,  those  flights  are  not  dis¬ 
patched. 
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♦  “Index  of  Blue  a/c  for  Red’s  planning”  column.  Only  one  aircraft  type  can 
be  selected  for  this  parameter.  In  Figure  6-3,  the  first  aircraft  type  is  des¬ 
ignated  as  Red’s  “planning  aircraft.”  This  aircraft  type  is  used  by  Red  to 
compute  its  expected  payoff  function  and  optimize  its  dispatch  choices. 
Experience  has  shown  that  Red  optimization  is  relatively  insensitive  to  the 
selection  of  the  Blue  type,  so  nonnal  practice  is  to  select  the  most  numer¬ 
ous  Blue  type  for  Red  planning. 

In  addition  to  the  aircraft  specific  parameters,  several  other  input  parameters  must 
be  selected  on  the  BlueOOB  worksheet. 

At  the  top  of  the  worksheet,  the  number  of  days  in  the  campaign  is  specified.  In 
Figure  6-3,  the  number  of  campaign  days  is  8.  The  longer  the  campaign,  the 
longer  it  will  take  SLAACM  to  run,  since  each  campaign  day  requires  the  full 
breadth  of  SLAACM  computations.  Initially,  the  user  may  need  to  start  out  with  a 
longer  campaign  and  then  decrease  its  length  as  a  steady  state  in  aircraft  losses  is 
reached.  Eight  is  a  good  initial  number  to  use  since  some  campaigns  will  com¬ 
plete  earlier  and  some  later,  depending  on  the  mixes  and  capabilities  of  Blues  and 
Reds  specified. 

Beneath  the  number  of  campaign  days  is  a  parameter  labeled  “BM  r|.”  This  pa¬ 
rameter  is  the  battle  management  efficiency  factor.  BM  r)  is  a  value  greater  than  0 
and  less  than  or  equal  to  1 .  Simply  stated,  a  BM  r)  less  than  1  indicates  the  frac¬ 
tion  of  time  that  Blue  defender  will  be  able  to  find  a  Red  attack  package.  In  Fig¬ 
ure  6-3,  the  efficiency  factor  is  shown  as  0.85,  which  means  that  15  percent  of  the 
not-smart  Blue  defenders  will  chase  “ghost”  attackers.  No  losses  will  result,  but  in 
those  cases,  Blue  will  not  be  confronting  an  actual  attacker,  so  its  efficiency  is 
reduced.  Since  smart  Blues  know  Red’s  attack  package  compositions  in  advance, 
we  assume  they  have  perfect  management. 

Next  to  the  campaign  days  and  the  efficiency  factor  is  a  drop-down  box  labeled 
“LINGO?”  If  “No”  is  selected,  as  shown  in  Figure  6-3,  SLAACM  will  use  its  in¬ 
ternal  heuristic  optimization  to  calculate  engagement  results.  If  “Yes”  is  selected, 
LINGO  will  be  used  to  perform  integer  optimization  (a  licensed  copy  of  LINGO 
software  is  required  to  be  installed  on  the  PC  for  this  feature  to  work).  Both  the 
heuristic  and  the  exact  integer  programming  solution  are  described  in  detail  in 
Chapter  5.  If  LINGO  is  used,  there  are  two  LINGO  input  worksheets:  “RedOpt” 
and  “BlueOpt.”  However,  those  sheets  are  already  complete  and  should  not  be 
altered  by  the  user. 

The  final  input  selection  on  BlueOOB  is  a  drop-down  box  labeled  “Alerts.”  The 
default  setting  for  this  feature  is  “Off.”  However,  if  “On”  is  selected,  message 
boxes  will  become  visible  as  the  analysis  runs. 
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ExRatios  Worksheet 

The  ExRatios  worksheet  contains  a  row  for  each  Red  aircraft  that  has  been  desig¬ 
nated  on  the  RedSupply  worksheet.  The  worksheet  has  several  parameter  col¬ 
umns,  discussed  below,  and  a  column  for  each  Blue  aircraft  labeled  “F  v.  Blue  x” 
where  “F”  indicates  that  Red  is  a  fighter  and  “x”  corresponds  to  the  label  provided 
on  BlueSupply.  Columns  are  repeated  for  each  Blue  aircraft,  except  they  are  la¬ 
beled  “B  v.  Blue_x”  where  “B”  indicates  that  Red  is  a  bomber  and  “x”  is  as  just 
described.  The  default  setting  in  SLAACM  assumes  that  loss  ratios  are  for  en¬ 
gagements  between  4  Blues  and  8  Reds,  with  the  Reds  and  Blues  fighting  to  anni¬ 
hilation.  However,  those  assumptions  can  be  changed  on  this  sheet  under  “Red 
Start,”  “Blue  Start,”  “Red  Quit,”  and  “Blue  Quit.”  Red  Start  and  Blue  Start  corre¬ 
spond  to  the  engagement  type,  e.g.,  8  and  4,  respectively.  Red  Quit  and  Blue  Quit 
denote  the  breaking  points,  e.g.,  0  and  0  means  that  both  sides  will  fight  to  annihi¬ 
lation.  The  user  will  want  to  specify  values  that  correspond  to  the  conditions  un¬ 
der  which  the  loss  ratios  were  obtained. 

To  keep  the  rows  and  columns  synchronized  with  the  BlueSupply  and  RedSupply 
worksheets,  it  is  imperative  that  no  rows  or  columns  on  ExRatios  be  inserted  or 
deleted. 


Caution! 

As  stated  earlier,  use  the  '  Add  a  New  Row”  and  “ Delete  a  Row”  buttons 
on  the  BlueSupply  and  RedSupply  worksheets  to  add  or  delete  rows. 


Figure  6-5  is  a  snapshot  of  the  ExRatios  worksheet. 
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Figure  6-5.  ExRatios  Worksheet 


The  first  column  on  the  ExRatios  worksheet  is  labeled  “Usage.”  For  each  Red  air¬ 
craft  type,  a  value  of  1,  -1,  or  0  needs  to  be  entered:  1  indicates  that  the  Red  air¬ 
craft  is  a  fighter,  -1  indicates  that  the  Red  aircraft  is  a  bomber,  and  0  indicates 
that  the  Red  aircraft  can  be  used  as  both  a  fighter  and  a  bomber.  The  second  col¬ 
umn  indicates  the  name  of  the  Red  aircraft  type;  it  is  linked  to  the  RedSupply 
worksheet.  The  third  column,  “Payload,  lb,”  is  used  to  enter  the  number  of  pounds 
of  bombs  that  the  associated  Red  aircraft  can  carry. 

The  fourth  column  is  labeled  “cepr,”  which  stands  for  the  circular  error  probable 
ratio  (CEP/CEPref).  The  cepr  is  a  value  greater  than  0  and  less  than  or  equal  to  1, 
where  1  indicates  that  the  bomber  cannot  deliver  smart  bombs  (guided  munitions). 
Values  of  cepr  less  than  1.0  indicate  increasing  accuracy  of  precision  guidance. 
Bomber  payoff  values  used  by  both  Red  and  Blue  are  calculated  from  the  payload 
tonnage  and  the  weapon  accuracy  as  defined  by  the  cepr  value.  The  tonnage  de¬ 
fined  as  smart  bombs  in  the  SLAACM  plotting  routine  is  based  on  a  threshold 
value  contained  on  the  worksheet  named  “RedBombs.” 

The  remaining  “F  v.  Bluex”  and  “B  v.  Bluex”  columns  are  filled  with  loss  ratio 
values.  Those  values  are  generally  determined  from  experience  or  simulation,  and 
are  beyond  the  scope  of  a  user’s  guide.  It  is  expected  that  users  will  have  access  to 
relevant  loss  ratio  data.  Once  the  loss  ratios  are  completed,  the  large  button  at  the 
bottom  of  the  screen  labeled  “Get  Kill  Rate  Ratios”  needs  to  be  clicked.  That  but¬ 
ton  will  run  a  macro  that  uses  the  loss  ratios  to  compute  kill  rate  ratios.  The  kill 
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rate  ratios  are  automatically  entered  onto  the  worksheet  labeled  “Rdata.”  As  an 
added  safeguard,  every  time  SLAACM  is  run,  the  kill  rate  ratio  macro  is  run  to 
ensure  that  the  most  recent  values  are  entered  onto  Rdata. 


SLAACM  will  not  run  if  the  ExRatios  worksheet  contains  any  blank  values.  If  a 
loss  ratio  is  left  blank  on  ExRatios,  the  error  message  shown  in  Figure  6-6  will 
appear. 


Figure  6-6.  Loss  Ratio  Error 


After  clicking  “OK,”  the  user  will  be  taken  to  the  ExRatios  worksheet  and  any 
blank  values  will  be  highlighted  as  shown  in  Figure  6-7. 


Figure  6-7.  Blank  Loss  Ratio  Flighlighted 


CM  Worksheet 


The  final  user  inputs  are  on  the  worksheet  named  “CM,”  which  stands  for  cruise 
missiles.  On  this  worksheet,  the  user  can  specify  the  number  of  cruise  missiles  a 
Red  aircraft  could  carry.  Figure  6-8  shows  that  column  C  is  used  to  input  the 
cruise  missile  capacity  of  each  aircraft. 
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Figure  6-8.  CM  Worksheet 


L^]  Microsoft  E/cs!  • 

-  Arftiyo 

LAACM_ 

Oct'/OO 

C.Az  L 

File  Edit  View 

i  help  Adobe  PDF 

'  Insert 

Format 

X°°ls 

Data 

Window 

_  &  x 

Z  i  B  L 

E  ‘  E 

Aliy  i 

J  1  ' 

- 1  M  750/0  - 

1 

C4 

▼ 

f* 

A 

c— 1 

: eh 

Fl 

G  j— | 

b 

DAY: 

0 

1 

2 

3 

B 

3 

AC  TYPE 

CM  Paili 

sad 

4 

Red  FI 

i 

0 

0 

0 

0 

5 

Red  F2 

1 

0 

0 

0 

0 

6 

Red  F3 

. 1 . 

0 

0 

0 

0 

7 

Red  F5 

1 

0 

0 

0 

0 

8 

Red  F6 

1 

0 

0 

0 

0 

9 

Red  B1 

1 

27 

19 

11 

2 

10 

Red  B2 

1 

51 

38 

21 

3 

11 

Red  B3 

1 

18 

12 

7 

. j.n 

4^ 

- aoo 

_ « _ 

_ fij 

A 

A 

14  < 

i  ►  h  \  CM  /  RedOpt  /  RedIP  / 

BlueC  |£Si|] 

0 1 

|  Ready 

NUM 

No  tradeoff  between  cruise  missiles  and  bomb  payload  is  required  because  the 
cruise  missile  computation  is  a  secondary  analysis  that  is  meant  to  provide  “what- 
if  ’  information  to  the  user.  In  this  computation,  the  number  of  Red  aircraft  or  each 
type  that  get  past  Blue  defenses  each  day  is  multiplied  by  the  number  of  cruise 
missiles  specified  on  the  CM  worksheet.  This  analysis  merely  gives  an  idea  of  the 
potential  number  of  cruise  missiles  that  Red  could  launch  if  its  aircraft  were 
armed  to  do  so. 

Running  SLAACM 


Once  all  of  the  inputs  have  been  entered,  the  user  needs  to  return  to  the  BlueOOB 
worksheet.  To  launch  the  analysis,  the  user  clicks  the  “Run  Campaign”  button. 

If  the  heuristic  is  being  used,  no  other  action  is  required  on  the  part  of  the  user. 

If  LINGO  is  being  used,  then  a  series  of  message  boxes  will  appear  to  let  the  user 
know  what  day  is  being  considered  and  which  optimization  is  taking  place.  The 
user  simply  needs  to  click  “OK”  to  proceed.  The  optimization  in  LINGO  usually 
takes  only  a  couple  of  seconds  and  the  screen  often  clears  automatically.  If  it  does 
not  clear  by  itself,  the  user  may  need  to  click  “Close,”  provided  LINGO  has 
solved  for  the  globally  optimal  solution  (noted  by  “Global  Opt”  in  the  “State”). 
This  case  is  shown  in  Figure  6-9. 
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Figure  6-9.  LINGO  Solution  for  Red  Optimization 


If  the  “State”  shows  “Feasible”  as  its  status,  then  the  user  should  wait  until  the 
global  optimal  is  reached.  If  the  solver  is  running  for  a  long  time,  and  the  differ¬ 
ence  between  the  “Best  Obj”  and  “Obj  Bound”  is  small  (e.g.,  less  than  0.5%), 
then  the  user  may  click  “Interrupt  Solver”  and  the  best  feasible  solution  at  that 
point  will  be  entered  into  SLAACM.  This  event  occurs  rarely  in  SLAACM,  since 
LINGO  is  usually  able  to  solve  these  integer  programming  problems  very  quickly. 

Figure  6-10  shows  LINGO  before  it  has  reached  optimality;  it  has  reached  a  fea¬ 
sible — but  not  necessarily  optimal — solution.  In  this  example,  the  difference  be¬ 
tween  the  best  objective  and  the  objective  bound  is  less  than  1/1,712  =  0.05%,  so 
one  could  certainly  interrupt  the  solver  at  this  point. 
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Figure  6-10.  Feasible  LINGO  Solution 


Outputs 


Once  the  analysis  has  been  completed,  the  user  has  several  output  results  to  re¬ 
view.  SLAACM  produces  daily  records  of  Red  and  Blue  orders  of  battle,  of  Red 
and  Blue  losses,  of  their  variances  due  to  engagements,  of  tons  of  smart  and  not- 
smart  bombs  delivered,  and  of  cruise  missiles  delivered. 

Six  worksheets  report  SLAACM ’s  output  results: 

♦  BlueOOB 

♦  RedOOB 

♦  RedBombs 

♦  Losses 

♦  MBCharts 

♦  CM. 

We  discuss  the  outputs — charts  and  worksheets  tables — below.  By  reviewing  the 
charts  and  tables,  the  user  can  gain  an  understanding  of  specific  campaign  out¬ 
comes. 


6-12 


Users  Guide  to  SLAACM 


Charts 


The  main  results  are  presented  graphically.  From  the  BlueOOB  worksheet,  the 
user  can  click  the  “View  Charts”  button,  or  simply  proceed  to  the  worksheet  tab 
labeled  “MB  Charts.”  Figure  6-11  shows  the  quad  chart  “dashboard”  of  the  main 
SLAACM  results. 


Figure  6-11.  Output  Charts 
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The  “Red  Forces”  chart  in  the  upper  left  comer  of  Figure  6-11  shows  the  quantity 
of  Red  aircraft  on  each  day  of  the  campaign.  In  the  example,  most  Red  aircraft  are 
annihilated  by  Day  4.  The  aircraft  that  have  some  quantity  remaining,  but  which 
reach  steady-state  quantities,  are  most  likely  not  being  sent.  For  example,  Red  FI 
is  a  weak  fighter,  so  Red  is  choosing  to  escort  its  bombers  with  more  capable  air¬ 
craft.  A  quick  check  of  the  “Tons  of  Bombs  Dropped”  chart  in  the  lower  right  cor¬ 
ner  shows  that  it  looks  like  the  bomb  levels  have  stopped  growing  by  Day  4, 
indicating  that  no  Red  bombers  are  getting  through  Blue  defenses  beyond  that 
day.  That  leads  us  to  believe  that  the  campaign  is  over  by  Day  4. 

The  “Blue  Forces”  chart  in  the  upper  right  corner  of  Figure  6-11  shows  the  daily 
quantities  of  Blue  aircraft.  We  see  that  Blue  sustains  some  losses,  although  not 
nearly  as  many  as  Red;  therefore,  we  can  assume  that  Blue’s  forces  are  superior  to 
Red’s.  In  addition,  in  the  first  few  days,  we  see  increases  in  some  quantities  of 
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Blue  aircraft.  That  is  due  to  Blue  receiving  reinforcements.  (See  the  BlueSupply 
worksheet  in  Figure  6-1.)  Blue  LADs  are  broken  out  separately  from  the  Blue 
Forces  and  are  shown  in  the  chart  in  the  lower  left  corner  of  Figure  6-11. 

The  Red  Forces,  Blue  Forces,  and  LAD  charts  show  the  quantity  of  each  aircraft 
type  that  is  available  on  each  day  of  the  campaign.  The  charts  do  not  indicate  the 
number  of  aircraft  sent  into  battle. 

Figure  6-12  shows  a  graph  of  the  potential  number  of  cruise  missile  launches.  As 
discussed  in  the  previous  section,  this  result  assumes  that  the  Red  bombers  that 
got  through  Blue  defenses  were  armed  with  cruise  missiles  in  lieu  of  bombs  and 
were  able  to  launch  their  cruise  missiles. 

Figure  6-12.  Cruise  Missile  Launches 


Worksheet  Tables 

Because  SLAACM  is  housed  in  an  Excel  workbook,  the  user  can  browse  the  data 
used  to  build  the  charts.  All  of  the  charts  in  Figure  6-11  are  built  from  data  in 
worksheet  tables: 

♦  The  Red  Forces  chart  uses  the  data  on  the  “RedOOB”  worksheet. 

♦  The  Blue  Forces  and  LAD  charts  use  data  from  the  “BlueOOB”  work¬ 
sheet.  (As  we  discussed  earlier,  the  BlueOOB  worksheet  is  both  an  input 
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and  an  output.  After  running  SLAACM,  available  aircraft  quantities  are 
loaded  into  BlueOOB  worksheet  for  each  day  of  the  campaign.) 

♦  The  Tons  of  Bombs  chart  uses  data  from  the  “RedBombs”  worksheet. 

♦  The  Cruise  Missile  chart  is  built  from  data  on  the  “CM”  worksheet,  where 
the  chart  resides. 

It  can  be  helpful  to  look  at  aircraft  losses  to  get  a  better  idea  of  how  the  campaign 
progressed.  SLAACM  has  a  worksheet  named  “Losses”  that  shows  the  loss  data. 
Figure  6-13  shows  the  Losses  worksheet.  Blue  and  Red  losses  are  shown  by  air¬ 
craft  type.  In  addition,  the  worksheet  shows  the  standard  deviation  of  the  losses 
for  Blue  aircraft. 


Figure  6-13.  Losses 
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Figure  6-14  shows  the  pounds  and  tons  of  bombs  dropped  by  Red  aircraft  each 
day  of  the  campaign.  In  addition,  the  bomb  types  are  broken  out  by  smart  and  not- 
smart  types.  These  are  the  data  used  to  build  the  “Tons  of  Bombs  Dropped”  chart 
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in  the  lower  right  comer  of  Figure  6-11.  Upon  inspection  of  the  data  in  Figure  6- 
14,  we  see  that  the  campaign  actually  ends  on  Day  5,  not  Day  4  as  it  had  appeared 
on  the  chart;  one  can  see  that  2  tons  of  bombs  were  dropped  on  Day  5. 

Figure  6-14.  Red  Bombs 


User  Note 

The  output  worksheets  in  SLAACM  are  dynamic  and  are  overwritten  dur¬ 
ing  each  model  run;  users  who  want  to  save  run  data  should  copy  and 
paste  the  values  into  an  external  file. 


Analyses  of  Sample  Runs 


In  this  section,  we  conduct  a  comparative  analysis  of  the  heuristic  and  LINGO 
using  two  SLAACM  runs  to  illustrate  some  of  the  subtle  differences  that  occur 
under  reasonably  similar  scenarios. 

First,  we  run  SLAACM  using  the  inputs  shown  in  Figures  6-1  through  6-5  using 
the  heuristic;  we  then  repeat  the  run  using  LINGO.  As  shown  in  Table  6-1,  the 
results  are  identical;  the  heuristic  performs  the  same  as  the  LINGO  optimization. 
This  case  is  not  uncommon,  and  it  gives  us  confidence  that  the  heuristic  can  be  a 
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good  approximation  for  the  optimal  solution  (the  heuristic  and  the  exact  integer 
programming  formulations  are  covered  in  detail  in  Chapter  4). 

Table  6-1.  Results  of  Initial  SLAACM  Run 


Item 

Heuristic 

LINGO 

Tons  dropped 

Smart  bombs 

i  1,241 

1,241 

Dumb  bombs 

174 

174 

Aircraft  losses 

Blue 

19 

19 

Red 

397 

397 

To  illustrate  a  situation  in  which  the  heuristic  and  LINGO  solutions  differ,  con¬ 
sider  the  same  scenario,  except  that  all  Blue  aircraft  are  smart.  This  change  re¬ 
quires  checking  the  boxes  on  the  BlueOOB  worksheet  (Figure  6-3).  The  results 
are  shown  in  Table  6-2.  For  both  model  runs,  the  tonnage  of  smart  bombs  dropped 
is  reduced  by  almost  50  percent,  demonstrating  the  impact  of  smart  Blues  prefer¬ 
entially  engaging  the  highest  value  Red  packages.  Table  6-2  also  shows  that  the 
heuristic  results  and  LINGO’s  exact  IP  results  show  significant  differences. 


Table  6-2.  Results  of  Setting  All  Blues  to  Be  Smart 


Item 

Heuristic 

LINGO 

Tons  dropped 

Smart  bombs 

640 

602 

Dumb  bombs 

289 

289 

Aircraft  losses 

Blue 

19 

17 

Red 

402 

404 

LINGO’s  optimal  solution  resulted  in  38  fewer  tons  of  smart  bombs  being 
dropped,  2  fewer  Blue  aircraft  losses,  and  2  additional  Red  aircraft  losses.  In 
short,  the  IP  solution  using  LINGO  resulted  in  Blue  reducing  the  amount  of  dam¬ 
age  while  sustaining  fewer  losses.  This  example  illustrates  what  we  have  seen  in 
other  results,  as  well  as  what  one  might  expect.  The  heuristic  and  the  exact  solu¬ 
tion  to  the  IP  problem  in  LINGO  generally  produce  qualitatively  the  same  results. 
However,  if  one  is  looking  for  very  precise  results,  the  heuristic  may  fall  short  of 
true  optimality  by  several  percentage  points. 
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Summary 


SLAACM  is  a  fast,  flexible,  robust  model  that  contains  the  key  input  parameters 
necessary  to  define  the  characteristics  of  a  realistic  air-to-air  campaign.  Model 
results  in  both  tabular  and  graphical  format  clearly  display  the  impacts  of  parame¬ 
ter  choices  and  provide  insight  into  campaign  scenarios.  The  Excel  workbook  im¬ 
plementation  allows  results  to  be  easily  copied  and  pasted  to  other  applications  for 
reporting  and  presentation. 

The  above  examples  show  differences  in  results  using  the  heuristic  optimization 
and  solving  the  exact  IP  formulation  with  LINGO.  In  the  many  practical  cases  we 
have  analyzed,  it  has  been  rare  for  the  heuristic  results  to  differ  from  the  LINGO 
results  by  more  than  10  percent,  but  careful  analyses  may  require  the  use  of  exact 
solution  methods,  especially  for  actual  campaign  planning. 
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Appendix  A 

Direct  Computation  of  Long-Time  Limiting 
Probabilities  of  Boundary  States 


The  structure  of  the  evolution  equations  for  probability  distributions  of  the  en¬ 
gagements  we  consider  makes  it  possible  to  find  the  long-tenn  limits  of  the  ab¬ 
sorbing  boundary  states  directly,  in  finitely  many  steps,  without  formally  solving 
those  differential  equations.  In  some  interesting  cases,  the  finite  steps  may  be  car¬ 
ried  out  by  straightforward  iteration.  This  appendix  explains  these  facts  and  gives 
some  examples. 

The  evolution  equations  for  the  state  probabilities  in  our  engagement  models  de¬ 
compose  naturally  into  a  set  of  equations  for  the  transient  states  and  a  set  for  the 
absorbing  boundary  states.  Equations  2-5  and  2-6  are  examples  of  this  for  the  ba¬ 
sic  M  vs.  N  probabilistic  engagement  model.  The  advanced  engagement  models 
of  Appendix  D  also  have  this  property,  except  for  the  low-observable  model, 
which  is  treated  completely  by  the  method  given  in  Appendix  D. 

If  we  make  the  probabilities  of  the  transient  states  the  components  of  a  vector  x, 
and  the  components  of  the  absorbing  boundary  states  into  the  components  of  a 
vector  y,  the  evolution  of  the  state  probabilities  can  be  described  by 


x  =  Ax  ,  x(0)  =  x0 
y  =  Bx 


[Eq.  A-l] 


To  illustrate  these  ideas  for  an  example  that,  unlike  the  one  of  Chapter  2,  has  re¬ 
current  states,  let  us  consider  treating  a  1  vs.  1  engagement  with  a  two-phase  kill 
model  with  break-lock.  We  use  the  state  description  (m,  i;  n,  j)  where  m  denotes 
the  number  of  Blue  aircraft,  i  the  number  of  Blues  tracking  opponents,  n  the  num¬ 
ber  of  Red  aircraft,  and  j  the  number  of  Reds  tracking  opponents.  Figure  A-l  is  a 
diagram  of  the  states  in  this  engagement  and  their  transitions. 
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Figure  A-1.  Diagram  of  1  vs.  1  Engagement 


The  evolution  equations  are 

P1010  =  — ' (kbj  +  krj)P1010  +bxP1011  +rxPino 
Pmo  =  — (kb2  +  kq  +rx)Piiio  +  kbjPjoj,,  +bxPnn 
Pion  =  — (kr,  +kbj  +bx)P1011  +  kr1P1010  +  rxPull 
Pnn  =  —  (kb2  +  kr2  +  rx  +bx)Pnn  +krjP1110  +kb1Plon 
Piooo  =  kb2(P1110  +Pnn) 

Pooio  =  kr2(pioii  +P  mi) 


[Eq.  A-2] 


The  first  four  of  these  describe  the  transient  states,  and  the  last  two  describe  the 
absorbing  boundary  states.  Defining  xi  as  Pioio,  X2  as  Pmo,  X3  as  P1011,  and  x4  as 
Pun,  and  defining  yi  as  Piooo,  yi  as  Pooio,  we  see  that,  in  this  case,  the  matrices  A 
and  B  of  Equation  A-1  are 


'-(kbj+krO 

kbj 

kr, 

rx 

-(kb2  +kr!  +rx) 

0 

bx 

0 

-(kr,  +kb!  +bx) 

0  ' 

bx 

rx 

[Eq.  A-3] 

l  0 

kfi 

kbj 

- 

(kr2+kb,+rx+bx); 

and 

B  = 

^0  kb,  0 
v0  0  kr2 

kb2  A 

kr2  y 

[Eq.  A-4] 
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Now,  the  solutions  of  systems  of  ordinary  differential  equations  of  the  fonn 
x  =  Ax  for  constant  coefficient  matrices  A  can  always  be  expressed  as  finite  lin¬ 
ear  combinations  of  generalized  exponential  functions  of  the  time,  that  is,  func¬ 
tions  of  the  fonn  tJe/a,  where  the  j  are  positive  integers  and  the  X  are  the 
eigenvalues  of  the  matrix  A. 1  If  the  transients  are,  in  fact,  transient,  then  the  real 
parts  of  the  X  are  all  strictly  less  than  0.  It  follows  that  the  x; ,  defined  by 


Xjs|xj(t)dt,  [Eq.A-5] 

o 

exist.  Then,  on  integrating  the  differential  equations  of  Equation  A-l  from  0  to  go, 
using  the  initial  condition  on  x  and  y,  and  remembering  that  the  xft)  tend  to  zero 
as  t  — >  go,  we  find 


-  x0  =  Ax 

y  lim  =  B*  ’ 


[Eq.  A-6] 


where  yum  denotes  the  limit  of  y  as  t  — »  oo.  It  follows  that 

y  Km  =  -BA~ 1  xo  •  [Eq-  a-7] 


Equation  A-7  expresses  the  long-time  limiting  values  of  the  absorbing  boundary 
state  probabilities  as  the  result  of  finite  operations,  that  is,  matrix  inversion  and 
multiplication. 

Let  us  continue  to  illustrate  these  concepts  with  the  1  vs.  1  example  considered 
above.  For  a  numerical  example,  we  take  kri  =  1,  kr2  =  2,  rx  =  1,  kbi  =  3,  kbi  =  4, 
and  bx  =  2.  Then  the  matrices  A  and  B  take  the  values 


and 


'-4  1  2 

3-6  0 

1  0 
0  1  3 


0  N 
2 
1 


B  = 


"0  4 

v°  0 


0  4^ 
2  2y 


[Eq.  A-8] 


[Eq.  A-9] 


It  is  as  well  to  check  that  the  eigenvalues  of  A  are  distinct  and  negative.  Direct 
calculation  shows  that  this  is  in  fact  the  case:  the  eigenvalues  are  approximately 


1  W.E.  Boyce  and  R.C.  DiPrima,  Elementary’  Differential  Equations ,  Third  Edition  (New 
York:  Wiley,  1976),  Sections  7-6  through  7-9. 
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-2.44,  -6.44,  -6.56,  and  -10.56.  Assured  by  this  that  the  x;  exist,  we  go  on  to 
find 


y  lim 


-BA_1x 


o  “ 


"0.827^ 

v0.173y 


[Eq.  A- 10] 


That  is,  the  probability  that  the  Blue  aircraft  defeats  the  Red  one  is  roughly 
83  percent. 

The  eigenvalues  of  the  plant  matrix  A  have  another  use:  they  tell  one  whether  or 
not  it  makes  sense  to  look  for  the  long-time  limiting  values  of  the  absorbing 
boundary-state  probabilities.  Those  probabilities  are  of  interest  if,  but  only  if,  an 
actual  engagement  can  continue  long  enough  for  the  system  to  be  in  boundary 
states  with  a  probability  near  one.  The  time  for  this  to  happen  is  the  time  for 
which  the  transient  states’  probabilities  are  all  much  less  than  one. 

When  the  eigenvalues  of  A  are  distinct  (and,  of  course,  have  negative  real  parts) 
so  that  the  x,(t)  are  linear  combinations  of  the  functions  cxp(A,t),  that  time  will  be 
no  greater  than  the  time  at  which  exp(A*t)  becomes  much  less  than  one,  where  A* 
is  the  eigenvalue  closest  to  zero.  The  time  to  make  that  happen  is  approximately 
2/A,*.  If  the  rate  parameters  kbi,  ...  of  our  example  are  in  inverse  minutes,  the 
mean  time  for  the  Blue  aircraft  to  acquire  a  target  is  20  seconds,  the  mean  time  for 
them  to  launch  a  successful  missile  is  15  seconds,  and  the  mean  time  for  them  to 
break  a  Red  lock  is  30  seconds,  while  the  mean  time  for  the  Reds  to  make  lock  is 
60  seconds,  the  mean  time  for  them  to  launch  a  successful  missile  is  30  seconds, 
and  the  mean  time  for  them  to  break  lock  is  60  seconds.  The  negative  eigenvalue 
of  the  plant  matrix  A  with  smallest  magnitude  is  -2.44  inverse  minutes,  which 
implies  that  the  time  for  the  transient  phase  of  the  1  vs.  1  engagement  to  be  over  is 
roughly  50  seconds.  Very  likely,  the  combatants  will  have  enough  fuel  to  fight 
that  long,  and  so  the  long-time  limit  is  meaningful  in  our  example. 

Although  we  have  found  numerical  evaluation  of  Equation  A-7  quite  helpful  in 
generating  insight  with  examples  of  modest  dimension,  determining  long-time 
limiting  probabilities  in  this  way  may  not  be  practical  for  systems  of  large  dimen¬ 
sion.  In  an  important  class  of  engagement  models,  however,  the  difficult  task  of 
computing  A"  xo  can  be  done  with  a  straightforward  (although  possibly  lengthy) 
iterative  scheme.  That  class  is  those  engagement  models  whose  diagrams  are 
acyclic,  like  the  example  of  Chapter  2.  For  these  models,  the  plant  matrix  A  is 
lower  triangular,  and  this  provides  the  iterative  scheme.  In  these  cases,  one  has 

anXi  =-l 

a21Xj  +  a22x2  =  0  [Eq.  A- 11] 
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so  that 


Xj  =  -1/ an 

x2  =  — Xj  a 21  /a22  [Eq.  A-12] 


With  the  x  determined,  the  long-time  limiting  values  of  the  absorbing  boundary 
states  follow  from  Equation  A-7. 

Determining  the  set  of  states  actually  occupied  in  a  given  engagement  can  be 
somewhat  tedious.  For  example,  finding  the  long-time  limit  of  a  missile-tracking 
engagement  model  for  4  vs.  4,  when  each  aircraft  has  6  missiles,  involves  more 
than  42,000  states.  Obviously,  it  is  not  practical  to  evaluate  these  cases  by  hand. 
In  these  cases,  we  generally  use  C++  code  to  find  both  the  set  of  states  occupied 
and  the  solution  (Equation  A-12). 
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Appendix  B 

Analyzing  Large  State-Space  Engagement 
Models  with  ASSIST  and  STEM 


Evolution  equations  like  Equations  2-5  and  2-6  of  Chapter  2  are  simple  in  principle. 
They  are  systems  of  linear  ordinary  differential  equations  with  constant  coefficients. 
The  solutions  of  the  initial  value  problems  formed  by  such  equations  and  specifica¬ 
tions  of  the  starting  probabilities  of  the  various  states  can  be  written  symboli¬ 
cally — and  perhaps  treated  practically —  with  matrix  exponentiation,  and  they 
generally  offer  no  difficulty  to  numerical  solution  by  difference  schemes,  but  they 
can  be  of  large  dimension. 

A  basic  4  vs.  4  engagement  has  24  states  and  32  transitions;  adding  6  missiles  to 
the  Blue  aircraft  produces  1,632  states  and  9,087  transitions;  adding  6  missiles  to 
both  Red  and  Blue  produces  108,534  states  and  2,416,252  transitions.  Also,  it  is 
somewhat  awkward  to  detennine  the  long-time  limits  of  the  absorbing  boundary 
states  by  integrating  the  differential  equations  numerically.  This  appendix  dis¬ 
cusses  NASA-developed  tools  that  we  have  found  helpful  in  dealing  with  both 
these  issues. 

The  decision  to  apply  the  NASA  tools  to  the  problem  of  fighter  combat  was  ser¬ 
endipitous.  One  of  us  was  engaged  in  Markov  analysis  of  fighter  combat,  and  an¬ 
other  of  us  was  using  the  NASA  tools  for  safety-related  reliability  analysis.  The 
need  for  a  tool  to  efficiently  develop  relatively  complex  fighter  combat  models 
surfaced  during  informal  discussions.  Using  the  NASA  tools,  we  can  develop  and 
run  complex  models  for  Blue  versus  Red  combat,  including  tracking  of  Blue  mis¬ 
sile  use,  in  relatively  short  order.  These  results  provide  insight  into  the  combat 
problem  and  support  development  of  the  SLAACM  algorithms. 

NASA  has  three  tools  that  can  be  applied  to  fighter  combat.  Two  are  Markov  compu¬ 
tational  analysis  programs.  The  third  is  a  sophisticated  utility  program  that  generates 
the  inputs  (models)  for  the  analysis  programs.  The  tools  were  developed  by  NASA  to 
estimate  failure  probabilities  in  highly  reliable,  reconfigurable  avionics  and  space 
electronics. 

The  two  analysis  programs  are  Scaled  Taylor  Exponential  Matrix  (STEM)  1  and 
Semi-Markov  Unreliability  Range  Estimator  (SURE).2  The  former  is  a  pure 
Markov  analysis  tool  in  which  all  failure  rates  are  constant  with  state  probabilities 


1  NASA,  The  PAWS  and  STEM  Reliability’  Analysis  Programs,  Technical  Memorandum 
100572,  R.  Butler  and  P.  Stevenson,  March  1988. 

2  NASA,  SURE  Reliability’  Analysis,  Program  and  Mathematics,  Technical  Paper  2764, 
R.  Butler  and  A.  White,  1998. 
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having  generalized  exponential  fonn.  The  latter  is  a  semi-Markov  model  that  al¬ 
lows  use  of  non-exponential  reconfiguration  probabilities.  STEM  and  SURE  use 
identical  input  files,  but  some  commands  are  processed  only  by  SURE.  For  rea¬ 
sons  discussed  below,  STEM  is  the  tool  used  for  fighter  combat  calculations.  The 
utility  program  ASSIST  (Abstract  Semi-Markov  Specification  Interface  to  the 
SURE  Tool)3  generates  the  STEM  and  SURE  input  files,  or  “models.”  ASSIST 
allows  the  straightforward  generation  of  extremely  complex  Markov  models.4 

Reliability  analyses  are  typically  conducted  for  short  model  times  compared  to  the 
failure  rates  of  the  components  studied,  such  as  a  10-hour  flight  and  a  failure  rate 
of  0.0001  failure  per  hour.  For  air  combat  models,  we  are  interested  in  running  the 
engagement  to  completion  so  we  will  run  model  time  units  on  the  order  of  the  re¬ 
ciprocal  of  the  lowest  kill  rate,  e.g.,  1,000  units  for  a  kill  rate  of  0.001.  STEM  has 
no  problem  running  such  models,  but  SURE’s  mathematical  algorithms  that  han¬ 
dle  non-Markov  recovery  rates  are  not  designed  to  run  to  completion.  Conse¬ 
quently,  the  remainder  of  the  discussion  will  focus  on  ASSIST  and  STEM. 

All  three  programs  were  developed  for  UNIX  platforms  and  have  been  converted 
to  Windows.  STEM  is  also  available  on  LINUX.  The  programs  are  available  at  no 
charge  from  NASA  Langley  Research  Center.  Documentation  includes  an 
ASSIST  users  guide,  a  report  on  modeling  techniques,  and  reports  on  SURE  and 
STEM  mathematics  and  performance.  Figure  B-l  shows  the  basic  relationships  of 
the  tools. 


Figure  B-1.  Basic  Relationships  of  the  Tools 


HOL:  Higher  Order  Language 

ASSIST:  Abstract  Semi-Markov  Specification  Interface  to  the  SURE  Tool  -  Input  model  generator 

Model:  Designation  for  the  input  files  for  the  Markov  analysis  programs,  containing  variables,  states,  transitions,  and  rates 
STEM:  Scaled  Taylor  Exponential  Matrix  -  Markov  analysis  program 

SURE:  Semi-Markov  Unreliability  Range  Estimator  -  Markov  and  Semi-Markov  analysis  program 

Fighter  Combat  Analysis 


In  this  section,  we  use  case  examples  of  increasing  complexity  to  describe  the  ca¬ 
pabilities  and  limitations  of  the  tools  for  analysis  of  fighter  combat.  All  the  analy- 

3  NASA  Langley  Research  Center.  ASSIST  User  Manual,  S.  Johnson  and  D.  Boerschlein, 
September  1993. 

4  NASA,  Techniques  for  Modeling  the  Reliability  of  Fault-Tolerant  Systems  with  the  Markov 
State-Space  Approach,  Reference  Publication  1348,  R.  Butler  and  S.  Johnson,  September  1995. 
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ses  use  the  common  structure  of  an  M  vs.  N,  Blue  vs.  Red  engagement  with  con¬ 
stant  kill  rates. 

Example  1 — Simple  2  vs.  2  Engagement 

The  2  vs.  2  engagement  demonstrates  the  basic  structure  of  the  models  without 
generating  extensive  output.  The  Blue  fighter  kill  rate  is  0.1  per  unit  time,  and  the 
Red  fighter  kill  rate  is  0.01  per  unit  time.  Listing  1  is  the  ASSIST  code  that  gener¬ 
ates  the  2  vs.  2  model. 

Listing  1.  ASSIST  Code  for  Generating  the  2  vs.  2  Model 

(*  ASSIST  model  for  Fighter  Combat  *) 

LIST  =  3;  (*  1  lists  accumulated  death  state  probabilities  only  *) 

(*  2  lists  all  state  probabilities  *) 

(*  3  lists  transitions  and  all  state  probabilities  *) 

ONEDEATH  OFF;  (*  OFF  enumerates  all  absorbing  (Death)  states  *) 

(*  comment  out  to  consolidate  absorbing  states  *) 

(*  equipment  list  *) 

nblue  =2;  (*  blue  fighters  *) 

nred  =2;  (*  red  fighters  *) 

(*  kill  probabilities  *) 
k_rate_blue  =  0.1; 
k_rate_red  =  0.01; 

SPACE  =  (blue:  0. nblue,  red:  O.nred); 

START  =  (nblue,  nred) ; 

DEATH IF  (blue  =  0)  AND  (red>0); 

DEATH IF  (red  =  0)  AND  (blue>0); 

(*  Blue  kills  *) 

IF  (blue  >  0)and  (red>0)  TRANTO  red  =  red-1  BY  blue*k_rate_blue; 

(*  Red  kills  *) 

IF  (red  >  0  )  and  (blue>0)  TRANTO  blue  =  blue-1  BY  red*k_rate_red; 

ASSIST  uses  a  higher  order  definition  language  that  allows  algebraic  manipula¬ 
tion  of  variables  and  compact  description  of  Markov  state  transfers.5  Several  fea¬ 
tures  are  noteworthy  in  the  listing  above.  The  listing  starts  with  the  editing 
commands  LIST  and  ONEDEATH  that  control  SURE/STEM  output.  Next,  the 
number  of  aircraft  and  their  kill  rates  are  input  parameters  defined  as  constant 
types.  The  SPACE  statement  defines  the  range  of  the  Markov  state  space,  i.e.,  in 
this  case,  the  first  state  can  vary  from  0  to  nblue  where  nblue  equals  2.  The 
START  statement  identifies  the  initial  state  populations.  The  DEATHIF  state¬ 
ments  define  the  absorbing  states  in  the  model;  in  this  case,  the  two  DEATHIF 
conditions  represent  Blue  and  Red  wins.  If  the  ONDEATH  OFF  command  is  en¬ 
abled,  STEM  and  SURE  list  all  the  individual  absorbing  states;  otherwise,  they 
combine  the  results  for  absorbing  states  into  the  appropriate  Death  states.6 


5  ASSIST  inputs  are  fully  defined  in  the  citations  of  Notes  1-3. 

6  SURE  mathematics  requires  the  existence  of  absorbing  states,  while  STEM  can  handle  mod¬ 
els  that  include  recovery  from  all  states. 
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We  now  come  to  the  powerful  TRANTO  statements,  which  define  the  conditions, 
nature,  and  rate  of  Markov  state  transfers.  In  this  example,  the  transfers  occur 
whenever  both  Blue  and  Red  fighters  are  available;  kills  are  made  one  at  a  time, 
and  the  kill  rates  are  proportional  to  the  number  of  killers.  The  two  TRANTO 
statements  above  are  adequate  to  generate  the  models  for  any  combination  of  Red 
and  Blue  aircraft. 

The  state  diagram  for  the  2  vs.  2  engagement,  Figure  B-2,  is  relatively  simple.  It 
shows  that  there  are  8  total  states,  of  which  4  are  absorbing  (Death)  states  and  8 
are  transitions.7 


Figure  B-2.  2  vs.  2  Engagement  Diagram 


Designates  an  intermediate  state  with  X  Blue  fighters  and  Y  Red  fighters 


Designates  an  absorbing  (Death)  state  with  X  Blue  fighters  no  Red  fighters  or  visa  versa 


kB  and 


Designate  Blue  and  Red  kill  rates 


Running  ASSIST  with  the  listing  above  generates  the  STEM/SURE  input  file  or 
“model”  shown  in  Listing  2.  Note  the  8  states,  the  4  absorbing  states  (0,  N  or  M, 
0),  and  the  8  transitions.  The  file  extension  is  “.mod.” 


7  There  is  no  (0,  0)  absorbing  state  because  we  chose  not  to  include  the  case  in  which  the  last 
two  opponents  simultaneously  shoot  each  other  down.  Such  a  state  could  be  modeled  with  an  ad¬ 
ditional  TRANTO  statement,  but  determination  of  the  controlling  kill  rate  would  take  some  care¬ 
ful  thought.  The  pure  statistical  rate  for  independent  events  suggests  that  it  would  be  equal  to 
KRATEBLUE  *  K  RATE  RED,  but  simultaneous  kills  can  happen  only  in  certain  conditions 
such  as  ramming,  head-to-head  gun  attacks,  and  head-to-head  missile  attacks.  These  are  only  a 
subset  of  the  configurations  that  make  up  the  basic  Blue  and  Red  kill  rates.  We  assessed  the  likely 
rate  to  be  sufficiently  low  to  justify  ignoring  the  (0,  0)  state  until  better  data  become  available. 
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Listing  2.  2  vs.  2  Engagement  Model 


LIST  =  3; 

TIME  =  100; 

NBLUE  =  2; 

NRED  =  2; 

MIN_B  =  0; 

MIN_R  =  0; 
K_RATE_BLUE  =  0.1; 
K  RATE  RED  =  0.01; 


1  (* 

2,2 

*) 

2  (* 

t — 1 

C\] 

*) 

* 

C\] 

II 

RATE 

BLUE ; 

1  (* 

2,2 

*) 

3  ( * 

1,2 

*) 

* 

CM 

II 

RATE 

_RED; 

2  (* 

t — 1 

C\] 

*) 

4  (* 

o 

C\] 

*) 

=  2*K 

RATE 

BLUE ; 

2  (* 

2,1 

*) 

5  (* 

1,1 

*) 

=  1*K 

RATE 

_RED; 

3  ( * 

1,2 

*) 

5  (* 

5 - 1 

5 - 1 

*) 

=  1*K 

RATE_ 

BLUE ; 

3  ( * 

1,2 

*) 

6  (* 

C\] 

o 

*) 

* 

CM 

II 

RATE_ 

_RED; 

5  (* 

1,1 

*) 

7  (* 

o 

t — 1 

*) 

* 

5 - 1 

II 

RATE_ 

BLUE ; 

5  (* 

1,1 

*) 

8  (* 

5 - 1 

o 

*) 

* 

5 - 1 

II 

RATE_ 

_RED; 

(*  NUMBER  OF  STATES  IN  MODEL  =  8  *) 

(*  NUMBER  OF  TRANSITIONS  IN  MODEL  =  8  *) 


The  model  file  contains  the  input  definitions  and  transition  descriptions  needed  by 
SURE  and  STEM.  Transitions  are  defined  by  source  state  number,  destination 
state  number,  and  rate  of  transfer  between  the  source  and  destination.  The  state 
descriptions,  e.g.,  (*2,2*),  are  included  as  comments;  “(*”  and  “*)”  are  comment 
delimiters.  Note  that  the  absorbing  states,  such  as  state  4,  (*2,0*),  appear  only  in 
the  destination  column.  The  DEATHIF  definitions  and  ONEDEATH  command 
are  not  printed.  Listing  3  is  the  STEM  output  file  for  the  model  above,  run  for  10 
time  units. 


Listing  3.  2  vs.  2  Engagement  STEM  Output 


Model  =  C:\Markovl\fighters\Report2x2.mod 
-  RUN  #1 

D-STATE  PROBABILITY  ACCURACY 


4  8. 6580085643 0E-0001 

6  1. 5151310349 5E-0002 

7  1.08204982318E-0001 

8  1.08204982318E-0002 


TOTAL  9.99977647330E-0001  11  DIGITS 
STATE  PROBABILITY 


1  2 . 78946809287E-0010 

2  9. 5861846700 9E-0009 

3  1.22878668130E-0006 

4  8. 6580085643 0E-0001 

5  2 . 11140182554E-0005 

6  1. 5151310349 5E-0002 

7  1.08204982318E-0001 

8  1.08204982318E-0002 


Figure  B-3  shows  the  STEM  results  for  increasing  units  of  combat  time.  As  dis¬ 
cussed  in  Chapter  2,  we  use  the  reciprocal  of  the  lowest  rate  as  a  rule  of  thumb  for 
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setting  the  model  run  time.  In  this  case,  the  lowest  rate  is  KRATERED  =  0.01 
and  its  reciprocal  is  100.  By  Time  =  70,  the  total  probability  of  being  in  an  ab¬ 
sorbing  state  is  0.9995,  and  at  Time  =  100,  it  is  0.99998,  indicating  that  our  rule 
of  thumb  is  good.  The  state  probabilities  at  Time  =100  for  the  individual  absorb¬ 
ing  states  are  (2B,  OR)  =  0.8658,  (IB,  OR)  =  0.1082,  (0B,  2R)  =  0.0152,  and  (0B, 
1R)  =  0.0108. 

Figure  B-3.  2  vs.  2  Engagement  Results  versus  Engagement  Time 


Elapsed  Time  Units 


Example  2 — 2  vs.  2  Engagement  with  Blue  Missiles 

This  example  adds  the  tracking  of  Blue  missiles  to  the  simple  2  vs.  2  engagement. 
We  assume  that  Blue  aircraft  each  carry  6  missiles,  each  having  a  known  single 
shot  probability  of  kill,  Pk.  It  is  important  to  understand  in  this  discussion  that 
Blue  and  Red  kills  are  still  dependent  only  on  the  kill  rate  ratio,  and  that  missile 
Pk  is  used  only  to  calculate  missile  consumption.  For  simplicity  and  clarity,  we 
assume,  for  this  example,  that  Red  aircraft  have  unlimited  missiles. 

Our  basic  approach  is  to  expand  the  state  space  to  (bO,  bl,  b2,  b3,  b4,  b5,  b6,  r)  to 
include  separate  states  for  Blue  aircraft  having  0  to  6  missiles.  We  still  need  only 
one  state  vector  for  Red  aircraft.  At  the  start  of  an  engagement,  we  have  2  Blue 
aircraft  in  the  b6  6-missile  state  and  2  Red  aircraft  in  the  r  Red  state. 


We  still  transition  from  Red  states  in  only  two  ways:  Red  can  kill  a  Blue,  or  vice 
versa.  Now,  however,  we  can  transition  from  a  Blue  state  in  three  ways:  Blue  can 
kill  a  Red,  Blue  can  fire  and  miss,  or  Red  can  kill  a  Blue.  To  formulate  the  trans¬ 
fer  statements,  we  need  to  determine  the  missile  miss  rate.  We  know  the  Blue  kill 
rate,  X,  and  we  know  the  missile  single  shot  probability  of  kill,  Pk-  We  want  to 
estimate  the  missile  usage,  including  misses.  We  assume,  for  now,  that  the  suc¬ 
cessful  missile  usage  is  one  missile  per  kill,  i.e.,  we  do  not  fire  salvos.  We  also 
assume  the  kill  rate  includes  both  the  missile  that  hit  and  any  that  missed. 
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Now,  we  define  /u  to  be  the  missile  firing  rate,  and  let  n  be  the  number  of  Blue 
aircraft  engaged.  The  successful  missile  rate  is  then:  n*Pk  *  |u,  and  the  unsuccess¬ 
ful  missile  rate  (miss  rate)  is  n*(l-Pk)*p. 

Based  on  our  assumptions,  the  successful  missile  rate  must  equal  the  kill  rate,  i.e., 


n*Pk*p  =  n*k 

and,  therefore,  the  missile  firing  rate  is 


n  *  p  =  n  *  — ,  and  p  =  —  . 


[Eq.  B-l] 


[Eq.  B-2] 


Substituting  for  p,  the  miss  rate,  n*(l-Pk)*p,  becomes 

n  *  (1  -  Pi. )  *  —  or  n  *  ^  ~  ^  *  A, .  [Eq.  B-31 

Pk  Pk 

Figure  B-4  is  the  state  transition  diagram  for  the  single  missile  per  kill  case. 

Figure  B-4.  State  Diagram  for  a  Single  Missile  per  Kill  Engagement 


Missile  Firing  Rate  =  p 

Missile  Kill  Rate  =  PkP 

Missile  Miss  Rate  =(1-Pk)p 

Engagement  Kill  Rates  =  A,b ,  A,r 
Aircraft  in  Engagement  =  n  Blue,  m  Red 


We  can  generalize  the  derivation  above  in  a  straightforward  manner  for  cases  in 
which  multiple  missiles  are  fired  per  engagement  by  noting  that  the  failure  prob¬ 
ability  for  x  failures  is  (1  -  Pk)x  and  substituting  [1  -  (1  -Pkf]  for  Rp 

♦  Define  //  to  be  the  missile  firing  rate. 

♦  Let  n  be  the  number  of  Blue  aircraft  engaged. 
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♦  Let  x  be  the  number  of  missiles  fired  per  kill. 

The  successful  missile  rate  is  n*[  1-  (1  -  Pk)A]  *  p,  and  the  unsuccessful  missile 
rate  (miss  rate)  is  n*(l  -  Pk)x*p. 

Based  on  our  assumptions,  the  successful  missile  rate  must  equal  the  kill  rate,  i.e., 

n*[l  -  (1  -  Pk)x]*p  =  n*L  [Eq.B-4] 

and,  therefore,  the  missile  firing  rate  is 


n  *  p  =  n  * 


X 


, and  p  = 


X 


i-(i-Pk)’ 


p-(i-pk)" 

The  miss  rate,  n*(l  -  Pk)x*p,  now  becomes 

-  (1-Pk)? 


n*(l-Pk) 


x  *  . 


or  n 


i-(i-Pi)" 


X. 


For  x  =  2  missiles  per  kill,  we  have  a  miss  rate  of 


n 


*  1  ~  2Pk  +  Pk  *  y 

2Pk-Pk2 


[Eq.  B-5] 


[Eq.  B-6] 


[Eq.  B-7] 


In  addition  to  missile  counting,  we  also  add  the  capability  for  either  Red  or  Blue 
to  exit  the  combat  when  losses  reach  a  preset  level  by  assigning  minimum  values 
(minb,  minr)  to  the  DEATHIF  conditions.  If  minb  =  0  and  minr  =  0,  the 
model  corresponds  to  fighting  to  annihilation;  higher  values  will  generate  models 
reflecting  disengagement  breakpoints.  Listing  4  is  the  ASSIST  code  for  the  2  vs.  2 
case  with  missile  tracking. 

Listing  4.  ASSIST  Code  for  a  2  vs.  2  Engagement  with  Blue  Missile  Tracking 

(*  ASSIST  code  for  2  Blue  vs.  2  Red  with  6  missiles  per  Blue  aircraft  *) 

LIST  =3;  (*  1  lists  accumulated  death  state  probabilities  only  *) 

(*  2  lists  all  state  probabilities  *) 

(*  3  lists  transitions  and  all  state  probabilities  *) 

ONEDEATH  OFF;  (*  OFF  enumerates  all  absorbing  (Death)  states  *) 

(*  comment  out  to  consolidate  absorbing  states  *) 

Time  =  100; 

(*  equipment  list  *) 
nblue  =2;  (*  blue  fighters  *) 

nred  =2;  (*  red  fighters  *) 

(*  kill  rates  *) 

k_rate_blue  =0.1;  (*  blue  kill  rate  *) 

k_rate_red  =  0.01;  (*  red  kill  rate  *) 

pk  =  0.85;  (*  missile  kill  probability  *) 

miss_rate  =  ( 1-pk) /pk*k_rate_blue;  (*  missile  miss  rate  *) 
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(*  minimum  number  of  combatants  for  state  pruning  *) 
min_b  =0;  (*  minimum  number  of  blue  aircraft  *) 
min_r  =0;  (*  minimum  number  of  red  aircraft  *) 

SPACE  = 

(nO : 0 . nblue, nl : 0 . nblue, n2 : 0 . nblue , n3 : 0 . nblue, n4 : 0 . nblue, n5 : 0 . nblue , n6 : 0 . nblue, 
red: 0 . nred) ; 

START  =  (0, 0, 0, 0, 0, 0, nblue, nred)  ; 

DEATH IF  (n0+nl+n2+n3+n4+n5+n6) <=min_b; ; 

DEATH IF  red  <=  min_r; 

(*  transition  cases:  Blue  kill  then  Red  kill*) 

IF  (n6>0)  AND  (red>0)  THEN 

TRANTO  n6=n6-l,  n5=n5+l,  red=red-l  BY  n6*k_rate_blue; 

TRANTO  n6=n6-l,  n5=n5+l  BY  n6*miss_rate; 

TRANTO  n6=n6-l  BY  (n6/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

IF  (n5>0)  AND  (red>0)  THEN 

TRANTO  n5=n5-l,  n4=n4+l,  red=red-l  BY  n5*k_rate_blue; 

TRANTO  n5=n5-l,  n4=n4+l  BY  n5*miss_rate ; 

TRANTO  n5=n5-l  BY  (n5/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

IF  (n4>0)  AND  (red>0)  THEN 

TRANTO  n4=n4-l,  n3=n3+l,  red=red-l  BY  n4*k_rate_blue; 

TRANTO  n4=n4-l ,  n3=n3+l  BY  n4*miss_rate; 

TRANTO  n4=n4-l  BY  (n4/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

IF  (n3>0)  AND  (red>0)  THEN 

TRANTO  n3=n3-l ,  n2=n2+l,  red=red-l  BY  n3*k_rate_blue; 

TRANTO  n3=n3-l,  n2=n2+l  BY  n3*miss_rate; 

TRANTO  n3=n3-l  BY  (n3/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

IF  (n2>0)  AND  (red>0)  THEN 

TRANTO  n2=n2-l ,  nl=nl+l,  red=red-l  BY  n2*k_rateJolue; 

TRANTO  n2=n2-l ,  nl=nl+l  BY  n2 *miss_rate ; 

TRANTO  n2=n2-l  BY  (n2/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

IF  (nl>0)  AND  (red>0)  THEN 

TRANTO  nl=nl-l,  n0=n0+l,  red=red-l  BY  nl*k_rate_blue; 

TRANTO  nl=nl-l ,  n0=n0+l  BY  nl*miss_rate; 

TRANTO  nl=nl-l  BY  (nl/ (n0+nl+n2+n3+n4+n5+n6) ) *red*k_rate_red; 

ENDIF; 

(*  This  line  is  commented  out  based  on  Blues  ability  to  disengage  after  firing 
all  his  missiles  *) 

( *  I F  (red>0)  AND  (n0>0)  TRANTO  n0=n0-l  BY 
(n0/ (n0+nl+n2+n3+n4+n5+n6 ) ) *red*k_rate_red; * ) 


Adding  the  missile  Pk,  missile  miss  rate,  Blue  and  Red  aircraft  minimums,  the 
expanded  state  space,  and  the  starting  conditions  is  straightforward.  The  expan¬ 
sion  of  the  transfer  statements  deserves  some  discussion. 

The  basic  transition  logic  is  as  follows: 

♦  Separate  sets  of  transfers  are  established  for  each  Blue  missile  state. 

♦  State  transfers  occur  only  when  Blue  aircraft  are  in  the  state  and  Red  air¬ 
craft  are  available. 
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♦  Blue  kills  at  a  rate  proportional  to  the  total  number  of  Blue  aircraft.  When 
Blue  makes  a  kill,  one  aircraft  is  removed  from  the  Blue  state,  one  aircraft 
is  added  to  the  Blue  state  having  one  fewer  missiles,  and  one  aircraft  is 
removed  from  the  Red  state. 

♦  Blue  misses  at  a  rate  proportional  to  the  number  of  aircraft  in  the  Blue 
state.  When  Blue  misses,  one  aircraft  is  removed  from  the  Blue  state,  one 
aircraft  is  added  to  the  Blue  state  having  one  fewer  missiles,  but  no  aircraft 
are  removed  from  the  Red  state. 

♦  Red  kills  Blue  at  a  rate  proportional  to  the  number  of  Red  aircraft  and  the 
fraction  of  Blue  aircraft  in  the  Blue  state.  When  Red  makes  a  kill,  one  air¬ 
craft  is  removed  from  the  Blue  state. 

♦  Only  Red  kills  are  possible  when  Blue  has  no  missiles  left.  This  TRANTO 
statement  is  commented  out  in  the  example  listing  using  “(*”  and  “*)”  op¬ 
erators  based  on  the  assumption  that  Blue  can  disengage  at  will.  Another 
alternative  would  be  to  give  Blues  with  no  missiles  a  “guns  only”  kill  rate. 


Running  the  ASSIST  program  above  generates  a  model  having  103  states,  includ¬ 
ing  34  absorbing  (Death)  states,  and  301  transitions.  Listing  5  shows  the  first  and 
last  few  lines  of  the  model  file.  The  state  notation  in  the  model  shows  the  state 
number  and  the  contents  of  each  element  in  the  state.  This  model  has  eight  ele¬ 
ments  in  each  state.  The  first  seven  elements  are  Blue  aircraft  with  0  to  6  missiles, 
and  the  last  element  is  Red  aircraft.  Thus,  State  1  is  1(0, 0,0, 0,0, 0,2, 2),  indicating  2 
Blues  with  6  missiles  and  2  Reds.  The  first  transition  is  to  2(0, 0,0, 0,0, 1,1,1),  indi¬ 
cating  that  one  Blue  has  fired  a  missile  and  killed  one  Red.  A  Blue  absorbing  state 
corresponds  to  0s  in  all  of  the  first  seven  elements,  and  a  Red  absorbing  state  has 
0  in  the  eighth  element. 

Listing  5.  Model  File  Segments  for  2  vs.  2  Engagement 
with  Blue  Missile  Tracking 


LIST  =  3; 

NBLUE  =  2; 

NRED  =  2; 

K_RATE_BLUE  =  0.1; 

K_RATE_RED  =  0.01; 

PK  =  0.85; 

MISS_RATE  =  (1-PK) / PK*K_RATE_BLUE ; 

MIN_B  =  0; 

MIN  R  =  0; 


1  (* 

CM 

CM 

O 

O 

O 

O 

O 

O 

*) 

2  (* 

o 

o 

o 

o 

o 

*) 

=  2*K  RATE  BLUE; 

1  (* 

0,0, 0,0, 0,0, 2, 2 

*) 

3  ( * 

0,0, 0,0, 0,1, 1,2 

*) 

=  2*MISS  RATE; 

1  (* 
*2  *' 

0,0, 0,0, 0,0, 2, 2 
K  RATE  RED; 

*) 

4  (* 

CM 

O 

O 

O 

O 

O 

O 

*) 

=  (2/ (0  +  0  +  0  +  0  +  0  +  0+2) 

2  (* 

0,0, 0,0, 0,1, 1,1 

*) 

5  (* 

0,0, 0,0, 0,2, 0,0 

*) 

=  1*K  RATE  BLUE; 

2  (* 

5 - 1 

o 

o 

o 

o 

o 

*) 

6  (* 

o 

CM 

O 

O 

O 

O 

O 

*) 

=  1*MISS  RATE; 

2  (* 
*  \  *' 

0,0, 0,0, 0,1, 1,1 
K  RATE  RED; 

*) 

7  (* 

O 

O 

o 

o 

o 

o 

*) 

=  (1/ (0  +  0  +  0  +  0  +  0  +  1  +  1) 

2  (* 

0,0, 0,0, 0,1, 1,1 

*) 

8  (* 

0,0, 0,0, 1,0, 1,0 

*) 

=  1*K  RATE  BLUE; 

2  (* 

o 

o 

o 

o 

o 

*) 

9  (* 

o 

o 

o 

o 

o 

*) 

=  1*MISS  RATE; 

2  (* 

0,0, 0,0, 0,1, 1,1 

*) 

10  r 

*  0,0, 0,0, 0,0, 1,1 

*■  ’ 

)  =  (1/ (0+0+0+0+0+1+1 
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*1*K  RATE  RED; 


3  ( * 

0,0, 0,0, 0,1, 1,2 

*)  , 

CT> 

O 

O 

O 

O 

O 

to 

O 

I—4 

=  1*K  RATE  BLUE; 

3  ( * 

0,0, 0,0, 0,1, 1,2 

*)  , 

.  11 (*  0,0, 0,0, 0,2, 0,2  *) 

=  1*MISS  RATE; 

3  ( * 
*2*: 

0,0, 0,0, 0,1, 1,2 
K  RATE  RED; 

*)  , 

.  12 (*  0,0, 0,0, 0,1, 0,2  *) 

=  (1/ (0+0+0+0+0+1+1) 

3  (* 

0,0, 0,0, 0,1, 1,2 

*)  , 

.  9 ( *  0,0, 0,0, 1,0, 1,1  *) 

=  1*K  RATE  BLUE; 

3  (* 

0,0, 0,0, 0,1, 1,2 

*)  , 

.  13  (*  0,0, 0,0, 1,0, 1,2  *) 

=  1*MISS  RATE; 

3  ( * 

0,0, 0,0, 0,1, 1,2 

*)  , 

.  4 (*  0,0, 0,0, 0,0, 1,2  *) 

=  (1/ (0+0+0+0+0+1+1) ) 

97 (*  1,0, 1,0, 0,0, 0,2  *),  99 ( *  1 , 1 , 0 , 0 , 0 , 0 , 0 , 1  *)  =  1*K_RATE_BLUE; 

97 (*  1,0, 1,0, 0,0, 0,2  *),  100 (*  1,1, 0,0, 0,0, 0,2  *)  =  1*MISS_RATE; 

97  (*  1,0, 1,0, 0,0, 0,2  *),  7  6  ( *  1 , 0 , 0 , 0 , 0 , 0 , 0 , 2  *)  =  ( 1 / ( 1+0+1+0+0+0+0 ) ) 
*2*K_RATE_RED; 

99 ( *  1, 1, 0, 0, 0, 0, 0, 1  *),  1 0 1  ( *  2, 0, 0, 0, 0, 0, 0, 0  *)  =  1*K_RATE_BLUE; 

99 ( *  1, 1, 0, 0, 0, 0, 0, 1  *),  102 (*  2, 0, 0, 0, 0, 0, 0, 1  *)  =  1*MISS_RATE; 

99 ( *  1,1, 0,0, 0,0, 0,1  *),  72  ( *  1,0, 0,0, 0,0, 0,1  *)  =  ( 1 / ( 1+1+0+0+0+0+0 ) ) 
*1*K_RATE_RED; 

100(*  1, 1, 0, 0, 0, 0, 0, 2  *),  102 (*  2, 0, 0, 0, 0, 0, 0, 1  *)  =  1*K_RATE_BLUE; 

100 ( *  1,1, 0,0, 0,0, 0,2  *),  103 (*  2, 0,0, 0,0, 0,0, 2  *)  =  1*MISS_RATE; 

100  ( *  1,  1, 0,0, 0,0, 0,2  *),  7  6  ( *  1,0, 0,0, 0,0, 0,2  *)  =  ( 1 / (1+1+0+0+0+0+0 ) ) 
*2*K_RATE_RED; 

(*  NUMBER  OF  STATES  IN  MODEL  =  103  *) 

(*  NUMBER  OF  TRANSITIONS  IN  MODEL  =  285  *) 


At  100  time  units,  the  engagement  is  0.99994  complete;  the  sums  of  the  absorbing 
missile  tracking  states  agree  within  four  significant  figures  with  the  sums  of  the 
absorbing  states  for  the  simple  2  vs.  2  example.  Figure  B-5  combines  the  prob¬ 
abilities  for  the  absorbing  states  to  show  the  marginal  probabilities  of  Red  victory 
and  the  marginal  probabilities  of  Blue  victory  with  specific  numbers  of  aircraft 
and  missiles  after  a  2  vs.  2  engagement  of  100  time  units. 

Figure  B-5.  2  vs.  2  Engagement  Aircraft  and  Missile  Configuration  Probabilities 

Marginal  Probabilities  of  Red  Win  and  of  Blue  Wins  for  N  Blue  with  X  Missiles  Remaining 
(Time=100,  Kb=0.1,  Kr=0.01,  Missile  Pk=0.85) 
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Example  3 — 4  vs.  4+4+4  Sequential  Engagement 

The  standard  SLAACM  campaign  scenario  has  4  Blue  defenders  engaging  4  Red 
long-range  escorts  (LEs),  4  close  escorts  (CEs),  and  4  bombers  sequentially.  This 
sequential  engagement  can  be  modeled  using  ASSIST  by  imposing  on  the 
TRANTO  statements  the  conditions  that  Blue  cannot  engage  the  close  escorts  or 
bombers  until  all  the  LEs  are  dead  and  cannot  engage  the  bombers  until  all  the 
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CEs  are  also  dead.  We  include  individual  kill  rates  for  the  three  pairs  of  combat¬ 
ants.  Listing  6  is  the  ASSIST  code  for  the  sequential  engagement.  The  resulting 
model  has  38  states  and  48  transitions  when  the  Blue  breakpoint  is  2  losses 
(minb  =  2  in  the  listing),  and  has  64  states  and  96  transitions  when  Blue  fights  to 
annihilation  (min  b  =  0).  Figure  B-6  shows  the  marginal  probabilities  for  the  out¬ 
comes  for  the  2  Blue  breakpoint,  and  Figure  B-7  shows  the  corresponding  ex¬ 
pected  number  of  surviving  aircraft  types. 

Listing  6.  4  Blue  Sequentially  Engaging  4  Red  LE,  4  Red  CE, 
and  4  Red  Bombers 


(*  ASSIST  model  for  Fighter  Combat  *) 

(*  This  is  multistage  combat  where  blue  aircraft  sequentially  engage*) 
(*  three  classes  of  red  aircraft  *) 

(*  This  version  does  not  include  missiles  *) 

LIST  =  3;  (*  1  lists  accumulated  death  state  probabilities  only  *) 

(*  2  lists  all  state  probabilities  *) 

(*  3  lists  transitions  and  all  state  probabilities  *) 

ONEDEATH  OFF;  (*  OFF  enumerates  all  absorbing  (Death)  states  *) 

(*  comment  out  to  consolidate  absorbing  states  *) 

(*"Time  =  5  to  100  By  5;"*) 

(*  the  time  is  set  here  -  no  matter  what  the  STEM  GUI  shows  *) 

Time  =  1000; 

(*  equipment  list  *) 


nblue  = 

4; 

(* 

blue 

:  fighters  *) 

nredl  = 

4; 

(* 

red 

fighters  -  lead  escorts  *) 

nred2  = 

4; 

(* 

red 

fighters  -  close  escorts  *) 

nred3  = 

4; 

(* 

red 

bombers  or  fighter  bombers  * 

min  b  = 

2; 

min  rl 

=  0; 

min  r2 

=  0; 

min  r3 

=  0; 

(*  kill  probabilities  *) 

(*  Note:  the  1st  12  characters  of  a  variable  must  be  unique  *) 
blue_redl_k_rate  =  0.1;  (*  0.1  *) 

blue_red2_k_rate  =  0.1;  (*  0.2  *) 

blue_red3_k_rate  =  0.1;  (*  0.4  *) 

k_rate_redl  =  0.02;  (*  0.05  *) 

k_rate_red2  =  0.01;  (*  0.01  *) 

k_rate_red3  =  0.005;  (*  0.005  *) 

SPACE  =  (blue:  O.nblue,  redl :  O.nredl,  red2 :  0.nred2,  red3 :  0.nred3); 

START  =  (nblue,  nredl,  nred2,  nred3) ; 

DEATH IF  (blue  =  min_b)  AND  ((redl>0)  OR  (red2>0)  OR  (red3>0) ) ; 

DEATH IF  (redl  =  min_rl)  AND  (red2  =  min_r2)  AND  (red3  =  min_r3)  AND  (blue>0); 

(*  Phase  1  combat  *) 

(*  Blue  kills  *) 

IF  (blue  >  min_b)  and  (redl  >  min_rl)  TRANTO  redl  =  redl-1  BY 
blue *blue_redl_k_r ate; 

(*  Red  kills  *) 

IF  (redl  >  min_rl  )  and  (blue>min_b)  TRANTO  blue  =  blue-1  BY  redl *k_rate_redl ; 

(*  Phase  2  combat  *) 

(*  Blue  kills  *) 

IF  (blue  >  min_b)  and  ( redl=min_r 1 )  and  (red2  >  min_r2)  TRANTO  red2  =  red2-l 
BY  blue*blue_red2_k_rate; 

(*  Red  kills  *) 
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IF  (red2  >  min_r2  )  and  (redl=min_rl )  and  (blue>min_b)  TRANTO  blue  =  blue-1  BY 
red2  *  k_rate_red2 ; 

(*  Phase  3  combat  *) 

(*  Blue  kills  *) 

IF  (blue  >  min_b)  and  (red3  >  min_r3)  TRANTO  red3  =  red3-l  BY 
blue*blue_red3_k_rate; 

(*  Red  kills  *) 

IF  (red3  >  min_r3  )  and  (redl=min_rl )  and  (red2=min_r2 )  and  (blue>min_b) 

TRANTO  blue  =  blue-1  BY  red3*k_rate_red3 ; 

Figure  B-6.  Sequential  Engagement  Marginal  Probabilities 

Aircraft  Count  per  Type  Probabilities 
4B  x  4R1 ,  then  x  4R2  then  x  4R3 
(B:R  Kill  Rates  -  5:1,  10:1,  20:1,  Blue  Min  =  2) 
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Figure  B-7.  Sequential  Engagement  Surviving  Aircraft 


Post  Engagement  Aircraft 
4B  x  4R1 ,  then  x  4R2  then  x  4R3 
(B:R  Kill  Rates  -  5:1,  10:1,  20:1,  Blue  Min  =  2) 
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Example  4 — 4  vs.  4+4+4  Sequential  Engagement  with  Missiles 

The  last  example  is  the  standard  sequential  engagement  with  missiles  on  the  Blue 
aircraft.  Listing  7  is  the  ASSIST  code  for  this  case.  This  listing  generates  a  model 
having  3,695  states  and  21,729  transitions  for  the  2  Blue  breakpoint  case  and 
3,823  states  and  23,454  transitions  for  the  Blue  fight  to  annihilation  case.  Figure 
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B-8  shows  the  marginal  probabilities  of  missile  counts  combined  with  infonnation 
on  how  many  aircraft  are  carrying  the  missiles  for  the  case  in  which  Blue  breaks 
after  two  losses.  Figure  B-9  shows  the  remaining  aircraft  count  by  type.  Table  B-l 
contains  the  summary  results  for  both  the  2  Blue  breakpoint  and  Blue  fight  to  an¬ 
nihilation  cases. 

Listing  7.  4  Blue  vs.  4  Red  LE,  4  Red  CE,  and  4  Bombers  Sequentially 

with  Blue  Missiles 


(*  ASSIST  model  for  *) 

(*  This  version  includes  sequential  combat  with  3  types  of  red  fighters  *) 
(*  and  includes  Blue  missiles  *) 

(*  Fighter  Combat  using  Dave  Lee's  construct  for  blue  missile  counting*) 

(*  7/20/05  This  version  includes  missiles  for  blues  *) 

(*  7/20/05  THIS  VERSION  INCLUDES  6  BLUE  MISSILES 

(*  THIS  VERSION  CAN  LIMIT  STATES  BASED  ON  A  MINIMUM  NUMBER  OF  COMBATANTS  *) 
(*  The  kill  rate  is  for  one-on-one  and  is  multiplied  for  N  killers  *) 

(*  Assume  one  missile  per  kill  even  with  N  killers  *) 

(*  Therefore,  the  kill  rate  includes  one  hit  and  all  misses  *) 

(*  Find  miss  rate  based  on  ratio  of  Pmiss:Pk  *) 

LIST  =3;  (*  1  lists  accumulated  death  state  probabilities  only  *) 

(*  2  lists  all  state  probabilities  *) 

(*  3  lists  transitions  and  all  state  probabilities  *) 

ONEDEATH  OFF;  (*  OFF  enumerates  all  absorbing  (Death)  states  *) 

(*  comment  out  to  consolidate  absorbing  states  *) 

(*"Time=10  TO+  100  BY  10;"*) 

Time  =  1000; 

(*  equipment  list  *) 
nblue  =4;  (*  blue  fighters  *) 

nredl  =4;  (*  red  fighters  *) 

nred2  =  4 ; 
nred3  =  4; 


(*  missile  equipage  *) 
b_msls  =6;  (*  blue  missiles  *) 

(*  the  first  12  characters  in  a  variable  name  are  significant  and  must  be 
unique  *) 


(*  kill  rates  *) 

redl_k_rate_blue  =  0.1;  (*  B:R1  =  5:1  *) 

red2_k_rate_blue  =  0.1;  (*  B:R2  =  10:1  *) 

red3_k_rate_blue  =  0.1;  (*  B:R3  =  20:1  *) 

k_rate_redl  =  0.02; 
k_rate_red2  =  0.01; 
k  rate  red3  =  0.005; 


(*  blue  missiles  miss  rates  *) 
bpk  =  0.85;  (*  missile  kill  probability  *) 

redl_bmiss_rate  =  ( 1-bpk) /bpk*redl_k_rate_blue;  (*  missile  miss  rate  against 
redl  *) 

red2_bmiss_rate  =  ( 1-bpk) /bpk*red2_k_rate_blue;  (*  missile  miss  rate  against 
red2  *) 

red3_bmiss_rate  =  ( 1-bpk) /bpk*red3_k_rate_blue;  (*  missile  miss  rate  against 
red3  *) 

(*  minimum  number  of  combatants  for  state  pruning  *) 
min_b  =2;  (*  minimum  number  of  blue  aircraft  *) 

min_rl  =0;  (*  minimum  number  of  redl  aircraft  *) 

min_r2  =  0; 
min  r3  =  0; 
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SPACE  =  (b : array [0 .b_msls] of  O.nblue,  redl :  O.nredl,  red2 :  0.nred2,  red3 : 

0 . nred3) ; 

(*  states  0  to  max  missiles  *) 

(*  THIS  STATEMENT  MUST  BE  CHANGED  WHEN  THE  MISSILE  COUNT  IS  CHANGED!  *) 

START  =  (0, 0, 0, 0, 0, 0, nblue, nredl, nred2, nred3)  ; 

DEATH IF  (sum (b) <=  min_b)  OR  (b [ 0 ] =4 ) ; 

DEATH IF  (redl  =  min_rl)  AND  (red2  =  min_r2)  AND  (red3  =  min_r3) ; 

(*  state  transitions  and  kills  *) 

(*  advanced  escorts  *) 

FOR  I  IN  [l.b_msls] 

IF  (b [I] >0)  AND  (sum (b) >min_b)  AND  ( redl>min_rl )  THEN 

TRANTO  b [ I ] =b [ I ] -1 ,  b[I-l]=b[I-l]+l,  redl  =redl-l  BY  b [ I ] *redl_k_rate_blue; 
TRANTO  b [ I ] =b [ I ] -1 ,  b [ I— 1 ] =b [ I— 1 ] +1  BY  b [ I ] *redl_bmiss_rate; 

TRANTO  b [ I ] =b [ I ] — 1  BY  (b [ I ] /sum (b) ) *redl *k_rate_redl ; 

ENDIF; 

ENDFOR; 

(*  no  missile  case  against  advanced  escorts  *) 

IF  (b [0] >0)  AND  (sum (b) >min_b)  AND  (redl  >  min_rl)  THEN 

(*  TRANTO  redl=redl-l  BY  b [ 0 ] *redl_k_rate_blue;  *)  (*  NO  BLUE  KILLS  WITHOUT 

MISSILES  *) 

(*  TRANTO  b [ 0 ] =b [ 0 ] — 1  BY  (b [ 0 ] /sum (b) ) *redl *k_rate_redl ;  *)  (*  BLUE  ESCAPES  IF 

NO  MISSILES  *) 

ENDIF; 

(*  close  escorts  *) 

IF  redl  =  min_rl  THEN 
FOR  I  IN  [l.b_msls] 

IF  (b [I] >0)  AND  (sum (b) >min_b)  AND  (red2  >  min_r2)  THEN 

TRANTO  b [ I ] =b [ I ] -1 ,  b [1-1] =b [ I  — 1 ] +1 ,  red2  =red2-l  BY  b [ I ] * red2_k_rate_blue ; 
TRANTO  b [ I ] =b [ I ] -1 ,  b [ I— 1 ] =b [1-1] +1  BY  b [ I ] *red2_bmiss_rate; 

TRANTO  b [ I ] =b [ I ] — 1  BY  (b [ I ] /sum (b) ) *red2 *k_rate_red2 ; 

ENDIF; 

ENDFOR; 

(*  no  missile  case  against  close  escorts  *) 

IF  (b [0] >0)  AND  (sum (b) >min_b)  AND  (red2  >  min_r2)  THEN 

(*  TRANTO  red2=red2-l  BY  b [ 0 ] *red2_k_rate_blue;  *)  (*  NO  BLUE  KILLS  WITHOUT 

MISSILES  *) 

(*  TRANTO  b [ 0 ] =b [ 0 ] — 1  BY  (b [ 0 ] /sum (b) ) *red2 *k_rate_red2 ;  *)  (*  BLUE  ESCAPES  IF 

NO  MISSILES  *) 

ENDIF; 

ENDIF; 

(*  bombers  *) 

IF  (red2  =  min_r2)  AND  (redl  =  min_rl)  THEN 
FOR  I  IN  [l.b_msls] 

IF  (b [I] >0)  AND  (sum (b) >min_b)  AND  (red3  >  min_r3)  THEN 

TRANTO  b [ I ] =b [ I ] -1 ,  b [1-1] =b [1-1] +1,  red3  =red3-l  BY  b [ I ] *red3_k_rate_blue; 
TRANTO  b [ I ] =b [ I ] -1 ,  b [ I— 1 ] =b [1-1] +1  BY  b [ I ] *red3_bmiss_rate; 

TRANTO  b [ I ] =b [ I ] -1  BY  (b [ I ] / sum (b) ) *red3*k_rate_red3 ; 

ENDIF; 

ENDFOR; 

(*  no  missiles  case  against  bombers  *) 

IF  (b [0] >0)  AND  (sum (b) >min_b)  AND  (red3  >  min_r3)THEN 

(*  TRANTO  red3=red3-l  BY  b [ 0 ] *red3_k_rate_blue;  *)  (*  NO  BLUE  KILLS  WITHOUT 

MISSILES  *) 

(*  TRANTO  b [ 0 ] =b [ 0 ] — 1  BY  (b [ 0 ] /sum (b) ) *red3*k_rate_red3 ;  *)  (*  BLUES  ESCAPES 

IF  NO  MISSILES  *) 

ENDIF; 

ENDIF; 

(*  FOR  NO  MISSILE  CASES:  *) 

(*  use  only  red  rates  if  blue  has  no  missiles  left  and  cannot  escape  *) 

(* —  comment  out  all  statements  if  blue  can  escape  *) 

(* —  include  blue  kill  statements  to  continue  combat  without  missiles  *) 

(*  possibly  add  rates  above  and  TRANTO  statments  here  for  Blue  guns  *) 
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Marginal  Probability 


Figure  B-8.  Sequential  Engagement  Missile  Marginal  Probabilities 
with  Carrier  Aircraft  Information 


4v12  Sequential  Engagement  Aircraft  and  Missiles  Remaining 
(BvLE-5:1,  BvCE-10:1,  BvBm-20:1,  Min  Blue  =  2,  T=1000  units) 
(Blue  with  0  missiles  escape) 
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Figure  B-9.  Sequential  Engagement  with  Missiles,  Remaining  Aircraft 


4v12  Sequential  Engagement  Aircraft  Remaining 
(BvLE-5:1,  BvCE-10:1,  BvBm-20:1 ,  Min  Blue  =  2,  T=1000  units) 
(Blue  with  0  missiles  escape) 
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Analyzing  Large  State-Space  Engagement  Models  with  ASSIST  and  STEM 


Table  B-1.  Summary  Engagement  Results 


Scenario 

Blue 

remaining 

Red 

remaining 

Missiles 

used 

2  Blue  breakpoint 

3.2 

2.1 

16 

Blue  fight  to  annihilation 

3.0 

0.8 

18 

Summary 


In  this  appendix,  we  have  shown  how  the  NASA-developed  Markov  tools 
ASSIST  and  STEM  can  be  used  for  analysis  of  increasingly  complex  air  combat 
engagements.  We  have  shown  through  examples  how  ASSIST  is  particularly  use¬ 
ful  for  supporting  parametric  analysis  involving  large  state-space  problems.  The 
NASA  tools  have  proved  useful  for  standalone  analyses,  for  prototyping  of 
SLAACM  engagements,  and  for  independent  confirmation  of  SLAAM  calcula¬ 
tions. 
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Appendix  C 

Alternatives  to  Computational 
Probabilistic  Tools 


In  this  appendix,  we  compare  results  of  two  familiar  models,  the  Lanchester 
force-on- force  model1  and  simulation  modeling,  with  the  results  of  the  probabilis¬ 
tic  engagement  model  applied  to  some  specific  problems. 

Lanchester  Force-on-Force  Model 


One  alternative  to  the  probabilistic  engagement  model  is  the  classic  Lanchester 
force-on-force  model.  In  this  engagement  model,  the  time  rate  of  change  of  the 
number  of  members  of  a  side  is  proportional  to  the  number  of  members  of  the 
other  side.  This  leads  to  the  initial  value  problem  Equation  C-l. 

nb  = -krnr 

nr=-kbnb  .  [Eq.  C-l] 

nb(0)  =  b0;  nr(0)  =  r0 

The  differential  equations  in  Equation  C- 1  can  be  solved  readily  in  closed  form 
with  exponential  Emotions.2 

The  Lanchester  force-on- force  model  is,  of  course,  entirely  deterministic.  In  fact, 
its  outcomes  are  determined  by  values  of  a  certain  ratio,  called  the  strength  ratio. 
To  see  this,  we  may  multiply  the  first  differential  equation  by  kt,nb  and  the  second 
by  krnr,  subtract  the  second  product  from  the  first,  and  integrate.  This  leads  to 

kb(b2-nj;)  =  kt(r02-n;).  [Eq.  C-2] 

Using  Equation  C-2,  we  can  develop  the  strength  ratio,  S,  which  is  defined  by 
Equation  C-3. 


s  =  Mbii-b2k) 

kr(r0  -rbk) 


[Eq.  C-3] 


Examining  Equation  C-3,  one  sees  that  if  Red  breaks  when  they  have  rbk  mem¬ 
bers,  and  Blue  breaks  when  they  have  bbk,  then  Blue  will  win  in  finite  time  when 

1  James  G.  Taylor,  Lanchester  Models  of  Warfare,  Volume  I  (Operations  Research  Society  of 
America,  March  1983),  p  83. 

2  See  Note  1. 
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the  strength  ratio  S  is  greater  than  one.  If  S  <  1,  Red  wins  in  finite  time.  If  S  =  1 
and  at  least  one  of  rbk  and  bbk  is  positive,  the  engagement  ends  at  finite  time  in  a 
tie;  otherwise,  the  engagement  continues  indefinitely,  with  nb  and  nr  tending  to 
zero.  The  fact  that  the  effectiveness  parameters  kb  and  kr  appear  linearly  in  the 
strength  ratio  S  while  the  sides’  numbers  appear  quadratically  is  called  the 
“Lanchester  square  law.”  That  law  has  been  used  to  argue,  not  always  convinc¬ 
ingly,  that  “quantity  trumps  quality”  for  the  strength  of  military  forces. 

If  the  Blue  side  wins  the  engagement,  its  losses  are 


b,-b  =  b„-  JA).  [Eq.C-4] 

V  kbb0 

Let  us  compare  the  results  of  analyzing  a  specific  M  vs.  N  engagement  with  the 
Lanchester  force-on- force  model  and  with  the  probabilistic  engagement  model. 
For  a  4  vs.  8  engagement  in  which  k  =  10,  both  sides  fighting  to  annihilation,  the 
Lanchester  model  predicts  that  the  Blue  side  will  win,  losing  0.902  aircraft  (inter¬ 
preting  fractional  aircraft  losses  is  one  difficulty  of  using  the  Lanchester  model). 
The  Reds  lose  8  aircraft. 

The  M  vs.  N  probabilistic  engagement  model  predicts  that  Blue  will  win  with 
probability  0.91,  with  the  expected  value  of  the  Blue  losses  equal  to  1.075  and  the 
expected  value  of  the  Red  losses  equal  to  7.595. 

The  expected  loss  values  are  not  too  different  from  those  of  the  Lanchester  result. 
Nevertheless,  the  information  that  Blue’s  victory  is  somewhat  in  doubt — Red  has 
nearly  a  10  percent  chance  of  winning — may  be  important  infonnation  for  plan¬ 
ners  to  have. 

The  distribution  of  Blue  losses,  including  the  probabilities  of  losing  0,  1,2,  3,  or  4 
Blue  aircraft,  shown  in  Figure  C-l,  is  also  infonnative.  There  is  a  43  percent 
probability  of  losing  no  aircraft,  a  29  percent  probability  of  losing  1  aircraft,  and  a 
28  percent  probability  of  losing  more  than  2  aircraft.  Thus  there  is  a  substantial 
probability  of  losses  almost  twice  the  expected  value  of  1.075.  This  sort  of  infor¬ 
mation  may  be  particularly  useful  to  loss-averse  planners. 

Certainly,  there  are  more  sophisticated  Lanchester  models  than  the  one  considered 
here.  There  are  also  more  sophisticated  computational  probabilistic  models  than 
the  one  considered  here.  We  believe  that  the  comparison  in  this  section  does, 
however,  indicate  the  relative  helpfulness  of  the  two  classes  of  models  to  military 
analysts. 
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Figure  C-1.  Distribution  of  Blue  Losses 


Blue  a/c  Lost 


Simulation  Modeling 


Simulation  models  have  significant  strengths.  They  enable  consideration  of  many 
factors  and,  like  probabilistic  engagement  models,  produce  outcome  statistics  that 
give  planners  insight  into  important  uncertainties.  However,  simulation  models 
also  have  defects  directly  corresponding  to  their  virtues.  Setting  up  and  operating 
large-scale  simulation  models  can  take  significant  time,  making  them  less  helpful 
to  planners  who  need  to  consider  wide  varieties  of  cases  quickly.  Moreover,  with 
many  inputs,  it  can  be  hard  to  tell  just  what  inputs  are  responsible  for  which  out¬ 
put  effects. 

Simulation  models  require  careful  consideration  of  the  numbers  of  runs  required 
to  establish  confidence  in  certain  results,  and  they  make  significant  demands  on 
random  number  generation.  Let  us  illustrate  these  points  with  an  example.  We 
consider  the  engagement  for  k  =  3.  Table  C-1  shows  the  probabilities  of  the  four 
outcome  states,  readily  computed  with  the  iterative  method  of  Appendix  A  for  the 
probabilistic  model. 


Table  C-1.  Outcome  State  Probabilities 


Outcome  state 

Probability 

(4,  0) 

0.6000 

(3,  0) 

0.2423 

(0,  3) 

0.0808 

(0,  4) 

0.0769 
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Figures  C-2  and  C-3  show  results  for  the  probability  of  state  (4,0)  and  for  the 
probability  of  state  (0,4),  respectively,  obtained  with  100-run  simulations. 

Figure  C-2.  P(4,  0)  from  100-Run  Simulations 


♦  P(4,  0)  by  Simulation  -Correct 


Figure  C-3.  P(0,  4)  from  100-Run  Simulations 
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♦  P(0,  4)  by  Simulation  —  "Correct" 

The  value  of  P(4,  0)  is  not  determined  very  well  with  100-run  simulations,  and  the 
value  of  the  rarer  outcome,  (0,  4),  is  detennined  badly.  These  results  for  such  a 
simple  case  may  serve  as  a  warning  to  users  of  simulations  who  pick  a  number  of 
runs  casually.  (“Oh,  a  hundred  runs  ought  to  be  enough.”) 

Either  analysis  or  experiment  show  that  at  least  1,000  runs  are  required  to  bring 
below  1 1  percent  the  coefficient  of  determination  of  values  of  P(0,  4)  inferred  by 
simulation.  Figure  C-4  shows  some  example  results. 


3 


Reference  for  analytic  determination  of  the  number  of  runs  required. 
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Figure  C-4.  P(0,  4)  from  1,000-Run  Simulations 


♦  P(0, 4)  by  Simulation  -Correct^ 


The  results  shown  in  Figures  C-2  through  C-4  were  obtained  with  a  random  num¬ 
ber  generator  that  we  believe  to  be  a  good  one.4  Careless  use  of  the  random  num¬ 
ber  generator  in  a  popular  application  gave  the  results  shown  in  Figure  C-5 

Figure  C-5.  P(0,  4)  from  1,000-Run  Simulations 
with  Poor  Random  Number  Generator 


♦  P(0, 4)  by  Simulation  -Correct 


Simulations  made  with  care  are  thoroughly  useful.  SLAACM  is  by  no  means  in¬ 
tended  to  compete  with  well-thought-out,  large-scale  simulations.  Rather,  we  be¬ 
lieve  it  is  a  useful  complement  to  such  models,  valuable  because  it  allows  rapid 
analysis  of  several  cases  in  the  PC  environment. 


4  It  is  the  routine  "rani"  described  by  W.H.  Press  and  others  in  Numerical  Methods  in  C,  Sec¬ 
ond  Edition  (Cambridge  University  Press,  1996),  p.  280. 
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Appendix  D 

Advanced  Probabilistic  Engagement  Models 


The  simple  M  vs.  N  probabilistic  engagement  model  offers  an  organized  approach 
to  modeling  the  statistics  of  air-to-air  engagements.  In  addition,  when  the  en¬ 
gagements  go  to  completion  (that  is,  to  the  long-time  limit),  it  is  straightforward 
to  obtain  its  one  parameter,  the  kill-rate  ratio,  k,  from  observations  of  loss  ratios 
obtained  from  combat  data,  training  exercises,  and  simulations.  However,  the  sim¬ 
ple  model  has  limitations.  Therefore,  in  this  appendix,  we  discuss  more  sophisti¬ 
cated  models  that  eliminate  some  of  those  limitations.  Straightforward 
modifications  to  SLAACM  will  allow  any  of  these  models  to  be  used  in  place  of 
the  M  vs.  N  probabilistic  engagement  model. 

Detailed  Engagement  Management 


One  limitation  of  the  M  vs.  N  probabilistic  engagement  model  is  the  assumption 
of  independent  events.  After  all,  the  Blue  and  Red  aircraft  are  attacking  fixed  sets 
of  airplanes,  and,  particularly  as  numbers  grow  small,  the  independence  may  seem 
unrealistic. 

One  way  to  remove  this  limitation  is  to  develop  models  with  more  specific  infor¬ 
mation  about  the  sides’  management  of  the  engagement.  For  example,  if  a  4  vs.  4 
engagement  splits  immediately  into  four  1  vs.  1  engagements,  the  independence 
assumption  is  more  realistic. 

Suppose  that  both  sides’  battle  management  has  the  4  vs.  4  engagement  split  into 
four  1  vs.  1  engagements,  each  of  which  goes  to  the  long-time  limit,  and  that  af¬ 
terwards  the  Red  and  Blue  survivors  depart  without  further  engagement.  Then 
there  are  five  outcome  states,  and  their  probability  distribution  is  binomial.  Spe¬ 
cifically, 


P(m,n)  =  B(m-n,  4,  k/(k  +  1)).  [Eq.  D-l] 

Equation  D-l  also  holds  when  the  time-to-kill  distributions  of  the  two  sides  are 
not  exponential;  in  such  cases,  the  quantity  k/(k  +  1)  is  replaced  by  the  probability 
that  the  Blue  side  makes  the  first  kill  in  a  1  vs.  1  engagement.  Figure  D-l  shows 
the  outcome  loss  distribution  when  k  =  3. 
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Figure  D-1.  Loss  Distribution:  4  vs.  4  Engagement 
Becomes  Four  1  vs.  1  Engagements 


If,  after  the  4  vs.  4  engagement  breaks  into  four  1  vs.  1  engagements,  each  side 
breaks  away  from  the  entire  engagement  upon  suffering  two  losses,  the  engage¬ 
ment  statistics  are  binary  through  the  first  two  kills.  If  these  result  either  in  two 
Blue  losses  or  in  two  Red  losses,  the  engagement  ends.  Otherwise,  the  first  two 
kills  have  left  three  aircraft  on  each  side,  and  the  engagement  ends  with  the  next 
kill.  Figure  D-2  shows  the  loss  distribution  for  this  battle  management  scenario, 
again  for  k  =  3. 

Figure  D-2.  Loss  Distribution:  Alternate  Battle  Management 


4 


Red  Loss 


Multiphase  Kill  Processes 


Another  concern  about  the  simple  M  vs.  N  probabilistic  engagement  model  is  that 
the  modes  of  the  assumed  exponential  time-between-kill  distributions  are  at  zero 
(Figure  D-3).  It  may  seem  unrealistic  to  have  zero  as  the  most  likely  time  between 
kills. 
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Figure  D-3.  Exponential  Distribution 


Time  in  units  of  mean  kill  time 

This  concern  may  be  addressed  by  modeling  the  kill  process  as  taking  place  in 
two  phases.  In  the  first  phase,  a  combatant  seeks  opponents.  This  phase  ends 
when  the  combatant’s  fire  control  system  locks  onto  an  opponent. 

The  second  phase  begins  with  lock-on  and  ends  when  the  combatant  has  fired  a 
successful  missile.  The  assumption  of  exponential  distributions  of  the  times  re¬ 
quired  for  these  two  phases  is  not  as  heavy  as  the  assumption  of  an  exponential 
distribution:  with  effective  guidance  from  C3ISR  systems,  a  combatant’s  lock-on 
might  occur  essentially  immediately  after  the  search  begins,  and,  with  a  firm  lock 
and  an  effective  missile,  firing  of  a  successful  missile  might  occur  essentially 
immediately  after  lock-on. 

The  total  time-between-kills  distribution  for  a  two-phase  process  is  the  convolu¬ 
tion  of  two  exponentials.  (It  is  not  always  the  Erlang  distribution;  that  happens 
only  when  the  parameters  of  the  two  exponential  distributions  are  the  same.)  Fig¬ 
ure  D-4  shows  an  example. 

Figure  D-4.  Time-between-Kills  Distribution  for  Two-Phase  Kill  Process 


Time  to  kill,  arbitrary  units 
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The  engagement  state  with  two-phase  kill  processes  requires  more  infonnation  to 
describe  it  than  the  two  numbers  required  when  single -phase  kills  were  used.  We 
use  the  state  variable  (m,  i;  n,  j),  where  m  is  the  number  of  Blue  aircraft  in  the  en¬ 
gagement,  i  is  the  number  of  Blues  that  have  locked  onto  targets,  n  is  the  number 
of  Red  aircraft  in  the  engagement,  and  j  is  the  number  of  Reds  that  have  locked 
onto  targets.  Instead  of  the  two-valued  state  transition  diagram  of  Figure  D-l,  the 
state  transition  diagram  for  the  two-phase  kill  process  has  six  elements,  as  shown 
in  Figure  D-5.  Text  on  the  arrows  in  the  figure  shows  the  rates  for  the  respective 
transitions. 

Figure  D-5.  State  Transition  Diagram,  Two-Phase  Kill  Process 


(m,  i  +  1;n,j)  Blue  makes  lock 


In  discussing  multiphase  kill  models  and  other  advanced  probabilistic  engagement 
models,  we  will  use  dimensioned  parameters,  because  that  seems  to  make  the 
sense  of  graphs  and  equations  easier  to  grasp.  In  Figure  D-5,  kbi  is  the  parameter 
of  the  exponential  distribution  of  the  time  required  for  Blue  aircraft  to  lock  on  to 
an  opponent,  and  kb2  is  the  parameter  of  the  exponential  distribution  of  the  time 
required  for  a  Blue  aircraft  to  make  a  kill  after  lock-on.  Parameters  kaJ  and  ka2  are 
the  homologues  of  kbi  and  kb2,  respectively,  for  the  Red  aircraft.  In  computing 
loss  distributions,  however,  we  non-dimensionalize  the  evolution  equations  with 
kir,  so  that  only  the  three  ratios  kbi/kri,  kb2/kri,  and  kr2/kri  are  used,  and  we  con¬ 
sider  only  long-time  limits.  Also  in  Figure  D-5,  no  aircraft  can  track  while  scan¬ 
ning,  so  the  rates  at  which  new  locks  are  made  is  (m  -  i)kbi  and  (n  -  j)kri. 

The  two-phase  kill  process  allows  much  more  detailed  modeling  of  air-to-air  en¬ 
gagements  than  does  the  single-phase  process.  Effects  of  track- while-scan,  of 
multiple-target  tracking,  and  of  sharing  targeting  infonnation  among  the  elements 
of  a  flight  can  be  modeled.  Figure  D-6  shows  how  the  state  transition  diagram 
changes  from  the  one  of  Figure  D-5,  when  elements  of  the  Blue  flights  can  track 
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while  scanning  and  share  targeting  information  among  flights.  The  state  descrip¬ 
tion  for  Figure  D-6  is  (m,  i;  n,  j),  where  now  i  is  the  total  number  of  targets 
tracked  by  the  Blue  flight. 

Figure  D-6.  State  Transition  Diagram,  Track-While-Scan 
and  Information  Sharing  for  Blue  Aircraft 


Blue  acquires  a  target  Blue  kills  a  tracking  red 

(m,  i+1;n,j)  (m,  i-1;  n-1,  j -1)  Blue  kills  a  non-tracking  red 


m  =  number  of  blues,  i  =  number  of  targets  tracked  by  blue  flight;; 
n  =  number  of  reds,  j  =  number  of  tracking  reds 

The  two-phase  kill  model  may  be  extended  to  account  for  a  tracked  aircraft  break¬ 
ing  lock,  as  shown  by  the  state  transition  diagram  of  Figure  D-7. 

Figure  D-7.  State  Transitions  with  Lock-Breaking 


(m,  i  ;  n,  j  - 1) 
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The  evolution  equations  for  the  state  probabilities  are  necessarily  more  complex 
for  two-phase  kill  processes.  For  no  track-while-scan  or  sharing  of  targeting  in¬ 
formation,  they  are 


=  -[(m  -  Okb!  +  ikb2  +  (n  -  j)kr,  +  jkr2  ]Pm4;nJ 
+  [i]+(m-i  +  l)kb1Pmi_1;nj 

+  (i  +  l)kb2Pm4+u+1>j-^lTJ 

+  (i  +  l)kb2Pmi+ln+1  j+1 

+  [j]+(n-j  +  l)kr1Pminj_1 

m  +  l-i 

+  (j  +  l)kr2Pm+1  i;n  j+1  •  — — 

+  (j  +  l)kr2Pm+li+ln  j+1 


[Eq.  D-2] 


The  symbol  [x]+  denotes  the  function  of  x  that  equals  1  if  x  >  0,  and  0  otherwise. 

Though  more  complex  than  the  differential  equations  for  the  evolution  of  state 
probabilities  with  single-phase  kill  processes,  Equation  D-2  is  still  a  system  of 
linear  ordinary  differential  equations  with  constant  coefficients.  The  iterative 
scheme  of  Appendix  A  still  allows  one  to  evaluate  the  long-time  limit  of  absorb¬ 
ing  boundary  states  with  a  finite-step  iterative  procedure,  and  NASA’s  ASSIST 
and  STEM  will  treat  large-dimensional  cases  effectively. 

Loss  statistics  for  long-term  limits  of  engagements  with  the  two-phase  kill  model 
do  not  look  different  from  those  with  the  one-phase  kill  model.  They  do,  however, 
fit  outcomes  of  certain  large-scale  simulations  better.  Table  D-l  shows  an  exam¬ 
ple  in  which  the  three  adjustable  rate  ratios  of  the  two-phase  kill  model  were  fitted 
to  give  a  match  of  the  four  outcome  probabilities  of  an  EADSIM  simulation.  The 
engagement  was  4  vs.  4,  with  both  sides  breaking  at  two  losses.  The  fitted  rate 
ratios  were  kbi/kri  =  4.94,  kb2/kri  =  4.90,  and  kr2/kri  =  1.00. 

Table  D-1.  Loss  Probabilities  from  EADSIM  and  from  Two-Phase  Kill  Model 


Model 

(0,2) 

(1,2) 

(2,  1) 

(2,0) 

EADSIM 

0.886 

0.10 

0.0084 

0.0047 

Two-Phase  Kill 

0.885 

0.10 

0.01 

0.005 

Figure  D-8  displays  the  outcome  probabilities. 
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Figure  D-8.  Loss  Probabilities  from  Two-Phase  Kill  Model 


We  were  encouraged  to  find  that  not  only  did  the  outcome  probabilities  of  the 
two-phase  kill  model  match  all  four  of  the  EADSIM  simulations  within  about  16 
percent,  but  the  time  evolution  of  the  absorbing  boundary  states  also  agreed  rea¬ 
sonably  well  with  those  of  EADSIM.  Figures  D-9  and  D-10  show  this  agreement. 

Figure  D-9.  Time  Evolution  of  Probability  of  Loss  State  (0,2) 

[0,  2] 
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Figure  D-10.  Time  Evolution  of  Probability  of  Loss  State  (1,2) 


[1.2] 


Since  the  initial  time  and  time  scale  are  arbitrary  in  EADSIM,  we  adjusted  those 
two  parameters  to  give  the  result  in  Figure  D-10,  but  made  no  further  adjustments. 

Missile-Tracking  Models 


The  probabilistic  models  that  we  have  considered  so  far  implicitly  assume  that  the 
combatants  never  run  out  of  missiles.  Missile  loads  may  limit  aircraft  effective¬ 
ness,  of  course,  and  we  have  made  probabilistic  models  that  account  for  this.  In 
this  section,  we  describe  a  probabilistic  engagement  model  in  which  each  aircraft 
is  loaded  with  6  missiles. 

An  engagement  state  description  for  this  case  is  (m0,  mi,  ...,  m6;  no,  ni,  ...,  n6), 
where  rrij  denotes  the  number  of  Blue  aircraft  with  i  missiles,  and  n,  gives  the 
number  of  Red  aircraft  with  j  missiles. 

Table  D-2  gives  the  end  states  and  transition  rates  for  this  model,  when  the  start¬ 
ing  state  is  (m0,  mi,  ...,  m6;  no,  ni,  ...,  n6).  In  the  table,  ntot  is  the  total  number  of 

6 

Red  aircraft,  given  by  ^rij  ,  and  mtot,  computed  analogously,  is  the  total  number 
j-<> 

of  Blue  aircraft.  The  quantity  Ub  =  (1  -  PkbVPkb,  and  ur  =  (1  -  PkrVPkr- 
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Table  D-2.  End  States  and  Transition  Rates,  Missile-Tracking  Model 


Event 

End  state 

Transition 

rate 

Blue  m,  kills  a  Red  nj, 

1  <  i  <  6,  0  <j  <  6 

(m0,  . . . ,  mM  +  1 ,  m,  -1 ,  ...  m6;  n0 . nj  - 1 ,  . . . ,  n6) 

rrii  kb  n/ntot 

Blue  m,  misses, 

1  <i<6 

(m0,  ...,  mM  +  1,  ny- 1,  ...;  n0,  n! . n6) 

mi  kb  ub 

Red  n,  kills  a  Blue  m,, 

1  <j  <  6,  0  <  i  <  6 

(m0,  ...,  171,-1,  ...,  m6;  n0,  ...,  n^  +  1,  nr1,  ..., 
n6) 

^  krmi/mtot 

Red  nj  misses, 

1  <j  <6 

(m0,  ...,  m6;  n0,  ...,  nj.!  +  1,^-1,  ...,  n6) 

IT  kr  ur 

The  evolution  equations  for  the  state  probabilities  of  the  missile-tracking  model 
are  again  a  system  of  linear  ordinary  differential  equations  with  constant  coeffi¬ 
cients,  and  again  the  iterative  scheme  of  Appendix  A  may  be  used  to  compute 
end-state  probabilities  with  finite-step  iterations.  As  Table  D-2  suggests,  however, 
the  number  of  equations  is  large;  tens  of  thousands  of  transient  states  would  be 
tracked  in  a  complete  model  of  a  4  vs.  4  engagement  in  which  each  aircraft  car¬ 
ried  6  missiles,  and  aircraft  left  the  engagement  only  by  being  destroyed.  We  usu¬ 
ally  treat  missile  tracking  models  with  NASA’s  ASSIST  and  STEM.  Appendix  B 
describes  this  method. 

Loss  probability  distributions  from  missile-tracking  models  have  the  same  general 
appearance  as  those  from  other  probabilistic  models. 

Radar  Duel  Models 


Engagements  may  not  begin  with  all  combatants  mutually  visible  by  radar. 
Rather,  they  may  start  with  each  side,  acting  on  guidance  from  their  own  battle 
management,  looking  for  opponents.  The  side  with  the  more  effective  radar,  and 
the  smaller  radar  cross  section,  may  detect  opponents  before  opponents  detect 
them.  In  this  case,  they  may  launch  missiles  and,  seeing  these  missiles’  effects, 
decide  whether  or  not  to  continue  the  engagement.  We  call  this  kind  of  engage¬ 
ment  a  “radar  duel.” 

In  our  models  of  radar  duel  engagements,  the  engagement  takes  place  in  two 
phases.  In  the  first  phase,  the  combatants  are  given  guidance  from  their  side’s 
C4ISR  and  battle  management  that  is  equivalent  to  an  interval  of  azimuth  A0  and 
an  interval  <E>A,  in  which  to  look  for  their  opponents.  This  phase  ends  when  one 
side  first  detects  its  opponents. 

If  the  side  that  first  detects  its  opponents  has  time  to  fire  one  or  more  rounds  of 
missiles  and  observe  their  effects  before  being  itself  detected,  that  side  may  de¬ 
cide  either  to  continue  the  engagement  or  to  break  it  off.  The  second  phase  ends 
when  the  action  after  the  first  detection  is  completed. 
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Ratios  of  Detection  Ranges 

1  2 

The  statistics  of  detection  by  pulse-Doppler  radars  are  complex.  ’  Nevertheless, 
fairly  straightforward  considerations,  which  we  summarize  here,  can  relate  ratios 
of  detection  ranges  to  properties  of  radar  systems  and  to  battle  management. 

The  signal  energy  Es  deposited  in  the  receiver  of  a  scanning  radar  from  a  target  of 
radar  cross  section  a  at  range  R  in  a  given  azimuth-elevation  box  is 


E 


S 


PG  a  A 
4ttR2  4ttR2  L  tbox' 


[Eq.  D-3] 


where  P  is  the  radar’s  average  power,  G  is  the  gain  of  its  transmitting  antenna,  A 
is  the  electrical  aperture  of  its  receiving  antenna,  and  the  factor  L  accounts  for 
losses.  The  time  tbox  is  the  equivalent  time  that  the  radar’s  main  lobe  dwells  in  the 
given  azimuth-elevation  box. 

The  time  W  is  equal  to  the  total  time  tscan  required  to  cover  the  search  region  es¬ 
tablished  by  the  guidance  parameters  A0  and  AO,  divided  by  the  number  of  the 
radar’s  azimuth-elevation  boxes  required  to  cover  the  search  region: 

^box  —  ^scan  /^box  •  [Eq-  D-4] 


The  parameter  nbox  will  be  given  roughly  by  the  product  of  the  number  of  azi¬ 
muthal  beam-widths  in  the  azimuth  guidance  interval  A0  with  the  number  of  ele¬ 
vation  beam- widths  in  the  elevation  guidance  interval  AO.3’4  Thus 


_  A0  AO 
a  p 


[Eq.  D-5] 


where  a  is  the  radar’s  azimuthal  beam- width  and  P  is  its  elevation  beam- width. 5 


The  time  tscan  is  bounded  above  by  the  need  to  complete  a  scan  before  the  range 
has  changed  by  a  significant  fraction  of  the  size  of  a  range  box  in  the  pulse- 
Doppler  radar  system. 


1  Many  standard  references  treat  this  subject.  See,  for  example,  G.W.  Stimson,  Introduction  to 
Airborne  Radar  (El  Segundo,  CA:  Hughes  Aircraft  Company,  1983)  and  M.l.  Skolnik,  Introduc¬ 
tion  to  Radar  Systems,  Third  Edition  (New  York:  McGraw-Hill,  2001). 

2  P.  Swerling,  "Probability  of  Detection  for  Fluctuating  Targets,"  IRE  Transactions  on  Infor¬ 
mation  Theory,  Volume  6,  April  1960,  pp.  269-308. 

3  This  assumes  that  battle  management's  guidance  is  good  enough  to  position  the  target  in  a 
region  that  the  radar  can  cover,  without  turning  the  aircraft.  If  that  is  not  so,  then  detection  times 
increase  greatly. 

4  Actual  radars  may  provide  fixed  choices  of  angular  box  widths.  For  example,  the 
AN/APG-66  provides  scans  of  ±10°,  ±20°,  or  ±30°,  with  1,  2,  or  4  elevation  bars. 

5  See  Note  3. 
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Thus  the  signal  energy  available  in  the  receiver  is 

_  PGgA  a  p 

s~(47rR2)2LtsCanA0A®  ' 


[Eq.  D-6] 


Noise  energy  from  various  sources  competes  with  signal  energy.  It  is  customary 
to  characterize  all  the  sources  of  noise  with  a  single  noise  temperature  T,  so  that 
the  noise  power  continuously  present  in  the  system  is  kTB,  where  B  is  the  system 
bandwidth  in  Hz  and  k  is  the  Boltzmann  constant,  1.38  x  10‘23  Joule/°K.  The 
bandwidth  B  needed  for  pulse-Doppler  processing  is  close  to  the  inverse  of  the 
time  on  target,  l/tot.6  Thus  the  noise  energy  deposited  over  time  tot,  which  com¬ 
petes  with  the  signal  energy  Es,  is  kT,  and  the  signal-to-noise  ratio  SNR  is 


SNR 


PGgA  a  p 

(47iR2)2kTL  scan  A0  AO 


[Eq.  D-7] 


Actual  detection  is  probabilistic,  and  involves  consideration  of  real  and  false  tar¬ 
gets.  A  simple  model,  which  nevertheless  captures  several  features  of  interest  for 
the  radar  duel,  is  to  assume  that  detection  is  certain  when  SNR  exceeds  a  critical 
value  SNR*.  With  this  assumption,  detection  is  certain  when  range  R  satisfies 


R  < 


PGgA  a  p 
SNR*  kTL  A0  AO  scan 


1 

4 


[Eq.  D-8] 


It  follows  that  the  ratio  of  detection  ranges  for  two  combatants  a  and  b  is  given  by 

Rb 

While  the  detection  range  ratio  of  Equation  D-9  was  derived  fairly  roughly,  it 
does,  we  believe,  give  useful  quantitative  means  of  detennining  which  side  is 
likely  to  win  the  first  phase  of  a  radar  duel,  from  specific  information  about  the 
two  sides’  radars  and  C4ISR/battle  management. 

As  an  example,  suppose  that  the  power-aperture  product  PA  for  the  Blue  side  is 
10  times  that  of  the  Red  side,  and  that  the  ratio  of  the  Blues’  radar  cross  section  to 
that  of  the  Red’s  is  0.01,  while  all  other  radar  and  battle  management  parameters 
are  the  same.  Then  the  ratio  of  Blue’s  detection  range  to  Red’s,  according  to 
Equation  D-9  ,  is  5.6  to  1.  Suppose  that  Blue’s  detection  range  is  56  nm.  Then 
Red’s  detection  range  is  approximately  9  nm. 


Pa  SNRb  Ga  Tb  Gb  Aa  Lb  aa  Pa  A@b  AQb  tscana 
Pb  SNRa  Gb  Ta  Ga  Ab  La  ab  Pb  A0a  AOa  tscanb 


[Eq.  D-9] 


6  M.I.  Skolnik,  Introduction  to  Radar  Systems,  Third  Edition  (New  York:  McGraw-Hill, 
2001). 
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Let  the  sides  each  have  4  aircraft  and  be  closing  at  600  kt.  Suppose  that,  on  de¬ 
tecting  their  adversaries  at  56  nm,  the  Blues  develop  a  firing  solution  and  launch  4 
missiles,  requiring  30  seconds  to  do  so.  When  Blue  launches,  the  flights  will  be 
5 1  nm  apart. 

Let  Blue’s  missiles  fly  at  Mach  2,  1,200  kt  at  the  altitude  flown.  If  the  Reds  are 
moving  toward  the  missiles  at  300  kt,  the  closing  rate  of  the  missiles  and  the  Red 
flight  will  be  1,500  kt.  The  missiles  will  reach  their  targets  in  a  little  over  2  min¬ 
utes,  when  the  opponents  are  3 1  nm  apart,  still  well  outside  Red’s  detection  range. 
If  surviving  Reds  continue  toward  the  Blues,  closing  at  600  kt,  the  Blues  have  2 
minutes  to  decide  whether  or  not  to  engage  them. 

Let  us  continue  this  example  to  determine  the  distribution  of  losses  in  the  en¬ 
gagement.  Suppose  that  the  4  missiles  that  Blue  launched  at  the  end  of  the  radar 
duel  phase  have  a  single-shot  kill  probability  of  0.3.  Then  Red’s  losses  from  this 
phase  have  a  binomial  distribution,  with  outcomes  shown  in  Table  D-3. 

Table  D-3.  Red  Loss  Distribution  from  Radar  Duel  Phase 


Red  aircraft  killed 

Probability 

4 

0.0081 

3 

0.0756 

2 

0.2646 

1 

0.4116 

0 

0.2401 

Further  suppose  that  the  kill-rate  ratio  kb/kr  is  2,  and  that  the  Blues’  rules  of  en¬ 
gagement  are  to  continue  to  a  close  engagement  only  if  the  probability  of  their 
sustaining  no  losses  while  defeating  all  their  opponents  is  greater  than  70  percent, 
and  that  they  will  break  off  the  near-range  engagement  on  sustaining  2  losses.  The 
M  vs.  N  probabilistic  engagement  model  described  above  gives  the  probabilities 
for  zero  Blue  losses  shown  in  Table  D-4.  In  view  of  those  probabilities,  the  Blues 
will  continue  the  engagement  only  if  the  Reds  lose  two  or  more  aircraft  in  the  ra¬ 
dar  duel  phase. 


Table  D-4.  Probabilities  for  Zero  Blue  Losses 


Engagement 

P  (no  Blue  losses) 

4  vs.  4 

0.345 

4  vs.  3 

0.517 

4  vs.  2 

0.711 

4  vs.  1 

0.889 
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Modeling  those  close  engagements  that  occur  with  the  M  vs.  N  probabilistic  en¬ 
gagement  model,  and  considering  losses  from  both  the  radar  duel  and  close- 
engagement  phases,  we  find  the  loss  distribution  of  Figure  D-l  1  for  this  radar 
duel  engagement.  In  the  figure,  we  show  Blue  win  states  with  blue  bars  and  Red 
win  states  with  Red  bars.  For  this  example,  the  outcome  (0,  0)  is  possible.  We 
show  the  bar  for  the  “tie”  outcome  in  purple.  We  also  show  the  bar  for  the  out¬ 
come  (0,  1)  in  purple,  since  there  seems  to  be  no  obvious  winner  in  that  case,  ei¬ 
ther. 


Figure  D-11.  Loss  Distribution  for  Radar  Duel  Example 


Model  for  Low-Observable  Aircraft 

If  the  probability  that  Red  radars  can  observe  a  Blue  aircraft  is  negligibly  small,  it 
is  reasonable  to  model  that  Blue  type  as  invulnerable:  it  is  never  destroyed  by  en¬ 
emy  action  in  combat.  Invulnerable  is  not  the  same  as  invincible,  however,  be¬ 
cause  the  Blue  aircraft  still  carry  only  a  finite  number  of  missiles,  with  single-shot 
kill  probabilities  less  than  one. 

Engagements  with  such  aircraft  against  vulnerable  foes  may  be  diagrammed  and 
their  outcome  probabilities  computed  straightforwardly.  Figure  D-12  shows  an 
example  when  invulnerable  aircraft  with  a  total  of  6  missiles  confront  4  vulner¬ 
able  Blue  aircraft  (it  does  not  matter  how  the  missiles  are  distributed  among  the 
Blue  aircraft).  The  engagement  state  is  the  pair  (i,  j)  where  i  is  the  number  of 
enemies  and  j  the  number  of  missiles.  The  parameter  p  is  the  single-shot  kill  prob¬ 
ability  of  the  Blues’  missiles,  and  q  is  1  -  p. 
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Figure  D-12.  Diagram  of  Engagement  with  Invulnerable  Blues 


Fire  until  either  all  targets  killed, 
or  no  missiles  left. 


B(0,6,p) 


3,0 


2,0 


1,0  B(4,6,p)-qB(4,5,p) 


B(1,6,p) 


B(2,6,p)  B(3,6,p) 


In  general,  if  the  starting  state  of  such  an  engagement  is  (e,  m),  and  e  <  m,  the  out¬ 
come  probabilities  are 

P(0 ,  m  -  e  -  j)  =  B(e,  e  + ),  p)  -  qB(e,  e  +  j  - 1,  p),  1  <  j  <  m  -  e 
P(0,m-e)  =  B(e,m,p)  .  [Eq.  D-10] 

P(j/0)  =  B(e  -  j,m,p) 

If  m  <  e,  all  m  missiles  will  of  course  be  fired,  and  the  outcome  probabilities  are 
P(e  -  j,0)  =  B(j,m,p),  0  <  j  <  m .  [Eq.  D-l  1] 

Figure  D-l 3  shows  a  loss  distribution  for  an  engagement  with  invulnerable  Blue 
aircraft.  The  bars  corresponding  to  Red  win  states  are  shown  in  red.  Although  no 
Blues  have  been  lost,  they  still  did  not  eliminate  their  opponents  before  running 
out  of  missiles. 


D-14 


Advanced  Probabilistic  Engagement  Models 


Figure  D-13.  Loss  Distribution,  Invulnerable  Blues 
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Applications  to  Analysis 


The  models  considered  in  this  appendix  offer  means  of  analyzing  the  impacts  of 
different  features  of  combat  aircraft  and  missiles  in  a  consistent  framework.  As  an 
example  of  this,  we  give  here  a  consistent  treatment  of  the  effects  of  radars,  com¬ 
bat  systems,  and  missile  single  shot  kill  probabilities. 

We  consider  a  two-phase  kill  model,  in  which  time-to-lock,  a  measure  of  radar 
effectiveness,  is  exponentially  distributed  with  parameter  ki,  and  time-to-fire  after 
lock,  a  measure  of  combat  system  effectiveness,  is  exponentially  distributed  with 
parameter  ki.  We  denote  the  missiles’  single-shot  kill  probability  by  p. 

We  assume  that  the  firing  discipline  is  lock-shoot-break  lock-lock-shoot...,  so 
that  the  probability  q  of  missing  directly  influences  the  time  to  make  a  kill.  This  is 
different  from  the  effect  of  p  in  a  single-phase  kill  process.  There,  p  influences 
only  the  rate  of  firing  required  to  give  a  specified  kill  rate. 

The  mean  time  to  make  lock  is  1/ki,  the  mean  time  to  fire  after  making  lock  is 
l/k?,  and  so  the  mean  time  to  fire  is  1/ki  +  l/k?.  The  probability  of  making  a  kill 
with  one  shot  is  p,  the  probability  of  making  a  kill  with  2  shots  is  qp,  where 
q  =  1  -  p,  and  the  probability  of  requiring  n  shots  to  make  a  kill  is  qn  lp.  Accord¬ 
ingly,  T  ,  the  expected  time  to  make  a  kill,  is 


[Eq.  D-12] 
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Since 


00  11 

Zjqrl=7r-T  =  — .  [Eq.D-13] 

i-i  0-q)  p 

we  have 

T  =  —  — +  — 

PLk!  k2 

where  TL  is  the  mean  time  to  make  lock,  and  TF  is  the  mean  time  to  fire  after 
making  lock. 

Thus  all  combinations  of  p,  TL  ,  and  TF ,  for  which  the  quantity  ( TL  +  TF  )/p  has  a 
given  value,  give  the  same  mean  time  to  kill,  and  Equation  D-14  provides  a  quan¬ 
titative  indication  of  the  way  in  which  radar  system  effectiveness,  combat  system 
effectiveness,  and  missile  effectiveness  can  be  traded. 

A  Comment  on  Distributions 


=1[tl+tf; 


[Eq.  D-14] 


We  have  used  exponential  distributions  in  most  of  our  models.  It  is  certainly  ap¬ 
propriate  to  note  that  these  distributions  assign  positive  probabilities  to  arbitrarily 
long  times,  and  that  this  is  not  physically  realistic  because  engagements  do  not 
last  more  than  a  few  minutes.  Nevertheless,  we  believe  that  the  exponential  distri¬ 
butions  are  useful,  because  their  logical  basis — that  the  probability  of  an  event 
occurring  in  a  sufficiently  small  interval  of  time  is  proportional  to  the  size  of  the 
interval,  and  independent  of  the  starting  time  of  the  interval,  and  the  probability  of 
simultaneous  events  is  negligibly  small — reflects  events  in  the  important  phases 
of  air-to-air  engagements.  We  would  reject  conclusions  that  involved  unrealisti¬ 
cally  long  times  in  any  significant  way. 

The  situation  is  analogous  to  statistical  mechanics.  There,  too,  distributions  with 
unbounded  support  are  used,  where  extreme  values  of  the  support  are  unrealistic. 
This  leads  to  such  paradoxes  as  positive,  but  very  small,  probabilities  that  all  the 
molecules  in  a  room  would  find  themselves  in  one  half  of  the  chamber,  so  that 
people  in  the  other  half  would  gasp  for  air.  Applied  scientists  dismiss  these  pre¬ 
dictions,  just  as  we  would  dismiss  predictions  involving  unrealistically  long  times. 

Over  the  last  40  years,  an  extensive  literature  on  stochastic  duels  has  developed, 
treating  engagements  (usually  with  very  small  numbers  of  combatants)  with  arbi¬ 
trary  time-to-kill  distributions.7  If  it  is  essential  to  use  distributions  with  bounded 
support,  SLA  ACM  can  be  adapted  to  use  stochastic  duel  results. 


7  See,  for  example,  A.V.  Gafarian  and  C.J.  Ancker,  "The  Two-On-One  Stochastic  Duel,"  Na¬ 
val  Research  Logistics  Quarterly,  Vol.  31,  1984,  pp.  309-324. 
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Appendix  E 

Accounting  for  Electronic  Attack 


The  method  of  Appendix  A  gives  a  readily-implemented  way  to  allow  present 
SLAACM  versions  to  include  impacts  of  “deception”  electronic  attacks  (EAs).  A 
two-phase  kill  engagement  model  with  a  “break-lock”  transition  gives  a  natural 
way  to  characterize  the  impacts  of  such  attacks.  Figure  E-l  diagrams  a  1  vs.  1  en¬ 
gagement,  with  two-phase  kill  and  break-away. 1  To  make  a  kill,  each  combatant 
must  make  a  radar  lock  on  his  opponent;  having  done  so,  he  must  fire  a  (success¬ 
ful)  missile.  A  combatant  who  perceives  that  his  opponent's  radar  has  locked-on 
may  use  deception  electronic  attack  to  break  the  lock. 

The  system  state  is  (m,  i;  n,  j),  where  m  is  the  number  of  Blue  aircraft,  i  the  num¬ 
ber  of  Blues  having  lock;  n  is  the  number  of  Red  aircraft,  and  j  is  the  number  of 
Red  aircraft  having  lock.  Labels  on  the  arrows  of  Figure  E-l  show  the  rates  (not 
probabilities)  for  each  transition. 

The  parameter  kbl  is  the  reciprocal  of  the  mean  time  for  Blue  to  make  lock,  and 
kb2  is  the  reciprocal  of  the  mean  time  for  Blue  to  launch  a  successful  missile, 
having  made  lock.  Parameter  kbx  is  the  reciprocal  of  the  mean  time  for  a  Blue 
aircraft  to  break  a  Red  lock.  The  parameters  krl,  kr2,  and  krx  have  the  homolo¬ 
gous  meanings  for  Red. 

Figure  E-1.  Diagram  of  Two-Phase  Kill  Engagement  With  Break-Lock 


1  The  discussion  of  the  1  vs.  1  engagement  with  breakaway  includes  material  from  Appendix 
A,  given  here  for  the  sake  of  completeness  and  ease  of  reading. 
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For  comparison,  Figure  E-2  shows  a  diagram  of  a  1  Vs  1  engagement  with  a  sin¬ 
gle-phase  kill  mechanism.  There,  the  state  is  simply  (number  of  Blue  a/c,  number 
of  Red  a/c). 


Figure  E-2.  Diagram  Of  Single-Phase  Kill  Engagement 


Obviously  from  Figure  E-2,  for  the  single-phase  kill  engagement  the  probability 
of  a  Blue  win  is  equal  to  kb/(kb  +  kr),  and  the  probability  of  a  Red  win  is  kr/(kb  + 
kr).  The  ratio  of  the  probability  of  a  Blue  win  to  the  probability  of  a  Red  win  is,  in 
this  case,  equal  to  the  kill-rate  ratio  kb/kr. 

This  fact  suggests  inferring  an  equivalent  kill-rate  ratio  from  more  complicated  1 
vs.  1  engagements,  as  the  ratio  of  the  probability  of  a  Blue  win  to  the  probability 
of  a  Red  win.  We  now  develop  equations  to  do  that,  for  the  two-phase  kill  en¬ 
gagement  with  break-lock,  diagrammed  in  Figure  E-l. 

By  referring  to  Figure  E-l,  one  may  write  down  the  evolution  equations  (forward 
Chapman-Kohnorogov  equations)  for  the  engagement.  Let  Pioio(t)  denote  the 
probability  that  the  system  is  in  state  (1,0;  1,  0)  at  time  t,  and  let  Pi  no,  Pion,  and 
Pi  in  be  the  probabilities  that  the  system  is  in  the  state  indicated  by  the  subscripts. 
Then,  representing  the  probability  of  a  Blue  win  at  time  t  by  Pnw(t)  and  the  prob¬ 
ability  of  a  Red  win  at  time  t  by  Prw(0,  the  evolution  equations  for  the  1  Vs  1  en¬ 
gagement  with  break-lock  are 
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P1010  _  — (kbl  +  krl)P1010  +  krxP1110  +  kxbP1011 

Pi  1 10  =  -(kb2  +  krl  +  krx)pi  110+  kblPioio  +  kbxPi  1 1 1 
P1011  =  -(kr2  +  kbl  +  kbx)P101 1  +  krlP1010  +  krxPj ,  j  t 

[Eq.  E-l] 

pim  =-(kb2  +  kr2  +  kbx  +  krx)Plin  +krlP1110  +kblP1011 
pbw  =  kb2(P1110  +P11U) 
prw  =  kr2(Pion  +Pllll) 

In  Appendix  A  we  explain  how  to  calculate  the  long-time  limiting  values  of  Pbw 
and  Prw  by  integrating  equations  F-l  from  zero  to  infinity,  using  the  condition 
that  all  the  probabilities  are  zero  at  t  =  0  except  Pioio,  which  has  the  value  1.  For 
completeness,  we  repeat  that  development  for  the  case  at  hand.  Defining  xi  as  the 
integral  over  all  positive  time  of  Pioio,  *2  as  that  same  integral  of  Pi  no,  X3  as  that 
integral  of  P 1011,  and  X4  as  that  integral  of  Pun,  one  finds  that  the  Xj  are  given  by 
the  solution  of  a  system  of  four,  linear  algebraic  equations,  the  augmented  matrix 
of  which  is 


f-  (kbl  +  krl) 
kbl 
krl 
0 


krx 

-(kb2  +  krl  +  krx) 
0 

krl 


kbx 

0 

-  (kr2  +  kbl  +  kbx) 
kbl 


0  -f 

kbx  0 

krx  0 

-  (kb2  +  kr2  +  kr  x  +  kbx)  0 


V 


J 


One  also  finds  that 


lim  PBW  =  kb2(x2  +x4) 
t  — »  00 

lim  PRW  =kr2(x3 +x4) 

t  — »  00 


[Eq.  E-2] 


[Eq.  E-3] 


Given  values  of  kbl,  kb2,  kbx,  krl,  kr2,  and  krx,  the  equivalent  kill-rate  ratio 
given  by  Pbw/Prw  follows  directly  from  Equations  E-2  and  E-3.  Varying  the  pa¬ 
rameters  kbx  and  krx  then  gives  means  of  inferring  the  impacts  of  deception  EA 
on  kill-rate  ratios. 


Numerical  experiments  suggest  that  this  impact  is,  in  fact  (and  perhaps  fortui¬ 
tously),  fairly  straightforward.  Figure  E-3  shows  the  variation  of  equivalent  kill- 
rate  ratios  with  the  ratio  of  kbx  to  krx,  for  two  cases.  In  the  first,  shown  with  tri¬ 
angle  plot  points,  the  mean  time  for  Blue  to  make  lock  is  20  seconds,  the  mean 
time  for  Blue  to  make  a  kill  given  lock  is  three  seconds,  and  the  mean  time  for 
Blue  to  break  lock  is  one  second;  the  mean  time  for  Red  to  make  lock  is  60  sec- 
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onds,  and  the  mean  time  for  Red  to  make  a  kill  given  lock  is  15  seconds.  In  the 
second  case,  shown  with  square  plot  points,  the  Blue  parameters  are  unchanged, 
while  the  mean  time  for  Red  to  make  lock  is  reduced  to  20  seconds,  and  the  mean 
time  for  Red  to  make  a  kill  given  lock  is  three  seconds. 

Figure  E-3.  Variation  Of  Kill-Rate  Ratio  with  Ratio  Of  Break-Lock  Parameter 


Ratio  of  Mean  time  for  Red  to  break  Blue  lock  to  mean  time 
for  Blue  to  break  Red  lock 

In  both  these  distinctly  different  cases,  reducing  the  ratio  of  the  mean  time  for 
Red  to  break  lock  to  the  mean  time  for  Blue  to  break  lock  from  a  value  of  10  to  a 
value  of  1  decreases  the  equivalent  kill-rate  ratio  by  a  factor  of  about  3. 

This  fact  suggests  that  a  crude  but  useful  method  of  accounting  for  deception  EA 
by  Red  aircraft  in  present  SLAACM  versions  is  to  allow  the  user  to  reduce  kill- 
rate  ratios  by  a  factor  varying  from  1  (no  change)  to  3  (significant  deception  EA). 

A  very  much  better  method  would,  of  course,  be  to  obtain  simulation  results 
(TAC  BRAWLER  or  EADSIM,  for  example)  or  results  from  mock  combat,  with 
sufficient  information  to  infer  effects  on  kill-rate  ratios  directly.  If,  for  example, 
complete  outcome  distributions  for  4  vs.  8  engagements  (sample  probabilities  of  i 
Blue  and  j  Red  surviving  the  engagement)  were  available,  then  the  five  ratios 
kbl/krl,  kb2/krl,  kbx/krl,  kr2/krl,  and  krx/krl  that  characterize  two-phase  kill 
engagements  with  break-lock  could  be  inferred  directly  for  cases  with  Red  EA 
and  without  it.  These  data  are  not  available,  however,  and  until  they  are,  the 
method  proposed  here  may  give  useful  indications. 
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Appendix  F 

Mathematical  Analysis:  Phased  Attack 
and  Defense 


In  air  defense,  it  is  sometimes  necessary  to  provide  continuous  combat  air  patrol 
coverage  of  the  battle  space.  The  current  SLAACM  includes  the  option  for  Blue 
defense  aircraft  to  be  time  phased  over  the  battle  space  in  order  to  provide  con¬ 
tinuous  defense  CAP.  The  inclusion  of  a  phased  Blue  defense  led  to  discussions 
with  TACAIR  about  whether  Red  could  phase  its  attack  to  take  advantage  of 
Blue’s  phasing.  In  response  to  these  discussions,  we  addressed  the  problem  in  two 
substantially  different  ways,  and  include  both  in  the  sections  that  follow.  Ap¬ 
proach  1  demonstrates  a  method  for  generating  probabilistic  results  for  a  well- 
defined,  specific  scenario.  Note  that  in  this  analysis  the  postulated  parameters  are 
selected  solely  to  demonstrate  the  approach.  While  these  parameters  are  generally 
reasonable,  they  do  not  represent  any  known  operational  scenario.  Approach  2  is  a 
generalized  combinatorial  analysis  that  explores  the  probability  space  for  random 
selections  of  phasing  strategies  for  both  Blue  and  Red. 

Approach  1 


We  are  interested  in  the  way  CAP  attrition  is  treated  in  SLAACM.  Presently, 
SLAACM  assumes  that  all  the  packages  Red  sends  in  a  “day”  are  simultaneously 
observable  by  Blue’s  battle  management,  which  may  include  target  identification 
onboard  certain  aircraft.  This  allows  certain  Blue  aircraft  to  selectively  engage 
high-threat  packages.  This  assumption  of  high-threat  Red  attacks  seems  plausible 
when  saturating  Blue’s  defenses  is  an  attractive  option  for  Red. 

But  for  a  scenario  in  which  Blue’s  target  identification/battle  management,  and 
the  two  sides’  available  forces,  would  allow  Red  to  saturate  Blue’s  defense  with 
relatively  low-value  attack  packages  while  Blue  remained  unaware  of  upcoming 
attacks  by  high-value  packages,  the  SLAACM  assumptions  are  optimistic  for 
Blue.  This  is  particularly  true  when  target  identification  is  onboard  aircraft  that  fly 
combat  air  patrols  from  bases  far  from  the  theater. 

Full  exploration  of  the  options  available  to  each  side  by  time  phasing  of  attacks 
and  defenses  will  involve  considerations  of  target  and  basing  geographies  and 
may  well  require  case-by-case  analysis.  Here  we  consider  some  specific,  simple 
examples  to  illustrate  the  potential  significance  of  these  effects. 

Suppose  that  Blue  has  a  fleet  of  192  aircraft,  which  must  accomplish  a  2-hour 
flight  to  the  theater.  Suppose  also  that  their  aircraft  can  remain  on  station  for  1 
hour,  that  they  carry  six  missiles,  that  they  will  keep  one  missile  for  self- 
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protection  while  returning  to  base,  and  that  their  turnaround  time  is  1  hour.  This 
gives  the  Blue  aircraft  1  hour  out  of  six  on  station,  so  that  six  flights  must  be 
available  to  keep  one  flight  continuously  on  station.  Suppose  that  Blue  aircraft 
always  fight  in  four-ship  packages. 

Suppose  that  Red  attacks  always  come  in  12-ship  packages,  made  up  of  advanced 
escorts,  close  escorts,  and  bombers,  and  that  the  Blue  aircraft  must  defeat  the  es¬ 
corts  in  order  to  engage  the  bombers.  Then  each  4-ship  flight,  having  20  usable 
missiles,  can  defeat  only  one  attacking  package  before  rearming. 

Now,  Blue’s  192  aircraft  provide  eight,  six-flight  groups,  so  that  Blue  can  keep 
eight  flights  on  station  continuously.  In  this  simple  example,  Blue  controls  only 
the  time  phasing  of  the  eight  groups’  CAPs. 

Figure  F-l  shows  an  example  of  the  build-up  of  Blue  forces  when  the  CAPs  relief 
times  are  evenly  spaced. 

Figure  F-1.  Arrival  Process  for  Evenly  Spaced  CAPs 
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By  sending  eight  attacking  packages,  Red  can  at  any  time  “reset”  Blue’s  strength 
to  the  start  of  the  buildup,  because  defeating  eight  packages  exhausts  the  missiles 
available  to  the  eight  Blue  flights  on  station. 1  In  this  way  the  Reds  can  give  them¬ 
selves  intervals  of  time  in  which  there  are  fewer  than  eight  Blue  flights. 

Let  us  look  first  at  a  bad  Blue  option  for  the  CAPs’  relief  times.  Suppose  that, 
rather  than  the  uniform  spacing  of  Figure  F-l,  the  eight  groups’  CAPs  are  exactly 
in  phase.  That  is,  the  flights  of  all  eight  groups  relieve  their  predecessors  at  the 
same  time.  Then,  by  sending  eight  low-value  packages,  Red  can  exhaust  the  mis¬ 
siles  available  to  the  Blue  force,  and  there  will  be  no  more  Blue  aircraft  for 
roughly  1  hour.  This  gives  Red  considerable  scope  to  dispatch  higher-value  pack¬ 
ages  unopposed. 


1  Blue  flights  might  regroup  so  that  the  remaining  missiles  available  to  two  or  three  flights 
could  be  used  against  another  Red  package,  but  this  appears  to  call  for  complicated  reorganization 
and  to  give  only  a  second-order  effect  anyway. 
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Presumably,  Blue  will  not  arrange  the  groups’  CAPs  in  this  way.  To  consider 
other  Blue  options  and  Red  tactics  in  more  detail,  we  must  be  more  specific  about 
the  time  required  for  Red  to  carry  out  missions  and  about  the  area  in  which  Red 
packages  are  vulnerable  to  Blue’s  flights. 

Let  us  say,  then,  that  Red’s  packages  are  vulnerable  for  10  minutes  ingress,  5 
minutes  delivering  bombs,  and  10  minutes  egress.  To  focus  strictly  on  time  phas¬ 
ing,  suppose  that  the  Blue  defenders  are  invincible,  so  that  any  Blue  flight  in 
Red’s  vulnerable  zone,  with  at  least  as  many  Red  packages  as  Blue  flights  are  also 
in  the  zone,  destroys  a  Red  package. 

For  further  simplification,  let  us  suppose  that  Red  has  eight  low-value  packages 
and  eight  high-value  packages  to  dispatch.  For  specificity,  suppose  that  bombers 
in  the  high-value  packages  carry  three  times  the  weight  of  bombs  as  those  in  the 
low-value  packages. 

Let  us  also  assume  that  all  eight  Blue  CAPs  are  in  place,  and  that  the  Blue  aircraft 
unfailingly  attack  high-value  packages  in  preference  to  low-value  packages. 

Now,  if  Red  dispatches  all  16  packages  at  once,  the  eight  defending  flights  will 
eliminate  all  eight  high-value  packages,  and  Red  will  deliver  eight  low-value 
bomber  loads  (we’ll  call  this  eight  “units”)  of  bombs.  No  aircraft  in  the  low-value 
packages  will  be  lost. 

For  completeness,  let  us  treat  the  above  “bad”  Blue  option  in  detail.  If  all  eight 
CAPs  were  in  phase,  it  seems  incredible  that  Red  would  not  know  the  times  at 
which  the  CAPs  were  relieved.  Then  Red  can  send  eight  low-value  packages  to  ar¬ 
rive  just  after  a  relief  time.  All  eight  will  be  destroyed.  But  the  missiles  of  all  eight 
CAPs  will  be  exhausted  too.  If  the  eight  high-value  packages  were  undetectable 
while,  say,  20  minutes  behind  the  low-value  packages,  they  could  ingress,  bomb, 
and  be  out  of  the  vulnerable  area  15  minutes  before  the  relieving  CAPs  arrived. 

Blue  would,  however,  want  to  arrange  the  CAPs’  phases  to  avoid  such  an  out¬ 
come.  One  option  would  be  to  space  the  CAPs’  relief  points  evenly  as  shown  in 
Figure  F-l. 

If  Red  sends  the  eight  low-value  packages  to  arrive  just  at  one  CAP’s  relief  time, 
and  delays  the  eight  high-value  packages  for,  say,  20  minutes  to  make  sure  that 
Blue  does  not  know  they  are  coming,  the  high-value  packages  will  face  only  two 
defending  flights,  with  a  third  defending  flight  joining  after  2.5  minutes.  One 
more  Blue  flight  will  arrive  just  as  the  five  Red  packages  start  bombing;  it  will 
destroy  a  Red  package,  and  four  high-value  Red  packages  will  drop  their  bombs. 

Two  defending  flights  will  arrive  before  the  four  packages  that  survive  ingress 
and  bombing  leave  the  vulnerable  area  (the  second  arrives  just  as  they  leave,  but  it 
seems  reasonable  to  give  Blue  the  benefit  of  a  tie),  so  that  the  outcome  of  the  at¬ 
tack  is  six  high-value  and  eight  low-value  packages  destroyed,  with  four  high- 
value  bomber  loads  of  bombs  delivered. 
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Thus,  spacing  the  CAPs’  relief  times  evenly  allows  Blue  to  destroy  six  more  high- 
value  packages,  and  allows  four  fewer  bomber  flights  to  "leak"  through  Blue's 
defenses,  than  would  have  been  the  case  if  the  relief  times  coincided.  Red,  on  the 
other  hand,  gets  four  high-value  bomber  flights  through,  delivering  50  percent 
greater  weight  of  bombs,  and  loses  two  fewer  high-value  packages,  than  he  would 
have  done  by  sending  all  the  packages  simultaneously  (as  SLAACM  now  as¬ 
sumes  he  would  do)  against  the  evenly  spaced  CAPs. 

As  a  final  example,  we  suppose  that  the  eight  CAPs’  relief  times  are  randomly 
shifted  in  time  and  that  Red  sends  eight  low-value  packages  at  a  random  time.  Let 
us  generalize  the  discussion,  assuming  that  the  defending  flights’  time-on-station 
is  S  minutes  and  that  Red  follows  up  eight  high-value  packages  D  minutes  later. 
Further  suppose  that  ingress  and  egress  require  I  and  E  minutes,  respectively,  and 
that  bombing  takes  B  minutes. 

Then,  the  number  j  of  Blue  defending  flights  arriving  while  the  high-value  pack¬ 
ages  ingress  and  drop  bombs  has  the  binomial  distribution  B(j,  8,  [D  +  I  +  B]/S). 
The  number  p  of  flights  that  do  not  arrive  during  that  period  is,  of  course,  distrib¬ 
uted  as  B(p,  8,  [S  -  D  -  I  -  B]/S). 

When  k  defending  flights  arrive  during  ingress  and  bombing,  then  8  -  k  flights 
arrive  in  the  interval  S  -  I  -  B,  and  their  arrival  times  are  uniformly  distributed 
over  that  interval.  Thus,  the  number  m  of  defending  flights  arriving  during  egress 
has  the  binomial  distribution  B(m,  8  -  k,  E/[S  -  I  -  B]),  for  0  <  m  <  8  —  k. 

These  results  give  the  distribution  of  the  number  p  of  bomber-loads  of  bombs 
dropped  as  B(p,  8,  [S  -  D  -  I  -  B]/S),  and  the  distribution  P(n)  of  the  total  num¬ 
ber  n  of  attack  packages  destroyed  as 

p(n)  =  ZB(n-j/8,pi)B(j/8-n  +  j/p2), 

j=o 

where  pi  =  (D  +  I  +  B)/S,  and  p2  =  E/(S  -  D  -  I  -  B).  With  this  information,  one 
can  plot  the  statistics  of  bomb  units  dropped  and  Red  packages  destroyed  (each 
high-value  package  delivers  three  bomb  units).  Figures  F-2  and  F-3  show  the  re¬ 
sults  for  S  =  60,  D  =  20, 1  =  E  =  10,  and  B  =  5. 


F-4 


Appendix  F:  Mathematical  Analysis:  Phased  Attack  and  Defense 


Figure  F-2.  Probability  Distribution  of  Bomb  Units  Dropped 


Figure  F-3.  Probability  Distribution  of  Red  Packages  Lost 


By  sacrificing  eight  low-value  packages,  Red  has  72  percent  confidence  of  deliv¬ 
ering  more  bomb  units  than  with  the  simultaneous  assault  now  assumed  in 
SLAACM,  and  a  90  percent  confidence  of  losing  fewer  high-value  packages 
(these  marginal  statistics  may  not,  of  course,  occur  together). 

Approach  2 


As  discussed  above,  we  are  interested  in  determining  if  Red  can  systematically 
phase  his  attack  to  maximize  the  success  of  his  high-value  attack  packages.  We 
assume  here,  as  above,  that  both  the  Red  attack  and  the  Blue  defense  can  be  time 
phased  within  the  engagement  period,  and  each  four-ship  Blue  flight  can  engage 
and  defeat  only  one  12-ship  Red  package.  We  further  assume  that  the  Red  attack 
is  numerous,  and  Red  is  actually  time  phasing  his  high-value  packages  within  a 
large  number  of  low-value  packages,  such  that  Blue  defense  flights  will  engage 
the  high-value  Red  packages  available  immediately  on  arrival  and  will  engage  the 
lower  value  packages  if  no  high-value  packages  are  available.  All  Blue  defense 
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flights  will  be  absorbed  by  low-value  Red  packages  whenever  high-value  pack¬ 
ages  are  not  available. 

Consider  the  example  in  which  eight  Blue  flights  can  arrive  during  a  single  en¬ 
gagement.  We  assume  the  engagement  is  divided  into  equal  periods  correspond¬ 
ing  to  the  number  of  Blue  flights,  eight  in  this  example.  Blue  can  distribute  his 
flights,  and  Red  can  distribute  his  high-value  packages  arbitrarily  among  the  eight 
periods.  Because  each  Blue  flight  engages  only  one  Red  package,  Blue  can  kill 
(and  Red  can  lose)  a  maximum  of  eight  high-value  Red  packages. 

We  wish  to  determine  if  there  is  a  way  to  phase  the  Blue  arrivals  to  ensure  maxi¬ 
mum  kills  of  high-value  Red  packages.  Our  assumptions  are  as  follows: 

♦  Blue  is  invincible  and  kills  Red  packages  with  100  percent  probability. 

♦  Blue  flights  can  preferentially  detect  and  attack  high-value  Red  packages. 

♦  All  Blue  flights  will  be  committed  (or  consumed)  in  each  period  either  by 
high-  or  low-value  Red  packages. 

♦  Red  must  commit  all  of  his  high-value  packages  during  the  engagement. 


Analysis 


By  inspection  we  can  deduce  that  Blue  can  guarantee  at  least  one  high-value  pack¬ 
age  kill  by  distributing  his  arrivals  equally.  Beyond  that,  the  results  are  not  obvious. 

To  address  the  Blue  and  Red  options,  we  need  to  determine  the  number  of  ways 
Blue  and  Red  can  distribute  eight  flights  over  eight  time  slots.  This  is  a  multinomial, 
combination-permutation  in  which  order  is  important  among  different  quantities, 
e.g.,  (1,  2)  ^  (2,1),  but  not  important  among  multiple  occurrences  of  the  same  quan¬ 
tities,  e.g.,  (2i,  22)  =  (22,  2i).  We  first  demonstrate  a  calculation  method  using  four 
balls  in  four  bins,  and  later  apply  the  method  to  eight  flights  in  eight  time  slots. 

Four  Balls  Distributed  Among  Four  Bins 

First,  we  note  that  there  are  4!  ways  to  distribute  four  distinguishable  items,  such  as 
4  colored  balls  among  4  bins,  1  to  a  bin.  In  our  problem,  however,  the  0,  1,2,  3,  and 
4  counts  of  balls  in  a  bin  are  distinguishable,  but  the  balls  themselves  and  multiple 
occurrences  of  the  counts  are  indistinguishable.  The  choice  of  placing  0,  1,2,  3,  or  4 
balls  in  a  given  bin  restricts  the  options  for  the  remaining  bins.  For  example,  if  4 
balls  are  placed  in  any  1  bin,  the  rest  of  the  bins  must  hold  0.  (This  is  same  as  having 
1  red  ball  and  3  blue  balls.)  The  order  of  the  bin  containing  the  4  balls  is  significant, 
so  there  are  4  ways  to  distribute  the  4-ball  set  among  the  4  bins.  There  is  only  1 ! 
way  to  arrange  the  single  bin  containing  4  indistinguishable  balls.  If  the  0s  were  dis¬ 
tinguishable,  there  would  be  3 !  ways  to  distribute  them  among  3  bins,  but  we  do  not 
distinguish  differences  in  order  among  the  bins  containing  0s.  When  we  reduce  the 
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4!  maximum  options  to  account  for  the  indistinguishable  bins  containing  the  same 
counts  of  balls,  i.e.,  the  three  Os,  we  get  4!/(3!  *  1 !)  =  4  unique  distributions.  Ta¬ 
ble  F-l  confirms  this  by  showing  the  complete  set  of  4  unique  states. 

Table  F-1.  State  Count  for  Four 
Balls  in  a  Given  Bin 


Distribution 

count 

Bins 

1 

2 

3 

4 

1 

4 

0 

0 

0 

2 

0 

4 

0 

0 

3 

0 

0 

4 

0 

4 

0 

0 

0 

4 

To  determine  the  total  number  of  unique  distributions  for  4  balls  in  4  bins,  we 
need  to  identify  all  the  positional  distributions  of  unique  ball  counts.  This  we  do 
by  hand.  Table  F-2  shows  the  unique  distributions  for  4  balls  in  4  bins  and  the 
corresponding  calculations  of  state  counts. 

Table  F-2.  Unique  Distributions  and  State  Counts  for  Four  Balls  in  Four  Bins 


Bins 

Distribution 

counts 

Distribution 

counts 

1 

2 

3 

4 

Unique 

distribu- 

tions 

1 

1 

1 

1 

4!/4! 

1 

2 

1 

1 

0 

4!/(1!*2!*1!) 

12 

2 

2 

0 

0 

4!/(2!*2!) 

6 

3 

1 

0 

0 

4!/(1!*1!*2!) 

12 

4 

0 

0 

0 

4!/(1  !*3!) 

4 

Total 

35 

As  a  check  on  the  method,  Table  F-3  shows  the  detailed  state  enumeration  for  4 
balls  in  4  bins. 
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Eight  Flights  in  Eight  T ime  Slots 

Now  we  consider  eight  flights  distributed  among  eight  time  slots.  The  unique  dis¬ 
tributions  and  distribution  counts  are  shown  in  Table  F-4. 

Table  F-4.  Unique  Distributions  and  State  Counts 
for  Eight  Flights  in  Eight  Time  Slots 


Time  slot 

Distribution 

Distribution 

1 

2 

3 

4 

5 

6 

7 

8 

counts 

counts 

1 

1 

1 

1 

1 

1 

1 

1 

8!/8! 

1 

2 

1 

1 

1 

1 

1 

1 

0 

8!/(6!  111!) 

56 

2 

2 

1 

1 

1 

1 

0 

0 

8!/(4!  2!  2!) 

420 

2 

2 

2 

1 

1 

0 

0 

0 

8!/(3!  3!  2!) 

560 

2 

2 

2 

2 

0 

0 

0 

0 

8!/(4!  4!) 

70 

3 

1 

1 

1 

1 

1 

0 

0 

8!/(5!  2!  1!) 

168 

3 

2 

1 

1 

1 

0 

0 

0 

8!/(3!  3!  1!  1!) 

1,120 

3 

2 

2 

1 

0 

0 

0 

0 

8!/(4!  2!  1!  1!) 

840 

3 

3 

1 

1 

0 

0 

0 

0 

8!/(4!  2!  2!) 

420 

Unique 

distribu¬ 

tions 

3 

3 

2 

0 

0 

0 

0 

0 

8!/(5!  2!  1!) 

168 

4 

1 

1 

1 

1 

0 

0 

0 

8!/(4!  3!  1!) 

280 

4 

2 

1 

1 

0 

0 

0 

0 

8!/(4!  2!  1!  1!) 

840 

4 

2 

2 

0 

0 

0 

0 

0 

8!/(5!  2!  1!) 

168 

4 

3 

1 

0 

0 

0 

0 

0 

8!/(5!  1!  1!  1!) 

336 

4 

4 

0 

0 

0 

0 

0 

0 

8!/(6!  2!) 

28 

5 

1 

1 

1 

0 

0 

0 

0 

8!/(4!  3!  1!) 

280 

5 

2 

1 

0 

0 

0 

0 

0 

8!/(5!  1!  1!  1!) 

336 

5 

3 

0 

0 

0 

0 

0 

0 

8!/(6!  111!) 

56 

6 

1 

1 

0 

0 

0 

0 

0 

8!/(5!  2!  1!) 

168 

6 

2 

0 

0 

0 

0 

0 

0 

8!/(6!  111!) 

56 

7 

1 

0 

0 

0 

0 

0 

0 

8!/(6!  111!) 

56 

8 

0 

0 

0 

0 

0 

0 

0 

8!/(7!  1!) 

8 

Total 

6,435 

From  Table  F-4  we  see  that  there  are  6,435  ways  to  distribute  eight  Blue  flights 
(and  eight  Red  packages)  among  eight  time  slots.  Now  we  want  to  see  the  poten¬ 
tial  impact  of  this  on  Blue  payoff. 

To  find  the  potential  payoff  for  Blue,  we  want  to  find  the  probabilities  of  Blue 
(and  Red)  experiencing  0  through  eight  flights  in  a  given  slot.  Table  F-5  shows 
the  flight  count  occurrences  for  each  unique  distribution. 
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Table  F-5.  Flight  Count  Occurrences  per  Unique  Distribution 


Unique  state  distributions, 
by  time  slot 

Flight  count  occurrences  per 
unique  distribution,  by  time  slot 

Distribution 

' 

3 

4 

3 

6 

3 

8 

0 

1 

2 

3 

4 

5 

6 

7 

8 

counts 

D 

D 

D 

D 

D 

D 

D 

D 

8 

1 

D 

D 

D 

D 

D 

D 

□ 

D 

□ 

D 

56 

0 

420 

0 

□ 

B 

B 

560 

□ 

4 

4 

70 

D 

D 

□ 

□ 

n 

□ 

1 

168 

0 

0 

0 

D 

D 

1120 

□ 

□ 

□ 

B 

D 

840 

0 

B 

2 

420 

Unique 

distribu¬ 

tions 

□ 

□ 

□ 

0 

D 

B 

168 

i 

0 

0 

0 

□ 

1 

280 

□ 

□ 

□ 

0 

□ 

D 

1 

840 

0 

2 

1 

168 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

D 

D 

D 

336 

0 

0 

0 

0 

0 

6 

2 

28 

D 

□ 

□ 

□ 

□ 

□ 

B 

1 

280 

□ 

D 

D 

1 

336 

□ 

□ 

□ 

□ 

□ 

6 

i 

1 

56 

0 

0 

0 

0 

0 

i 

m 

1 

168 

□ 

H 

□ 

□ 

□ 

6 

i 

1 

56 

□ 

D 

1 

56 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

□ 

7 

1 

8 

6,435 

Dividing  the  flight  count  occurrences  in  a  given  row  of  Table  F-5  by  8  gives  the 
conditional  probability  of  the  count  occurrences  given  the  distribution  corre¬ 
sponding  to  the  row.  Dividing  the  distribution  counts  by  the  total  count  gives  the 
probability  of  each  distribution.  Multiplying  the  distribution  probabilities  by  the 
conditional  count  occurrence  probabilities  generates  the  probabilities  for  count 
occurrences  shown  in  Table  F-6. 

Table  F-6.  Flight  Occurrences  Probabilities  for  Eight  Flights  in  Eight  Time  Slots 


Distribution 

probabilities 

Flight  count  probabilities 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0.00016 

0.0000 

0.0002 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.00870 

0.0011 

0.0065 

0.0011 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.06527 

0.0163 

0.0326 

0.0163 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.08702 

0.0326 

0.0218 

0.0326 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.01088 

0.0054 

0.0000 

0.0054 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.02611 

0.0065 

0.0163 

0.0000 

0.0033 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 
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Table  F-6.  Flight  Occurrences  Probabilities  for  Eight  Flights  in  Eight  Time  Slots 


Distribution 

probabilities 

Flight  count  probabilities 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0.17405 

0.0653 

0.0653 

0.0218 

0.0218 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.13054 

0.0653 

0.0163 

0.0326 

0.0163 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.06527 

0.0326 

0.0163 

0.0000 

0.0163 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.02611 

0.0163 

0.0000 

0.0033 

0.0065 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.04351 

0.0163 

0.0218 

0.0000 

0.0000 

0.0054 

0.0000 

0.0000 

0.0000 

0.0000 

0.13054 

0.0653 

0.0326 

0.0163 

0.0000 

0.0163 

0.0000 

0.0000 

0.0000 

0.0000 

0.02611 

0.0163 

0.0000 

0.0065 

0.0000 

0.0033 

0.0000 

0.0000 

0.0000 

0.0000 

0.05221 

0.0326 

0.0065 

0.0000 

0.0065 

0.0065 

0.0000 

0.0000 

0.0000 

0.0000 

0.00435 

0.0033 

0.0000 

0.0000 

0.0000 

0.0011 

0.0000 

0.0000 

0.0000 

0.0000 

0.04351 

0.0218 

0.0163 

0.0000 

0.0000 

0.0000 

0.0054 

0.0000 

0.0000 

0.0000 

0.05221 

0.0326 

0.0065 

0.0065 

0.0000 

0.0000 

0.0065 

0.0000 

0.0000 

0.0000 

0.00870 

0.0065 

0.0000 

0.0000 

0.0011 

0.0000 

0.0011 

0.0000 

0.0000 

0.0000 

0.02611 

0.0163 

0.0065 

0.0000 

0.0000 

0.0000 

0.0000 

0.0033 

0.0000 

0.0000 

0.00870 

0.0065 

0.0000 

0.0011 

0.0000 

0.0000 

0.0000 

0.0011 

0.0000 

0.0000 

0.00870 

0.0065 

0.0011 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0011 

0.0000 

0.00124 

0.0011 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0002 

Totals 

0.4667 

0.2667 

0.1436 

0.0718 

0.0326 

0.0131 

0.0044 

0.0011 

0.0002 

Table  F-7  repeats  the  occurrence  probabilities  from  Table  F-6  and  adds  the  Blue 
payoff  probabilities,  which  are  calculated  by  multiplying  the  probabilities  by  their 
corresponding  flight  counts.  The  expected  payoff  shown  in  Table  F-7  is  1  com¬ 
pared  to  the  maximum  payoff  of  8. 

Table  F-7.  Probabilities  and  Payoff  for  Eight  Flights  in  Eight  Time  Slots 


Probability  per  state 

0 

1 

2 

3 

4 

5 

6 

7 

8 

Total 

0.4667 

0.2667 

0.1436 

0.0718 

0.0326 

0.0131 

0.0044 

0.0011 

0.0002 

1.000 

Blue  payoff/Red  loss  probability 

0 

1 

2 

3 

4 

5 

6 

7 

8 

Total 

0 

0.2667 

0.2872 

0.2154 

0.1305 

0.0653 

0.0261 

0.0076 

0.0012 

1.000 

Look-Ahead  Target  Identification 

The  case  in  which  Blue  can  see  beyond  the  current  period,  and  thus  avoid  wasting 
flights  on  low-value  packages,  can  be  modeled  simply  by  reducing  the  number  of 
time  slots.  Using  the  same  approach  followed  above,  we  find  that  the  expected 
payoff  for  eight  flights  in  four  time  slots  is  2,  and  the  expected  value  for  eight 
flights  in  two  time  slots  is  4.  The  simplicity  of  these  numbers  suggests  that  there 
may  be  a  more  fundamental  way  to  derive  them  than  we  have  applied.  That  said, 
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the  results  demonstrate  the  value  of  long-range  identification  to  Blue  and  the 
value  of  preventing  such  identification  to  Red. 

Approach  2  Summary 

We  note  that  these  results  represent  the  case  in  which  Blue  and  Red  each  ran¬ 
domly  select  one  of  an  exhaustive  set  of  strategies.  It  does  not  represent  the  case 
in  which  Blues  and  Reds  arrive  randomly.  The  case  of  random  phasing  of  Blue 
arrivals,  and  random  arrival  of  the  Red  attack,  is  considered  in  Approach  1 . 

At  this  time,  we  do  not  know  of  any  gaming  strategy  that  can  reliably  improve  the 
results  for  either  Blue  or  Red. 


F-12 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (MM-YYYY)  2.  REPORT  TYPE 

05-2007  Final 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 

Stochastic  Lanchester  Air-to-Air  Campaign  Model  GS-00F-0026M 


3.  DATES  COVERED  (From  -  To) 


Methods  Used  to  Generate  Model  Outputs 
and  a  Users  Guide:  2007 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Hemm,  Robert  V. ;  Author 

Lee,  David  A.;  Author 
Eckhause,  Jeremy  M.;  Author 
Dukovich,  John  A. ;  Author 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 

PA603 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

LMI 

2000  Corporate  Ridge 
McLean,  VA  22102-7805 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

LMI-PA603T2 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

OSD/Program  Analysis  and  Evaluation  -  TACAIR 
Mr.  Frank  Lewis 
Pentagon,  Room  2C281 
Washington,  DC  20301-1800 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

A  Approved  for  public  release;  distribution  is  unlimited. 


14.  ABSTRACT 

The  report  describes  latest  version  of  the  Stochastic  Lanchester  Air-to-Air  Campaign  Model 
(SLAACM) .  The  report  has  two  major  sections  and  appendices.  The  first  section  includes 
several  chapters  that  describe  the  mathematics  and  analysis  methods  used  for  engagement 
modeling,  bomb  damage  calculations,  campaign  planning,  and  integer  program  optimization.  The 
second  section  is  a  user's  guide  for  the  latest  SLAACM  release.  The  appendices  describe  the 
unique  and  efficient  engagement  algorithm  used  in  SLAACM,  supporting  Markov  analysis  tools, 
mathematics  to  support  additional  combat  scenarios,  and  results  of  supplemental  analyses. 


15.  SUBJECT  TERMS 

Multiple  air  combat,  campaign  models,  Lanchester  models 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Nancy  E .  Handy 

a.  REPORT 

UNCLASSIFIED 

b.  ABSTRACT 

UNCLASSIFIED 

c.  THIS  PAGE 

UNCLASSIFIED 

Unclassified 

Unlimited 

58 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

703-917-7249 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


